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A STUDY  OF  THE  ELECTROMAGNETIC  FIELDS  DISTRIBUTION 
INSIDE  BUILDINGS  WITH  APERTURES  EXCITED  BY  AN  EXTERNAL  SOURCE 

M.  T.  Ma  and  M.  G.  Arthur 


Two  special  cases  of  the  penetration  of  electromagnetic  fields 
into  a cavity,  building  or  box  are  formulated  and  analyzed.  One 
is  to  consider  the  case  of  a lossy  cavity  with  small  apertures  in 
free  space,  based  on  an  application  of  the  equivalence  principle 
and  the  use  of  a generalized  network  formulation.  It  is  found 
that  the  field  strength  at  the  aperture  center  is  approximately 
inversely  proportional  to  the  square-root  of  the  conductivity  of 
the  cavity  walls  and  that  high  field  levels  can  exist  inside  the 
cavity  under  certain  physical  conditions.  The  second  case  is  to 
treat  the  problem  of  large  buildings  with  large  apertures  on  a 
practical  lossy  ground  by  a combination  of  theoretical  approach 
and  limited  measurement  data.  Field  levels  inside  the  building 
for  this  latter  case  depend  on  the  transmitter  power,  the 
transmitter-to-building  distance,  the  ground  conductivity,  and  the 
measurement  antenna  height  relative  to  that  of  the  transmitter. 

Keywords:  aperture;  cavity;  equivalence  principle;  field 

distribution;  slot. 


1.  INTRODUCTION 

The  capability  of  electromagnetic  waves  to  penetrate  into  a building  or  a 
box  which  houses  important,  sensitive  devices  and  equipment  has  always  been  of 
concern  to  industries  and  military  installations.  Generally  speaking,  it  is 
essential  to  know  the  level  of  the  field  strength  and  its  component  amplitudes 
and  phases  inside  the  building  before  an  assessment  of  hazardous  effects  on 
equipment  performance  and  human  health  can  be  made.  Genuine  efforts  have  been 
made  to  make  on-site  spot  measurements  inside  an  empty  building  for  a known 
transmitter  at  a specified  distance  away  from  the  building.  It  is  still 
impossible  from  these  spot  measurements  to  make  any  meaningful  predictions  of 
the  highest  field  level  and  its  polarization  inside  the  building,  which  are 
the  ultimate  goals  of  this  study.  Thus,  to  achieve  these  goals,  systematic 
measurements  of  component  field  amplitudes  and  phases  along  three  orthogonal 
directions  at  an  interval  distance  of  no  more  than  one  quarter  wavelength 
would  be  required.  These  systematic  measurements  are,  of  course,  very  time- 
consuming.  For  a transmitter  operated  at  150  MHz  (A - m 1 , which  is  typical  or 
a land  mobile  radio,  and  a practical  storage  building  ot  the  si_e  of  \ \ 

X iox , a huge  number  of  measurement  data  in  the  order  of  one  third  of  a 
million  data  points  have  to  be  taken  and  analyzed.  Furthermore,  these 
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measurement  data  are  only  applicable  to  the  specific  building  and  depend  on 
construction  material,  building  size  and  shape,  the  aperture  (doors  or  windows) 
size  and  shape,  the  ground  and  surrounding  conditions,  the  operating  frequency  of 
the  transmitter,  and  the  transmi tter-to-bui 1 di ng  distance  under  study. 

To  circumvent  this  measurement  impracti cal i ty , a theoretical  study  of  the 
generalized  three-dimensional  problem  was  deemed  desirable.  However,  in  view 
of  the  complexity  of  the  problem  involved  and  the  enormous  demand  on  the 
computing  facility  for  a practical  size  of  building,  only  special  cases  for  a 
CW  excitation  subject  to  certain  physical  constraints  are  considered  in  this 
report. 

In  Section  2,  the  problem  of  a three-dimensional  closed  conducting 
rectangular  box  in  free  space  illuminated  by  a linearly  polarized,  normally 
incident,  plane  wave  is  considered  with  two  different  approaches.  One  is  by 
the  magnetic  field  integral  equation  with  appropriate  boundary  conditions  [1- 
4],  and  the  other  is  through  the  concept  of  a generalized  impedance  matrix. 

Since  there  is  no  opening  on  the  box  surface  in  this  case,  no  field  can 
penetrate  into  the  box.  Numerical  results  for  the  induced  surface  current 
densities,  on  a limited  basis  are  given  therein,  based  on  which  the  field 
scattered  by  the  box  may  be  computed. 

The  problem  is  then  extended  to  the  case  of  lossy  rectangular  buildings  or 
boxes  with  small  slot  apertures  in  free  space,  which  is  presented  in  Section 
3.  Field  strengths  inside  the  box  are  shown  to  be  dependent  on  the  aperture 
and  box  sizes  relative  to  the  operating  wavelength  and  the  material  of  the  box 
wal 1 s . 

It  is  found  that  the  field  strength  at  the  aperture  center  is 
approximately  inversely  proportional  to  the  square-root  of  the  conductivity  of 
the  wall  material  and  that  high  levels  of  the  field  strength  can  exist  inside 
the  box  at  natural  resonances  and  under  a set  of  special  physical  conditions. 

In  Section  4,  large  buildings  and  large  apertures  located  on  a practical 
lossy  ground  are  analyzed  with  a combination  of  theoretical  consideration  and 
limited  set  of  measurement  data.  Requirements  on  the  measurement  equipment  and 
computing  facilities  for  this  case  are  outlined  in  Section  5 with  the 
measurement  and  computing  procedures  presented  in  Section  6.  Numerical  results 
of  the  approximate  induced  current  densities  and  the  approximate  field 
distribution  inside  the  building  for  a known  transmitter  frequency,  power,  and 
the  transmi tter-to-bui 1 di nq  distance  are  given  in  Section  7.  Physical 


2 


interpretations  of  these  numerical  results  and  the  limiting  conditions  under 
which  the  results  have  meanings  are  discussed  in  the  concluding  section.  A 
relevant  FORTRAN  program  used  to  obtain  the  numerical  results  is  listed  in 
Appendix  A.  A simplified  version  using  the  HPL*  suitable  to  Model  HP  9825* 
capable  of  making  a quick  on-site  check  with  a much  reduced  matrix  size  is 
presented  in  Appendix  B. 

2.  A THREE-DIMENSIONAL  CONDUCTING  RECTANGULAR  BOX 

The  geometry  depicted  in  figure  1 may  be  used  to  represent  a theoretical 
model  for  the  problem  of  the  penetration  of  electromagnetic  waves  into  a 
rectangular  cavity  (building  or  box).  The  source  located  outside  the  cavity 
may  be  a CW  or  pulsed  signal  of  arbitrary  polarization  and  frequency,  incident 
on  the  cavity  with  an  arbitrary  angle.  The  cavity  walls  are  generally  made  of 
non-perfectly  conducting  material,  and  the  cavity  itself  is  situated  on  a 
homogeneous  flat  lossy  earth  represented  by  the  ground  constants  of 
permeability  (y),  permittivity  (e)  and  conductivity  (a). 

The  objective  of  formulating  and  analyzing  this  problem  is  to  determine 
the  spatial  field  distribution  and  locations  of  the  highest  levels  of  field 
strength  inside  the  cavity  as  a function  of  many  parameters  such  as  the  sizes 
of  the  cavity  and  the  aperture,  the  ground  constants,  etc.  Obviously,  this 
type  of  problem  is  too  complicated  to  be  treated  rigorously  using  Maxwell's 
equations  with  appropriate  boundary  conditions.  To  simplify  the  problem,  a 
three-dimensional  perfectly  conducting  rectangular  box  without  any  apertures 
located  in  free  space,  illuminated  by  an  external  source,  as  shown  in  figure  2 
may  be  considered  as  a first  approximate  model  to  determine  the  induced  surface 
current  densities,  from  which  the  scattered  field  outside  the  box  may  then  be 
calculated.  The  numerical  results  so  obtained  can  be  compared  with  the 
existing  measurement  data  on  radar  cross-secti ons  to  assess  the  validity  of  the 
formulation  and  the  accuracy  of  the  numerical  computations. 


♦"Certain  commercial  equipment,  instruments,  or  materials  are  identified  in 
this  paper  in  order  to  adequately  specify  the  experimental  procedure.  Such 
identification  does  not  imply  recommendation  or  endorsement  by  tie  ie’v. 
Bureau  of  Standards,  nor  does  it  imply  that  the  materials  or  equipment 
identified  are  necessarily  the  best  available  for  the  purpose. 
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2.1  Magnetic  Field  Integral  Equation 

Two  approaches  have  been  available  for  deriving  the  necessary  integral 
equations  to  determine  the  current  densities  induced  on  the  surfaces  of  the 
box  in  figure  2.  The  first  is  to  formulate  the  electromagnetic  boundary-val ue 
problem  in  terms  of  the  magnetic  field  integral  equation,  which  takes  the 
following  form  [5]: 

j J (r)  = n X H (r)+T^nX}  j (?')  X V'G(r.r')  ds'  (1) 

L 1 4tT  S 

where  H.  (r)  is  the  known  incident  magnetic  field,  J (r)  is  the  unknown 
tangential  current  density  induced  on  the  six  box  surfaces, 

G (r~r')  = exp  ( - j kR ) /R  (2) 

R = \r  - r'\  = [(x  - x')2  + (y  - y'f  + (z  - z')2]1/2  (3) 

r = xx+yy  + zz  (4) 

r'  = x x'  + y y'  + z z'  (5) 

x,  x-',  y,  y' , z,  z'  = coordinates  on  the  box  surfaces 
x,  y and  z = unit  vectors  in  x,  y,  and  z directions 


V" 


(6) 


n = a unit  vector  in  the  outward  normal  direction  to  each  box 
surface 

S = the  area  of  each  of  the  six  surfaces, 

and  the  horizontal  bar  across  the  integral  sign  in  (1)  denotes  that  only  the 
principal  value  will  be  taken. 


4 


Since  there  are  two  components  of  J (r)  on  each  surface  of  the  box, 
equation  (1)  results  in  12  coupled  scalar  integral  equations  for  the  component 
current  densities.  As  an  example,  in  reference  to  figure  2,  we  have 

H.  (r)  = x H.  (r)  + z H.  (r) , (7) 

1 IX  1Z 

for  the  wave  propagating  in  the  y -direction,  normal  to  Face  I (right).  Tnus, 
for  faces  I and  II  (right  and  left), 

3 (r" ) = x Jx  ( r' ) + z Jz  ( r" ) , (8) 

for  faces  III  and  IV  (top  and  bottom), 

3 {?')  = x J (?' ) + y J (r-),  (9) 

^ y 

for  faces  V and  VI,  (front  and  back) 

j {?' ) = y Jy  (?')  + z Jz  (^)»  (10) 


and 


where 


and 


V"  G (r,  r')  = (x 


3x" 


+ y 


9y^ 


+ z 


3Z 


:)  exp  ( - j kR)/R 


x it  + y it  + Z\i) 
yx  y z 


(11) 


ip  = ( jk  + i)(x  - x^)  exp  (- jkR)/R2  , (12) 

X 

\py  = (jk  +{)(y-yO  exp  ( - j kR ) / R2  , (13) 

ipz  = (jk  +-^)(z-z")  exp  (-jkR)/R2 . ^14^ 
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The  integrand  of  (1)  then  becomes 
for  Faces  I and  II, 


x (-Jz  fy)  + y Ux  Jz  " Jx  ^ + z (Jx  ^y^’ 


(15) 


for  Faces  III  and  IV, 


X ( JyXpz  ) + y (“Jx  + Z (Jx  Vy  “ Jy)> 


(16) 


and  for  Faces  V and  VI 


x (Jy^z  " Jz  ^y)  + y (Jz  ^x)  + z (_^x  Jy)‘ 


(17) 


After  taking  the  principal  value  of  (1),  we  finally  have, 
for  Face  I where  n = -y. 


7 J (r) 

l V 


and 


Jx(r')yls'+  4!/  J (r')lMis'  (18) 

II, III, IV  y III, IV, V, VI  y 


T J (r)  = H.  (r)  - — / J (r')  ij;  ds"  + — J J (r")  ^ ds"  (19) 

2 z lx  4tt  J z y 4tt  j y z 

1 1 , V , VI  y III, IV, V, VI  y 


Equations  similar  to  (18)  and  (19)  can  be  obtained  for  Faces  II  through  VI, 
resulting  in  exactly  12  coupled  scalar  integral  equations.  If  the  method  of 
moments  [6]  is  used  to  solve  for  this  system  of  equations  by  subdividing  each 

O O 

surface  into  n x n patches,  we  will  end  up  1 2n ^ x 12n  simultaneous  linear 
algebraic  equations.  When  the  value  of  n is  sufficiently  large  to  insure 
meaningful  results,  the  inversion  of  a matrix  size  involved  may  be  prohibitive 
even  with  the  availability  of  a large  computer  capacity.  The  results  of  |J  | 
and  | J z | on  Face  I for  a (0.5x)^  cube  illuminated  by  a broadside  incidence 
plane  wave  H^z  with  unity  strength  are  presented  in  figure  3,  which  were 
obtained  with  8x8  patches  [5].  Similar  results  for  other  box  surfaces  are 
also  available  [5].  Based  on  the  results  of  surface  current  densities,  the 
scattered  fields  outside  the  cube  may  then  be  computed.  For  larger  sizes  of 
cube  or  arbitrary  incidence,  the  numerical  procedures  are  more  involved. 
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2.2  Generalized  Impedance  Approach 


An  alternative  approach  to  solving  the  problem  depicted  in  figure  2 is  by 
using  the  concept  of  the  generalized  impedance  matrix.  Basically,  the  electric 
field  scattered  by  the  box  shown  in  figure  2 can  be  determined  by 

E = - jwA  - V<}> , (20) 

where  A and  <j>  are,  respectively,  the  magnetic  vector  potential  and  the  electric 
scalar  potential,  satisfying  the  following  equations: 

A = /s  3 (r)  G (r,  r ) ds'  (21) 

and 

4>  = - 3^4^  Js  v'«  j C?')  G(r,  r' ) ds'  (22) 

where  J and  G have  been  explained  or  defined  in  (1)  and  (2). 

On  the  surface  of  a perfectly  conducting  box,  the  tangential  component 
of  the  total  electric  field  vanishes,  or 


i ,tan 


- - E 


's  ,tan 


(23) 


where  E;  is  the  tangential  component  of  the  electric  field  at  the  box 

I 5 L a H — 

urface  due  to  the  incident  signal,  and  Es  tan  is  the  tangential  component  of 
lat  in  (20). 

Combining  (20)  and  (23),  we  obtain 


E. 
i , 


tan 


[ ^ fs  j (f') G (’r-  "r") ds' 


+ 


l 

4tt  jwe 


J (r) 


. V'  G (r. 


(24) 


= L [J  (r')] 


where  the  linear  operator  L represents  the  necessary  double  integrals. 
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Using  the  standard  method  of  moments  [6],  we  subdivide  the  surface  S into 
patches  and  expand  the  components  of  the  unknown  surface  current  density 
function  J (r')  as  follows: 

n 

J.  (r")  = j-th  component  of  J lr')  = V If  (25) 

j L , n n 

n=l 

where  f is  a set  of  chosen  basis  functions,  and  I ' s are  complex  coefficients 
to  be  determined. 

Substitution  of  (25)  into  (24)  yields 


where 


(26) 


(27) 


' <wi’ 

Lf  i> 

<Wr 

Lf  2> 

...  <wx. 

LV' 

rz  1 = 

L mnJ 

<w2, 

Lfl> 

<w2. 

Lf  2> 

• • • ^ j 

LfN> 

(28) 

, <WN» 

Lfl> 

• • • 

<WN- 

Lf2> 

• • • » 

LfN>. 

E.  .>  ' 
1 ,J 


<WM’  •> 

N 1 ,J 


(29) 


E-  ■ = the  j-th  component  of  E.  in  (24), 

1 l ,tan 
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W2’  *•*  constitute  a set  of  testing  functions  which  are  tangential  to 
S,  and  <a,b>  represents  an  inner  product  defined  as: 

<a,b>  = / ab  ds'  (30) 

The  elements  of  [Zmr)]  in  (28)  may  be  called  generalized  impedances;  those  of 
[Ip]  in  (27)  generalized  currents;  and  those  of  [V  ] in  (29)  generalized 
voltages  representing  the  sources.  When  [Z  ] is  nonsingular  which  is  true  in 
general,  its  inverse  [Zmn]-^  exists.  Then,  we  have 

[IJ  = [Z  ]_1  [V  ] (31) 

nJ  L rnnJ  L mJ  v ‘ 

which  is  the  set  of  coefficients  required  in  (25)  for  obtaining  an  approximate 
component  of  j(r'). 

Specific  numerical  results  of  | J /H^ z | and  |JZ/H^Z|  on  Face  I of  figure  2 
for  (O.U)3  cube  with  the  same  normal  incidence  are  given  in  figures  4 and 
5.  From  figure  4,  it  is  clear  that  I J x / H -j z I 1S  again  the  dominant  response 
for  the  particular  plane-wave  incidence  being  considered.* 

3.  LOSSY  BUILDINGS  WITH  SMALL  APERTURES 

For  electromagnetic  sources  at  VHF  impinging  on  a practical  building 
size,  the  box  represented  in  figure  2 will  be  electrically  large,  and  thus, 
the  methods  outlined  in  Section  2 will  be  very  time  consuming,  if  not 
impossible.  In  this  section,  a different  approach  to  solving  the  problem  of 
large  cavities  is  considered.  In  addition,  apertures  representing  doors  or 
windows,  and  the  case  of  non-perfectly  conducting  building  walls  are  also 
included.  For  the  sake  of  simplicity,  only  the  small  aperture  in  the  form  of 
a horizontal  narrow  slot  in  an  infinite  ground  plane  as  shown  in  figure  6 is 
studied  first. 

The  approach  is  based  on  an  application  of  the  equivalence  principle  L:] 
and  the  use  of  a generalized  network  formulation  [8].  By  this  approach  the 
small  aperture  can  be  replaced  by  an  equivalent  magnetic  current  defined  by: 

M = y x E , (32) 

d 

*Private  communication,  D.K.  Cheng,  Syracuse  University,  Syracuse,  New  York. 
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where  y is  the  unit  vector  in  the  direction  of  propagation  of  the  incident 


The  field  in  region  a of  figure  6 is  that  produced  by  the  incident  wave 
and  that  due  to  M over  the  aperture  with  the  aperture  covered  by  a perfect 
conductor.  The  field  in  region  b is  that  produced  by  -M  over  the  aperture 
with  the  aperture  covered  by  a perfect  conductor.  The  basic  problem  here  is 


component  of  the  magnetic  field  is  continuous  across  the  aperture.  This 
eventually  leads  to  another  integral  equation  involving  the  unknown  M.  The 
technique  similar  to  that  presented  in  Section  2.2  namely, 
approxiinati  ng  M by  with  unknown  coefficients  Vn , may  then  be  utilized. 

For  narrow  slots  of  a length  not  greater  than  three-quarters  of  a 
wavelength,  there  is  evidence  that  M in  the  slot  can  be  closely  approximated 
by  an  equiphasal  half-cosine  wave  [9].  This  suggests  that  we  may  simplify  the 
problem  at  hand  considerably  by  using  a one-term  expansion  for  M and  write 


where  LY  is  the  length  of  the  slot  under  the  condition  Lv  < 3/4  a.  The  matrix 

A X 

equation  is  then  reduced  to  a single  scalar  equation  representing  the 
generalized  Ohm's  law: 


the  determination  of  M,  with  the  boundary  condition  that  the  tangential 


M = V M 


o 


= X V cos  (ttX/L  ) 
X 


(33) 


( ya  + Yb)  v = I 


(34) 


where  V is  a generalized  voltage  to  be  determined,  Yd  and  Yb  are  generalized 
admittances  (in  the  si ot-aperture)  of  the  half-space  and  the  cavity 
respectively,  and  I is  the  generalized  current  source  given  by: 


1 = 2 f M . H . ds , 


J 0 1 
apert. 


(35) 


with  H representing  the  incident  magnetic  field  over  the  si ot-aperture  in 
free  space. 
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In  (34)  the  generalized  half-space  admittance  in  region  a is  [8]. 

Y3  = -2  / Mq  . H (Mq)  ds  (36) 

apert . 

where  H (M  ) is  the  magnetic  field  produced  and  radiated  to  free  space  by 
In  region  b,  the  admittance  may  be  expressed  in  terms  of  the  mode 
characteri sti cs  in  the  cavity  [9], 

Yb  = -j  z V Y.  cot  k.d  (37) 

where  and  k-  denote,  respectively,  the  characteri stic  admittance  and  the 
wavenumber  of  the  ith  mode,  and  the  coefficient  A.j  is  given  by: 

A.  = f M . y x e . ds , (38) 

apert . 

with  e.  being  the  normalized  ith-mode  transverse  electric-field  vector. 

If  there  exists  only  one  transmission  mode  in  the  waveguide  in  region  b; 
that  is,  if  all  modes  except  the  i=l  mode  are  cut-off,  then  all  and  k-  are 


imaginary  except  Y^  and  k-p  We  write  [9a] 

Yi  = j6i 

(i  * 1) 

(39) 

ki  = - j“i 

(i  * 1)  . 

(40) 

For  i = 1 , Y-j  is  the  (real)  characteri  sti  c admittance  of  the  dominant  mode,  and 

kl  = 3i  - Jai  , l4l) 

where  ^ is  the  phase  constant  which  is  related  to  the  guide  wavelength, 
of  the  dominant  mode 


3l  = 27r/Agl 

and  oj  is  the  attenuation  constant  due  to  power  loss  in  cavity  walls. 


(42) 


Based  on  the  above,  we  can  construct  a generalized  equivalent  network  for 
the  cavity-backed  si ot-aperture  shown  in  figure  6.  The  network  consists  of  a 
current  source  I in  parallel  with  Yd,  together  with 

j B = j z Bi  coth  a.d  (43) 

i*l 

and  a terminated  transmission  line  coupled  through  an  ideal  transformer  with  a 
turns-ratio  Ap  This  is  given  in  figure  7. 

In  order  to  study  cavity  fields,  it  is  essential  to  consider  that  the 
equivalent  transmission  line  is  terminated  by  a nonzero  surface  impedance 


zs  - RS  + * d + j)  T 


= (1  + j) 


TT  f y 

a 


(44) 


where  a is  the  conductivity  of  the  cavity  material,  and  for  the  dominant 
TE^q  mode  in  an  a x b rectangular  waveguide  is  given  by 


ai 


nb 


/ 1 - x/2a) 


[1 


(45) 


n being  the  intrinsic  impedance  of  the  medium  inside  the  cavity.  For  air,  n 
is  approximately  377  n. 

For  cavity  walls  made  of  a good  conductor  the  reactive  part  of  Zs,  X$ , 
has  an  insignificant  effect  on  the  resonant  frequency  of  the  cavity  and  may  be 
neglected.  The  real  part,  R$,  however  influences  the  Q and  field  calculations 
in  a major  way;  it  must  be  retained,  as  was  done  in  figure  7. 

The  input  admittance  of  a lossy  transmission  line  of  characteri  Stic 
admittance  Yp  wavenumber  kj  = - jap  length  d,  and  terminated  in  Rs  is 


Y. 

in 


r 1 + tan  ( 3 1d  - j ot x d ) 

\ Y1RS  + j tan  (3^  - ja^)  . 


(46) 
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At  natural  resonance. 


0 id  - n-n,  n = 1 ,2, . . . 


(47) 


Equation  (46)  becomes 


< 1 + Y.R  tanh  aid  \ 

Y = Y s l 

in  1 { Y1ks  + j tan  a]d 


(48) 


The  quantity  ot^d  is  usually  very  small  for  highly  conducting  walls,  and  it  is 
safe  to  approximate  tanh  cqd  by  cqd  and  to  neglect  the  term  Y^RS  tanh  a^d  in 
comparison  with  1.  We  have,  from  (48) 


Y. 

• i n 


(49) 


where 


rd  - ♦ aid 


(50) 


is  a dimensionless  quantity.  Thus  at  natural  resonance,  Y^n  is  purely  real. 

Tne  input  impedance  of  the  low-loss  transmission  line  terminated  in  R$  in 
the  neighborhood  of  natural  resonance  is  (^d  not  exactly  equal  to  pn) 


where 


Z. 


in 


Y^—  = R + jX 
' i n 


(51) 


R = 


rd/Yl  ‘ rdZl 


(52) 


and 


X = Zjtan  B^d. 


(53) 


The  Q of  the  cavity  can  be  obtained  from  the  definition 


dX 

duj 


(54) 
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Using  (52)  and  (53)  for  the  dominant  mode,  we  have 


Q = 7^ — Hr)2  secVd  (55) 

rdAgl  A 1 

where  the  relation 

»gl  * * 1 ' (56> 

has  been  used. 

At  natural  resonance. 


d 


n=l ,2 , 


(57) 


Equation  (55)  becomes 


Q = 


TT  d 

Vgl 


A 

(- 


hTT 

2rd 


<r,2 


(58) 


In  order  to  determine  the  fields  in  the  si ot-aperture  and  inside  the 
cavity,  we  must  find  V from  (34).  In  other  words  we  need  to  find  I,  Ya  and 
Yb.  Assuming  an  incident  plane  wave  with 


H. 

l 


xH  e 
o 


-jky 


-jky 


(59) 


where  k = 2tt/a,  I can  be  evaluated  easily  for  a narrow  Lx  x Ly  slot  by  using 
(33)  and  (59)  in  (35). 


I 


_4 
tt  n 


L L E 
x y o 


(60) 


The  equivalent  half-space  admittance  Ya  can  be  computed  from  (36)  in 
terms  of  an  electric  vector  potential  F and  a magnetic  scalar  potential  <f>m  due 
to  the  magnetic  current  MQ.  That  is, 

Y3  = 2 / M . (jwF  + v$  ) ds  (61) 

o m 

apert . 
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where 


F = — f m exP(-jk^  - r 
4tt  ^ o - - 

apert.  r - r' I 


ds- 


(62) 


and 


J 


4tt  p go 

o apert. 


/ (V-  • M ) 


exp(-jk|r  - r 


ds 


(63) 


These  integrals  have  to  be  evaluated  numerically. 

The  cavity  admittance  is  obtained  by  combining  (37),  (43),  and  (49) 


b 2 

Y = A,  Y.  + jB. 
1 in  J 


(64) 


It  is  also  necessary  to  find  the  turns-ratio  of  the  ideal  transformer  shown 
in  figure  7 for  the  dominant  mode  as  well  as  the  A-  (i  = 2,3...)  for  the 
higher-order  modes  in  (43)  from  (38).  For  an  a x b rectangular  waveguide,  the 
normalized  ith  mode  function  for  TEmn  modes  (m+n/O)  is  [10], 


where 


e . = 

i 


eomeon 


ab 


[(f)  + 02rl/2 


/- nir \ /muxf  . s 

. [x  (^-)  cos  (- ) sin  (— g-) 


~ /mii>  • ,mirX''x  <nTiZ'n 

+ y (— ) sin  (-g— ) cos  (-g— )], 


(65) 


eom 


/ 1,  m = 0 
■2,  in  = 1 ,2, . . . 


(66) 
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The  normalized  ith  mode  function  for  TMnm  modes  (m/0  and  n/0)  is  [10] 


ei 


■J  ab 


t(^)  + 02rl/1 


r /II  ITT  \ / 1 

.[  X (— ) cos  (- 


\ • t iittz  \ 

-)  si n (-—=-) 


(67) 


-z  (t-~)  sin  (— — ) cos  ( — K— )] » 


In  (65)  and  (67),  x'  = x + a/2  and  z'  = z + b/2.  Substitution  of  (33),  (65) 
and  (67)  in  (38)  and  integration  yield  the  values  of  for  all  TE  and  TM 

modes. 

The  wavenumber,  k.j , of  the  ith  mode  can  be  expressed  in  terms  of  its  cut- 
off wavelength  X-j c and  the  correspondi ng  cut-off  wavenumber  k^, 


where 


k/l  - <kAic>2  • * < *ic 


-jk.Jl  - (x.c/x)\  \>xu 


(68) 

(69) 


(70) 


The  characteristic  admittance  for  the 


ith  mode  is 


k.j  /oou 

5 

TE  Modes 

(71) 

we/k . 

9 

1M  Modes. 

(72) 

Having  found  the  values  of  An- , k^  and  for  the  dominant  and  higher-order  TE 
and  TM  modes,  we  are  then  in  a position  to  determine  Yb  from  (64)  in 
conjunction  with  (43)  and  (49).  With  Yd  and  Yb  known,  the  generalized  voltage 
V can  be  found  from  (34)  using  the  value  of  I in  (60). 
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Tne  transverse  electric  field  inside  the  cavity  for  the  dominant  TE1Q 
mode  is 


where 


E = zVA  C cos  (— ) 
It  1 a 


sin  k (d-z) 

ttX  , 1' 


sin  k ^d 


(73) 


C 


(74) 


The  field  distribution  for  the  higher-order  modes  has  the  following 
expressions : 


E. 

it 


E VA.e. 


i ^ 1 


l l 


sin  k . (d-z) 

i 

sin  k d . 


(75) 


These  field  distributions  can  be  computed  for  a rectangular  cavity  of  given 
dimensions,  made  of  a material  with  known  conductivity  and  having  a slot- 
aperture  of  specified  width  and  length. 

We  now  present  some  numerical  results  for  the  following  aspects  of  a 
cavity-backed  si ot-aperture  problem  at  natural  resonance:  cavity  Q,  field 
amplitude  in  the  slot-aperture,  and  field  distribution  inside  the  cavity.  The 
effects  of  the  cavity  material  on  these  results  are  considered.  For 
computational  purposes  we  choose  the  following  parameters: 


Frequency 
Wavel ength : 

Cavity  dimensions: 

Aperture  dimensions: 


f = 100  MHz 
A = 3 m 
a = 0.762a 
b = 0.339a 

d = A g i / 2 = 0.663a 
Lx=  0.25a 
Ly=  0.02a 


17 


The  cavity  quality  factor  for  the  dominant  TE-^q  mode,  Q,  is  calculated  for 
three  different  wall  materials;  namely  copper,  aluminum  and  iron,  using  (58) 
with  n = 1,  which  are  given  in  Table  1. 


Table  1 

Cavity  quality  factor,  Q,  for  dominant  mode  at 
natural  resonance,  (x  = 3 m ) 


Cavity  material 

Copper 
A1 umi num 
I ron 


Conductivity,  g 

5.80x1 07  S/m 
3.54xl07  S/m 
0.95x1 07  S/m 


Quality  factor,  Q 

9. 19xl04 
7 . 1 8x1 04 
3. 72xl04 


The  electric-field  distribution  in  the  si ot-aperture  at  z = 0 is 
approximately  a half-cosine  whose  normalized  amplitudes,  |E  |/E  , are 
calculated  in  accordance  with  (73)  for  the  three  wall  materials  and  tabulated 
in  Table  2. 


Tabl 

Normalized  amplitude  of 
si ot-aperture  (x=y=z=0) 

Cavity  Material 

Copper 
A1 umi num 
Iron 


_2 

electric  field  in  the 
at  natural  resonance. 

Normalized  electric-field 
ampl i tude , | E$ | /EQ 

8.1 6xl0“5 
10.45xl0“5 
20. 17x10“'* 


The  amplitude  of  the  electric  field  in  the  aperture  is  seen  to  be  several 
orders  of  magnitude  smaller  than  that  of  the  incident  wave.  It  is  most 
interesting  to  note  that  the  field  there  is  almost  exactly  inversely 
proportional  to  the  square-root  of  the  conductivity  of  wall  material. 
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The  amplitude  of  the  electric  field  inside  the  cavity  at  y = d/2  has  also 
been  calculated  for  different  wall  materials.  For  the  present  example,  it  is 
several  orders  of  magnitude  larger  than  that  at  the  si ot-aperture ; it  is  more 
than  three  times  that  of  the  incident  wave;  and  it  is  practically  independent 
of  the  conductivity  of  the  wall  material.  Table  3 is  included  to  emphasize 
these  points. 


Table  3 

Normalized  amplitude  of  electric  field  at  the  cavity  center, 
x=z=0,  y=  d/2  at  natural  resonance 

Normalized  electric-field 

Cavity  Material  amplitude,  |E  |/E 

Copper 

Aluminum  3.29 

I ron 

Figure  8 shows  a plot  of  |Elt|/Eg  along  the  center  axis  x = z = 0,  indicating 
that  it  is  at  least  possible  to  produce  a high  field  strength  inside  the 
cavity.  Note,  however,  that  the  specific  numerical  results  given  thus  far  are 
good  only  at  the  first  natural  resonance  when  the  dominant  TE^g  mode  exists 
and  the  slot  is  not  too  long.  For  larger  slots,  more  terms  should  be  added  in 
(33)  and  the  computations  will  be  much  more  complicated. 
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4.  LARGE  BUILDINGS  AND  APERTURES  WITH  LOSSY  GROUNDS— A HYBRID  APPROACH 


The  analysis  given  in  the  preceding  section  is  valid  only  for  small 
apertures  and  for  the  building  (cavity)  in  free  space.  With  buildings  of 
practical  sizes  for  the  purpose  of  storage,  the  buildings  and  apertures  (doors 
or  windows)  are  usually  large  at  VIH.  Furthermore,  the  buildings  and  the 
exterior  transmitters  are  above  an  imperfect  ground  and  the  building  walls  may 
be  made  of  some  reinforced  material  which  can  hardly  be  considered  as 
perfectly  conducting.  In  view  of  the  outlines  presented  in  Section  2,  a 
rigorous  analysis  of  a practical  three-dimensional  model  based  on  the 
Maxwell's  equations  and  appropriate  boundary  conditions  is  certainly  not  a 
simple  task.  For  this  reason,  a hybrid  approach — a combination  of  equations 
and  measured  data — together  with  some  necessary  simplifications  and  practical 
restrictions  is  proposed  in  this  section. 

The  field  distribution  inside  an  empty  building  is  considered  to  be 
produced  by  a direct  component  due  to  the  source  being  operated  outside  the 
building  and  a scattered  component  re-radiated  by  an  equivalent  set  of  surface 
current  densities  induced  on  the  building  walls  (including  the  floor  and 
ceiling)  which  would  exist  if  the  wall  surfaces  were  infinitely  thin  and 
perfectly  conducting. 

The  geometry  and  coordinates  used  here  are  shown  in  figure  9,  where  the 
building  is  hollow,  situated  on  an  imperfectly  conducting  ground,  with  only 
one  door  (aperture)  appearing  on  Plane  m^m^m^m^  (xz-plane)  or  Plane  3.  The 
floor  plane  (m^m^nigm,-)  is  designated  as  Plane  1,  the  front  (m^m^nigm,-)  as  Plane 
2,  the  ceiling  (m^m^m^mg)  as  Plane  4,  the  back  (m^m^m^mg)  as  Plane  5,  and  the 
left  (rTigirigm^mg)  as  Plane  6.  There  is  a vertically  polarized  signal  from  a 
transmitter  located  at  x = Dc,  y = -yt,  and  z = h^ , where  Dc  is  the  door 
center  at  floor  level  as  indicated  in  figure  9. 
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Tne  basic  equation  needed  here  for  calculating  the  scattered  field  is  the 
combination  of  (20),  (21),  and  (22): 


Es(0  = - / J (r1 ) G (r,  r' ) ds ' 

+ jdbl  v / v'  • J (^')  G (r»  r')  ds* 

= . j j (r- ) G (f } f- ) ds' 

- — ^ V / J (r1)  . v‘  G (r,  r')  ds 1 . (76) 


After  approximating  the  double  integral  in  (76)  with  a double  summation  by 
subdividing  each  of  the  six  planes  into  m x n rectangular  patches,  assuming  a 
constant  unknown  value  for  the  surface  current  density  on  each  patch  and 
expressing  the  final  result  into  components,  we  obtain 


j4„E  (x,y,z)  = mill  4mn  Gmnl  AXAy 
m n 


in  n 


(3) 
xmn  ~rnn  3 


+ coy  l l J'  ' G „ AZAX 


- — III 

I 


(1)  1+JkRmnl 


. xmn  n2 
m n ’ Rmnl 


+ r(x-x’)  + (y-y')<f  *da  ,Ig  .&x&y 

LV  nr  xmn  ' ny  ymn  inn  1 j mnl  ^ 


_ J_  y y r(y-y')  j(2)  + (z-z1)  ] a 0 G 0 AyAz 

we  L L uymn  v n;  zmn  J mn2  mn 2 


m n 


(j0£ 


m n 


nr  ymn 
1+jkR 


-Z7ll\  4mn  :^!l  + C(z-^Jz™n+(x-^)Jx™>inn3 


R 


mn  3 


Gmn3  AZAX> 

(77) 
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J 4TTEy  (x  ,y 


and 

j 4tt E z (x, 


Z)  = “P  I I jiy2  Gmnl 
m n J 

* mil  Jy2  Gmn2  AyAZ 
m n 

- — l l \ J(1)  1 + [(x-xO  J(1)  +(y-y^)J(1)]b  Ig  -I AxAy 

we  L L | ymn  v nr  xmn  J u'  ymnJ  mnlj  inn  1 J 

m n K mnl 

- h { J)mn  ^2  mn2  + C(y-ym>Jimn  + < z' zn  > J imn  3bmn  2 \ Gmn2AyAZ 

m n mn2 

- ~ l l T (z_z')  + (X_XJ)  b 7]azax,  (78) 

we  m n m zmn  V n xmn  mn3j  v ' 


'■z>  - mil  42„>  Gmn2  AyAz 


m n 


(3) 

+ mil  Grnn3  AZAX 

m n 


- ~ y y ((x-x'')  ^ + iy-y')  c , G n AXAy 

008  m n ^ xmn  J n ^mn  innl  mn^  ^ 


l + r(y_r)J(2)  + (z_z0J(2) 

we  t £■  1 zmn  ^2  J ymn  v n' 


m n 


(2)  UJkRmn2 


zmn 


zmn 


mn2 


]cmn2j 


Gmn2AyAZ 


— I l jj 

we  m n I 


(3)  1+jkRmn3 
zmn  p 2 

Kmn3 


|(z-z^)J^3)  + (x-x")  J ^ 3 c 01g  qazax, 
O nr  zmn  v n'  xmn/1  mn3j  mn3 


(79) 
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where 


amnl 

= (x  - 

X') 

ITT 

si 

3 

^mnl 

= (y  - 

yn> 

si 

3 

cmnl 

= zSx 

5 

amn2 

= xS^ 

3 

^mn2 

= (y  - 

S2 

3 

Cmn  2 

= (z  - 

ZS> 

S2 

3 

dmn3 

= (X  - 

S3 

3 

^mn3 

■ ys3 

3 

cmn3 

= (Z  - 

zm> 

S3 

3 

2 2 
k R . 
mm 

- 

3 - j 

3kR  . 
mm 

R . 
mm 


(80) 


(81) 


k 


fMHz 


2n  2tt  f 
" X '300  MHz’ 

= operating  frequency  in  MHz 


(82) 
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and 


x'  = (m  - 1)  a x + x 
m v ' o 


x;  = (n  - 1)  A x + xQ 


y‘  = (m  - 1)  a y + y 
J m v ' o 


y n = (n  - 1)  A y + y0 


m,n  = 1,2,3 


(83) 


xQ,  yQ,  and  zq  being  the  possible  distances  from  which  the  patches  are 
counted.  Note  that,  for  simplicity,  we  have  dropped  the  subscript  s in  (77) 
through  (79). 

A detailed  explanation  for  (77)  is  now  in  order.  The  quantity  Ex  on  the 
left-hand  side  represents  the  x-component  of  the  scattered  complex  field 
strength  inside  the  building  at  a specified  location  (x,  y,  z)  contributed  by 
the  induced  surface  currents.  In  reality,  this  is  supposed  to  be  known  as  it 
can  be  obtained  by  measurements.  Let  the  amplitude  and  phase  be  denoted  by 
V , and  ^ for  the  measured  voltage  due  to  the  x-component  of  the  total  field 
at  (x,  y,  z).  The  amplitude  of  the  x-component  of  the  total  field  itself  is 
then  |E.  | = V /i  , where  the  effective  length  of  the  dipole,  % , is 
approximately  equal  to  (0.6-0. 9)  of  the  actual  dipole  length,  depending  on  the 
frequency.  The  phase  information  can  be  secured  by  using  a reference 

antenna  and  a vector  voltmeter.  This  implies  that  the  x-component  of  the 


V 


total 


If  E.  represents  the  incident  field  from  the 


source,  and  E scattered  field  in  (76),  the  total  field  becomes 
s v ' 


E 


t 


E.  + E , or  E = E - E. 
is  s t 1 


(84) 
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The  x-component  of  (84)  is,  therefore. 


E = 
x 


(85) 


where  E^x  has  been  neglected  in  view  of  the  fact  that  the  incident  field  is 
assumed  z-polarized. 

Separating  (85)  into  the  real  and  imaginary  parts,  we  have 

V 

Re  Ex  = cos  *x  » (36) 

X 

and 


V 

Im  E = — si  n . 

X tx 


(87) 


The  first  term  on  the  right-hand  side  in  (77)  is  the  x-coinponent  field 

contributed  by  the  x-component  of  J (r')  on  the  floor  (xy-plane  where  z1  = 0, 

or  Plane  1)  with  representing  the  x-component  of  the  source  concentrated 

t h 

at  the  center  of  the  mnLn  patch  in  Plane  1,  Gmnl  = exp(-jkRmnl)/Rmnl , 

2 2 2 1 /2 

Rmni=  [(*  - x')  + (y  - y„)  + z ] and  (xj,  y'.  0)  representing  the 

location  of  the  mnL  patch.  The  second  term  on  the  right-hand  side  is  the 

m 

contribution  due  to  the  x-component  of  the  surface  current,  J^,  on  Plane  3 

(xz-plane)  where  Gmn3  = exp(-jkRmn3)/Rmn3,  Rfnn3  = [(x  - x^)2  + y2  + 

2 1/2 

(z  - z')  1 7 . The  first  two  terms  on  the  right-hand  side  of  (77)  are, 

' nr 

therefore,  directly  from  the  first  term  of  (76). 

The  remaining  three  terms  in  (77),  derived  from  the  second  term  of  (76) 
after  performing  the  necessary  mathematical  operations  of  divergence  ana 
gradient,  represent  contributions  from  various  components  of  the  surface 
currents  on  Planes  1,  2 (yz-plane),  and  3 (zx-plane).  Similar  explanations 
for  the  various  terms  in  (78)  and  (79)  can  also  be  made. 

In  (78),  we  note,  for  the  left-hand  side,  that 


E 


y 


9 


(88) 


25 


yielding 


V 


Re  E = 


_ y 


— COS  tf, 

*y  y 


and 


Im  E 


— — sin  \ii 
ix, 


(89) 


where  and  ^ are  respectively  the  measured  amplitude  in  volts  and  phase  in 


degrees  due  to  the  y-component  of  the  total  field,  and  % = % . Once  again, 

y x 


E.  has  been  omitted. 
iy 

In  (79),  we  note,  for  the  left-hand  side,  that 


V 


E = E.  - E.  = 
Z tz  1 z 


i,  - E. 

rz  1Z 


(90) 


where  Vz  and  ^ are  the  respective  measured  voltage  amplitude  in  volts  and 


'z  ~x 

the  z-component  of  the  incident  field  above  a lossy  ground  is  approximately 


phase  in  degrees  due  to  the  total  z-component  of  the  total  field,  and 

the  z-comi 
given  by: 

730  PG 


E.  =•!  1 + R exp 
iz  1 v K 


[-jMrrr)]]' 


exp  ( — j k r)  exp  (jkr0) 


(91) 


where 


r1 

h 


input  power  to  the  transmitter,  in  watts, 
transmitting  antenna  gain  (numeric) 

distance  between  the  transmitter  at  (D  , -y^,  hfc)  and  the 

measurement  point  at  (x,  y,  h ) 

2 2 m 2 1/2 

[(x-Dc)  + (y  + yt)  + (hm  - ht)  J 7 , in  meters, 

distance  between  the  transmitter  image  and  the  measurement  point 

2i  1 /2 


(< 


r = 


R,,  = 


x-E>c)  + (y  + yt)  + (hm  + ht)  J 7 , in  meters, 
feed  point  of  the  measuring  antenna  above  ground,  in  meters, 
distance  between  the  transmitter  and  the  reference  antenna,  in  meters, 
Vertically  polarized  ground  reflection  coefficient 
(cose  - n' )/ (cose  + n' ) , (92) 


and 


90°  - tan' 


h 

\ m 


ht)/  /(x  - Dc)2  + (y  + y t ) 2 J . 


(93) 
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In  (92),  n is  a complex  factor  involving  the  ground  constants,  which  may  be 
determined  empirically  by  measuring  the  z-component  field  strength  called  E'  at  a 
convenient  point.  This  convenient  point  may  be  at  the  door  aperture  (D  , 0,  h ) when 
the  door  is  electrically  large.  Thus, 


E'  = | E 1 1 - 

z 1 z 


no 


Az  = E1  + Rvexp(-jkrJ  + jkr')]7^—^  exp(-jkr')exp(jkro),  (94) 


where 


/ 2 2 

'yt  + (hm  - ht)  , 


r'  = Jy 


i Jh  + <hm  f ht>2 


(95) 


From  (94)  an  approximate  value  for  Rv  (here  called  Rj)  at  this  particular  convenient 
point  may  be  determined.  The  angle  e at  this  particular  point  becomes 


S'  = 90»  . tan-1((hm  - ht)/yj 


(96) 


Substituting  (96)  and  R^  obtained  from  (94)  into  (92),  we  can  compute 


n'  = cose'  (1 


R')/(l  + R0  • 

V 7 v V7 


(97) 


With  n1  so  computed,  the  ground  reflection  coefficient  Rv  for  a general  measurement 

location  can  be  determined  from  (92)  and  (93). 

Note  that  in  (77)  to  (79),  the  correspond!- ng  contributions  from  Planes  4, 

5,  and  6 have  been  ignored  in  order  to  reduce  the  number  of  equations 

involved.  This  can  be  justified  as  long  as  the  measurement  data  to  be  used 

for  E , E , and  E are  taken  from  the  region  near  the  origin  m1  (see  figure 
x y z i 

9).  Of  course,  when  the  measurement  data  taken  from  the  region  near  the 

corner  m?  are  to  be  used.  Plane  5 (back)  will  be  more  important  than  Plane  2 
L (2) 

(front).  Under  this  ci rcumstance,  we  replace  in  (77)  to  (79)  by 


,(5)  i(2) 

zmn 

b*  cmn5  = <z 


Jymn’  Jzmn  by  Jzmn’  amn2  by  a™n5  “ W'S5’  bmn2  by  bmn5  ' y.n ^ S5 *c-nn2 


2 2 


zn)S5  5 S2  by  S5  Rmn 5 3 " j 3kRmn5^ /Rmn5 ’ and  Rnn0  by 


R 


mn5 


= ((X  - W)2  + (y  - y')  + (Z  - Z')J 


1/2 


where  W is  the  width  of  the 


building  (see  figure  9). 
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When  the  measurement  data  are  taken  from  spaces  where  contributions  from 

(4  \ 

all  the  planes  are  considered  equally  important  , the  contributions  by  J^'and 

from  Plane  4 (ceiling),  and  from  Plane  5 (back),  and 

ymn  x ymn  zrnn  v ' 

and  from  Plane  6 (left)  can  be  easily  added  to  (77)  through  (79), 

with  similar  definitions  for  amni , bmni  , cmni , Rmni , and  Smni , i = 4,  5,  and 

6.  Of  course,  the  number  of  unknowns  will  then  be  doubled. 

When  the  assumption  made  in  (77)  through  (79)  is  justified,  there  will  be 
six  different  complex  surface  currents  with  two  from  each  plane,  yielding  1 2mn 
real  unknowns.  The  left-hand  side  of  (77)  through  (79)  gives  six  data  values 
(ReE.j  and  ImE^,  i = x,  y,  and  z)  at  each  measurement  point.  Thus,  if  p is 
designated  as  the  number  of  locations  at  which  measurements  of  three 
orthogonal  components  (both  amplitude  and  phase)  are  made,  there  will  be  6p 
equations.  Then,  we  require 

6p  > 12  inn  or  p > 2 mn  , (98) 

in  order  to  be  able  to  solve  for  the  system  of  equations  by  matrix  inversion. 
In  (98),  we  have  assumed  that  the  same  number  of  patches  (mn)  are  considered 
for  each  of  the  three  planes.  On  the  other  hand,  should  a different  number  of 
patches  be  taken  into  account  for  different  planes,  then  we  have, 

6p  > 4 jm  (1)n(1)  + m(2)n(2)  + m(3)n(3)]  , 
or 

p >|  +m<2>n(2>  +m(3>n<3>]  , (99) 

where  m^and  n^)  are  the  number  of  patches  considered  from  Plane  i, 
i=  1,  2,  and  3. 
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It  should  be  emphasized  that  although  (77)  through  (79)  are  important  for 
computing  piecewisely  the  unknown  surface  currents  with  a set  of  measurement 
data  taken  at  some  points  near  the  planes  of  interest,  the  same  equations  are 
equally  useful  for  computing  the  field  strengths  at  points  other  than  those 
data  points  as  soon  as  the  unknown  surface  currents  are  determined.  In  this 
latter  process,  the  computation  is  strai ghtforward  and  does  not  involve  matrix 
inversions.  For  this  reason,  this  latter  computation  is  designated  as  the 
forward  computation. 

Naturally,  the  larger  the  number  of  patches  to  be  considered  at  a time, 
the  more  accurate  are  the  numerical  results.  The  actual  number  of  patches 
used  is  ultimately  determined  by  compromising  the  required  accuracy  and  the 
computing  capacity  and  cost.  Since  there  is  no  a priori  knowledge  about  the 
accuracy  of  the  results,  we  recommend  to  begin  with  a small  number  of  in  and  n, 
and  then  proceed  to  obtain  a solution  for  the  surface  currents.  These 
currents  can,  in  turn,  be  used  to  perform  a forward  computation,  also  using 
(77)  through  (79),  for  the  field  at  different  locations  where  the  measured 
data  are  also  available.  A comparison  of  the  computed  and  measured  results  on 
fields  should  reveal  the  degree  of  accuracy  by  the  procedure  so  outlined.  If 
the  accuracy  is  at  an  acceptable  level,  the  originally  assumed  small  number 
for  m and  n is  considered  satisfactory.  Otherwise,  a large  number  for  m and  n 
has  to  be  used  to  improve  accuracy. 

This  process  may  be  verified  by  on-site  measurements  if  a small  computer 
with  reasonable  capacity  is  available. 

In  principle,  whenever  the  equal  sign  holds  in  (98)  or  (99),  a unique 
solution  is  always  obtainable.  In  practice,  however,  we  do  prefer  to  have  p 
larger  than  the  minimum  required  such  that  the  number  of  equations  is  larger 
than  the  number  of  unknowns.  Under  this  condition,  the  solution  obtained  in 
the  sense  of  least  square  fit  will  minimize  large  deviations  which  may  result 
because  of  the  potential  errors  contained  in  the  measured  data. 
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With  an  accurate  and  approximate  set  of  piecewise  solutions  for  the 
surface  currents  obtained  from  (77)  through  (79),  we  then  use  them  to  compute 
the  field  strength  everywhere  else  inside  the  building  so  that  the 
el ectromagnetic  environment  of  the  entire  building  may  be  assessed.  In 
general,  the  space  near  the  door  (when  it  is  open)  for  a large  building  and 
door  should  experience  relatively  high  field  strengths.  Therefore,  careful 
and  regular  measurements  should  be  made  for  the  space  close  to  Planes  2,  3, 
and  5 in  order  to  obtain  more  accurate  results  for  the  surface  currents  on 
Planes  1,  2,  3,  and  5 near  the  door,  which  are  important  to  compute  the  field 
strengths  for  the  front  space  near  the  door.  With  such  a hybrid  analytical 
procedure  one  can  avoid  a large  scale  of  tedious  and  systematic  measurements, 
which  would  normally  be  required.  Sample  numerical  results  using  this  hybrid 
procedure  are  given  in  Section  7. 
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5.  MEASUREMENT  REQUIREMENTS  AND  PROCEDURES 


The  equipment  required  for  measuring  the  field  strengths  consists  of 
three  orthogonal  dipoles,  a pole  for  supporting  the  antennas,  an  antenna 
selector  switch  box,  a vector  voltmeter,  a three-channel  digital 
voltmeter/printer,  and  a reference  antenna  (monopole),  connected  as  given  in 
figure  10.  A movable  platform  on  which  the  antennas  and  the  selector  switch 
box  are  mounted  is  also  required.  The  digital  vol tmeter/pri nter  prints  out 
three  amplitudes  in  volts  and  three  phases  in  degrees. 

The  measurement  procedures  may  be  outlined  as  follows: 

(1)  Select  an  origin  for  the  building  as  shown  in  figure  9.  Identify 
x,  y,  and  z directions. 

(2)  Measure  the  building  width  W,  depth  D,  and  height  H}all  in 
meters.  These  may  be  obtained  from  the  base  where  the  building  is 
being  investigated.  When  H is  larger  than  5 meters,  it  may  be 
ignored  since  the  contribution  from  Plane  4 (ceiling)  may  be 

del eted. 

(3)  Measure  the  width  and  height  in  meters  of  the  door.  Determine  Dc 
(see  figure  9)  in  meters. 

(4)  Divide  the  floor  (Plane  1)  into  regular  grids  by  drawing  straight 
lines  parallel  to  the  x-  and  y-  axes  such  that  ax  and  Ay  are  x/4  or 
less,  where  x is  the  operating  wavelength.  Let  ax  = Ay  although 
they  are  not  required  to  be  equal  in  principle.  Record  ax  and  Ay 
in  meters.  Observe  the  intersecting  points.  Also  choose  az. 

(5)  Connect  the  measurement  equipment  as  shown  in  figure  10,  with  one 
antenna  in  x-direction,  another  in  y-direction,  and  a third  in  z- 
di recti  on  (vertical). 

(6)  Record  the  feed  height  h of  the  measurement  antennas  above  the 
floor,  in  meters. 

(7)  Set  the  transmitter  outside  the  building  at  x = Dc,  y = -yt,  and  z 
= ht.  Record  yt  and  ht  in  meters. 
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(8)  Set  the  reference  antenna  at  a convenient  location  outside  the 
building.  Record  the  distance  rQ  in  meters  between  the  transmitter 
and  the  reference  antenna. 

(9)  Record  the  frequency  f in  MHz  of  the  transmitting  antenna  which  is 
assumed  to  be  vertically  polarized  (E^z  only),  the  input  power  P in 
watts,  and  the  antenna  gain  G (numeric). 

(10)  Determine  xQ,  yQ,  and  zQ  in  meters  required  in  (83).  This  is 
necessary  because  of  the  size  of  the  moving  platform  on  which  the 
measurement  antenna  is  mounted. 

(11)  Mark  the  first  intersecting  point  [obtained  in  (4)]  near  the  origin 
where  measurements  are  possible  as  position  (1,1).  Number  the  next 
intersecting  points  toward  the  +x-direction  as  (2,1,),  (3,1)  etc. 
to  form  the  first  row.  Similarly,  number  the  next  intersecting 
points  toward  the  +y-direction  as  (1,2),  (1,3),  etc.  to  form  the 
first  column. 

(12)  Make  as  many  columns  as  possible  between  Planes  2 and  5.  Make  as 
many  rows  as  possible  up  to  y = D/2. 

(13)  One  person  is  assigned  to  the  transmitter  to  turn  on  and  off  of  the 
transmitting  signal  as  instructed. 

(14)  Another  person  is  assigned  to  record  the  signal  from  the  printer. 

(15)  A third  person  is  assigned  to  align  the  feed  point  of  the 

orthogonal  dipoles  to  position  (1,1).  Once  done,  he  asks  the 

transmitter  operator  outside  the  building  to  key  on  the  signal, 

walks  away  from  the  measuring  antennas,  and  asks  the  printer 

operator  to  record  the  amplitudes  and  phases  of  the  three  component 

voltages  (V  , V , V , ^ , 4/,  41  ) in  volts  and  degrees 
a y z x y z 

respectively,  together  with  the  measurement  position. 
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(16)  Once  the  signals  are  recorded,  the  person  assigned  to  move  the 
measuring  antennas  asks  the  transmitter  operator  to  turn  off  the 
signal.  He  then  moves  the  measuring  antenna  to  position  (2,1), 
signals  the  transmitter  to  be  turned  on,  and  repeats  the  process. 

(17)  Take  and  record  at  least  10  rows  and  5 columns  each  along  Planes  2 
and  5. 

(18)  Take  and  record  a few  points  around  the  center  of  the  building  (x 
W/2,  y = 0/2). 

(19)  Take  and  record  a few  points  at  the  door  aperture  (y  = 0). 

(20)  Raise  the  antenna  feed  point  by  x/4,  and  repeat  the  same  process 
starting  position  (1,1). 
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6.  COMPUTING  REQUIREMENTS  AND  PROCEDURES 


This  section  contains  a discussion  of  four  computer  programs  used  to 
solve  (77),  (78),  and  (79)  given  in  section  4.  The  main  part  of  the 
discussion  deals  with  two  FORTRAN  IV  programs,  one  that  computes  the  currents, 
J , from  the  known  measured  field  strengths,  and  another  that  computes  field 
strengths,  E,  from  given  currents.  A second  pair  of  programs,  using  HPL 
language,  is  based  on  the  FORTRAN  programs  but  is  designed  for  a smaller 
desktop  computer. 

When  computing  currents,  a 6p-by-12mn  matrix  is  formed,  which  can  become 
large  if  p,  m,  and  n are  large  (p,  m,  and  n are  defined  in  section  4).  Thus 
the  computer's  central  memory  storage  capacity  requirement  will  be  large.  For 
example,  let  m and  n be  equal  to  2.  This  defines  a 2 x 2 x 2 -point 
measurement  cell  which  we  consider  to  be  the  smallest  data  set  that  will 
produce  meaningful  results.  The  matrix  size  is  48  x 48  elements.  To  compute 
the  currents  in  Planes  1,  2,  and  3 requires  a central  memory  containing 
approximately  25,000  storage  locations.  If  m and  n are  equal  to  6 (i.e., 
a 6 x 6 x 2 -point  data  set),  the  matrix  size  is  216  x 120  elements  and 
approximately  250,000  storage  locations  are  required.  For  a computer  that  is 
too  small  to  execute  this  program  even  for  the  smallest  (48  x 48)  matrix,  the 
program  is  divided  into  two  parts  as  described  in  section  6.3  below. 

When  computing  field  strength,  no  matrix  operations  are  required,  and  the 
required  amount  of  storage  does  not  strongly  depend  upon  the  data  set  size. 
Approximately  20,000  storage  locations  are  adequate  for  a 6 x 6 x 2 -point 
data  set. 

The  actual  computer  size  required  will  depend  upon  the  individual 
computer  system's  archi tectural  details  such  as  compiler  design,  word  length 
capacity,  exponent  range,  on-line  storage  capability,  etc.  The  large  computer 
used  to  execute  these  programs  at  NBS  has  the  following  characteri sti cs : 


e. 


a. 


d. 


c . 


b. 


NOS  operating  system 

400K  words  (octal)  central  memory  capacity 

60-bit  word  size 

±322  exponent  range 

~107  instructions  per  second 
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6.1.  CURRENT  Program  (FORTRAN) 

The  FORTRAN  program  used  to  solve  (77),  (78),  and  (79)  for  the  induced 
currents  in  Planes  1,  2,  and  3 is  given  in  Appendix  A,  section  A.l.  In  its 
present  form,  this  program  contains  a large  amount  of  redundancy  as  a result 
of  the  decision  to  mimic  the  original  algebraic  equations  of  the  problem. 

This  allowed  the  program  to  be  proof-read  for  accuracy  by  project  personnel 
who  have  only  a casual  acquaintance  with  FORTRAN  language.  For  a large 
program  such  as  this  one,  much  debugging  cost  can  be  saved  by  this  approach. 

A more  efficient  program  could  be  written,  but  little  overall  machine  savings 
would  result  because  the  vast  majority  of  computer  memory  is  consumed  by  the 
6p-by-12mn  matrix.  The  program  as  shown  is  considered  adequate  for  the 
present  purposes. 

The  following  background  discussion  will  be  helpful  for  understanding 
this  program. 

Equations  (77),  (78),  and  (79)  can  be  written  in  matrix  notation  as 
fol 1 ows : 


[E]  = [A]  [J]  , 


(100) 


where  [E]  is  a u-vector,  [3]  is  a v-vector,  and  [A]  is  a u-by-v  matrix.  In 
our  case,  [E]  is  a 6p  - element  column  matrix  (u  = 6p), 


[E] 


P 


E (x,  y,  z) 

X 

E (x,  y,  z) 

y 

Ez  (x,  y,  z) 


[101) 


and  [3]  is  a 12mn  - element  column  matrix  (v  = 12mn),  The  elements  in  matrix 
[A]  are  the  coefficients  of  the  currents  J as  found  in  (77),  (78)  and  ('9'. 

The  solution  of  (100)  for  the  currents,  [3],  makes  use  of  the  "OR 
factorization"  [11]  of  matrix  [A],  thus  converting  (100)  into  a for  1 that  c a 
be  solved  directly  by  back-substitution  [ 1 2 j . In  the  QK  f ac  t.  or  i z a t.  i o e t a d , 
an  orthogonal  matrix,  [Q],  and  an  upper  triangular  matrix,  [R].  are  e 
such  that 

[Q]  [R]  = [A].  (102) 
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An  orthogonal  matrix  has  the  property. 


[Q]t  [Q]  = [I], 


(103) 


where  [Q]T  is  the  transpose  of  [Q],  and  [I]  is  the  unit 


(or  identity)  matrix. 


[I] 


10  0..  . 
0 10... 

0 0.... 

....00 
...010 
...001 


(104) 


An  upper  triangular  matrix  has  the  form. 


[R] 


r r r r r . 

o r r r r . 

o o r r r . 

ooo  r r . 

etc. 


J 


(105) 


in  which  elements  r^j  are  equal  to  zero  when  i > j. 

From  (102)  and  (103) 

[R]  = [Q]T  [A],  (106) 

Multiplying  (100)  through  by  [Q]f  and  applying  (106)  gives 


[Q]T  [E]  = [Q]T  [A]  [J] 

= [R]  [J]. 


[107) 


Because  [R]  is  upper  triangular,  (107)  is  easily  solved  for  the  elements 
of  [J]  by  back-substitution,  starting  with  the  last  (12nm-th)  element.  Thus, 


J12mn 


(q  E) 

r r 
6p, 


6p 

12inn 


(108) 
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is  the  last  element  in 


where  J12mn  is  the  last  element  in  [3],  (q  E)6p 
[Q]  [E],  and  ^2mn  is  the  correspondi nq  element  in  [R].  Having  found 

we  next  eiemeirt  ^I2mn-1’  ^12mn-2  anci  so  on  un1:ii  ai  1 elements 

in  [3]  have  been  found. 

The  solution  thus  obtained  is  a least  square  solution  [13]  of  the  system 
of  equations. 

The  program  as  shown  in  section  A.l  is  based  upon  the  following  assumptions 

a.  The  field  strength  measurement  data  are  stored  in  a computer  file  in 
the  specific  ordered  sequence  as  shown  in  Tables  1 and  2,  section  7.  Each 
data  value  occupies  a specific  position  in  the  data  grid  that  corresponds, 
poi nt-for-poi nt , with  the  measurement  grids  discussed  in  section  5. 

b.  The  intervals  ax  and  Ay  , between  the  data  points,  are  constant  in 
both  the  X-  and  Y-di recti ons . They  are  designated  DX  and  DY  in  the  program, 
and  may  be  either  equal  or  not  equal  to  each  other. 

c.  Measurement  data  are  obtained  at  two  levels  above  Plane  1 (floor). 

The  lower  level  height  above  Plane  1 is  identified  as  ZSTART,  and  the  upper 
level  is  OZ  above  ZSTART. 

d.  Out  of  the  complete  set  of  measurement  data,  the  number  of  data  point 
locations  in  the  X-direction  at  the  lower  (first)  level  that  are  used  in  these 
computations  is  2.  This  is  identified  as  NXTOTH(l)  = 2,  which  freely 
translated,  reads,  "the  TOTal  J^umber  of  _X-di  recti  on  data  locations  in  the 
first  level  Held  for  computation  use  is  equal  to  two."  Similarly, 

NXT0TH(2)  = 2 identifies  that  there  are  two  X-direction  data  locations  at  the 
upper  (second)  level  used  in  the  computations.  NYTOTH(l)  and  NiT0TH(2^  rsfer 
to  the  number  of  Y-direction  data  locations  in  the  first  and  second  levels, 
respectively,  used  in  the  computation.  Bear  in  mind  that  there  are  six  data 
numbers  at  each  location;  viz.,  two  each  (the  real  and  imaginary  partsl  t e 
three  (X,  Y and  Z)  orthogonal  directions. 

e.  Only  the  currents  in  Planes  1,  2 and  3 are  calculated.  his  is 
keeping  with  the  discussion  in  section  4. 
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Assumption  a allows  the  data  at  any  desired  grid  location  to  be  used  in 
the  computation.  The  computer  selects  the  desired  data  according  to  program 
statements  beginning  with  the  letters  IXH  and  IYH.  This  is  discussed  further 
below. 

The  consequence  of  assumption  b is  that  the  coordinates  of  the  data  used 
in  the  computations  are  calculated  by  the  computer  for  use  in  the  program 
rather  than  being  provided  to  the  computer  as  input  data.  The  formulae  for 
coordinate  computation  has  the  form, 

x = (m  - 1)ax  + x(start)  , (109) 

where  m is  an  integer-index  (m  > o)  that  specifies  the  desired  data  point. 

Assumptions  c and  d are  arbitrary  and  can  easily  be  changed  if  desired. 

We  have  executed  this  program  for  data  sets  as  large  as  six  (6)  in  each  of  the 
X-  and  Y-directions  (i.e.,  6 data  point  locations  in  the  X-direction,  6 in  the 
Y-direction,  and  2 in  the  Z-di recti  on ) . If  measurement  data  had  been  taken  at 
more  than  two  levels,  a simple  modification  of  the  program  would  allow  it  to 
be  used.  However,  the  maximum  size  of  the  data  set  is  limited  by  the 
computer's  central  memory  storage  capacity  as  stated  above. 

The  consequence  of  assumption  e is  that  the  program  simply  omits 
statements  that  would  be  used  to  calculate  currents  in  Planes  4,  5,  and  6. 
These  statements  can  be  added  as  desired,  using  the  same  basic  code  as  used 
for  Planes  1,  2,  and  3.  We  have,  for  example,  a second  program,  like  the  one 
shown  in  section  A.l,  that  calculates  currents  in  Planes  1,  3,  and  5. 

For  discussion,  the  program  may  be  divided  into  five  parts  as  follows: 

1.  Lines  1 through  56  contain  basic  information  such  as  variable 
definitions,  dimension  statements,  values  of  constants,  and  selectable 
parameter  values. 

2.  Lines  57  through  237  accomplish  the  computation  of  the  elements  of 
coefficients  matrix  [A].  These  computations,  for  the  most  part, 
involve  the  geometry  of  the  problem,  i.e.,  the  locations  of  the 
currents  relative  to  the  locations  of  the  points  at  which  field 
strength  data  were  measured. 

3.  Lines  238  through  244  deal  with  subroutine  CQRDCT  by  which  the  QR 
factori zation  of  matrix  [A]  is  accomplished. 


38 


4.  Lines  245  through  309  construct  matrix  [E]  by  reading  in  the 
measurement  data  from  the  data  file  and  modifying  it  according  to  (90), 
et  seq.,  of  section  4. 

5.  Lines  310  through  359  (the  end  of  the  program)  accomplish  the  solution 
of  (100)  for  the  currents  by  use  of  subroutine  CQRSLT,  which  is  the 
back-substitution  procedure  for  solving  (107). 

4 finely  detailed  discussion  of  this  program  will  not  be  given  because 
familiarity  with  FORTRAN  programming  is  assumed  on  the  part  of  the  reader. 
However,  certain  features  should  be  briefly  described.  These  are  grouped 
according  to  the  five  parts  of  the  program  listed  above. 

Part  1 (lines  1 through  56).  The  statements  that  begin  with  the  letters 
IXH  and  IYH  determine  which  data  are  used  in  the  computation  of  currents. 

They  do  this  by  defining  an  index  that  is  used  at  various  places  throughout 
the  program  to  select  the  desired  data  values  (from  the  data  file)  and  to 
compute  their  coordinate  locations. 

These  statements  have  the  form, 

IXH(a,b)  = c,  (110) 

where  a,  b and  c are  integers.  The  letter  I means  that  this  is  an  index.  X 
means  it  pertains  to  the  X-direction.  H means  that  the  index  is  held  in 
memory  for  use  as  needed.  Integer  a is  an  ordinal  number  signifying  the 
sequential  order  of  this  index.  For  example,  if  a = 1,  then  this  is  the  " i rst 
index  of  its  particular  kind  that  will  be  used  in  later  do-loops.  Integer  b 
refers  to  measurement  grid  level.  If  b = 1,  the  index  refers  to  data  in  the 
first  (lower)  measurement  grid;  b = 2 refers  to  the  second  (upper)  grid. 
Integer  c is  an  ordinal  number  that  refers  to  the  row  or  column  in  the  chosen 
measurement  grid  where  the  desired  data  point  is  located.  For  example,  if  c = 
1 in  eq.  (110),  the  index  refers  to  the  first  data  point  in  the  X-direction. 
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Here  are  two  examples.  In  words. 


IXH  (1,2)  = 3 

reads,  "the  first  data  point  used  by  the  do-loop  from  the  second  (upper) 
measurement  qrid  is  the  third  point  in  the  X-direction  in  that  measurement 
grid."  Similarly, 

IYH  (2,1)  = 9 

reads,  "the  second  data  point  used  by  the  do-loop  from  the  first  (lower) 
measurement  grid  is  the  ninth  point  in  the  Y-direction  in  that  measurement 
grid." 

The  user  may  change  the  integer  c to  select  any  particular  sets  of  data 
that  he  desires  to  be  used  in  the  computation  of  currents.  The  computer  will 
calculate  the  coordinates  of  those  data  and  use  them  whenever  they  are 
required  in  the  program. 

As  mentioned  above,  the  program  printed  in  section  A.l  is  written  to  call 
two  data  points  from  each  level  in  each  of  the  X-  and  Y-di  recti ons . If  more 
data  points  are  desired,  suitable  IXH  or  IYH  statements  are  added.  For 
example,  if  a third  data  point  is  to  be  added  from  the  first  level  measurement 
grid  in  the  X-direction,  this  statement  would  be  added: 

IXH( 3, 1)  = c,  (111) 

where  c is  defined  above.  In  this  case  the  NXTOTH  (1)  statement  must  be 
changed  to  read 

NXTOTH ( 1 ) = 3, 

since  now  there  are  a total  of  three  data  points  in  the  X-direction  that  will 
be  incorporated  by  the  program.  The  number  to  the  right  of  the  equals  sign  in 
the  NXTOTH  and  NYTOTH  statements  must  be  the  same  as  integer  a in  the 
correspondi ng  IXH  and  IYH  statements. 
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— rt  2 (lines  57  through  237).  Coordinates  of  field  locations  within  the 
building  and  of  current  locations  in  Planes  1,  2,  and  3 are  calculated  in  this 
part  of  the  program.  The  former  are  designated  by  the  letters  X,  Y,  and  Z, 
and  the  latter  by  XP,  YP,  and  ZP. 

A;->o  calculated  here  are  the  constants  a • , b • c • G • R anH 

mm  5 mm  5 Lmm  » hnni’  Kmni 5 dna 

Si#  which  are  defined  in  section  4. 

Using  this  information,  the  computer  calculates  current  coefficients  CJX, 
,.JY,  and  CJZ.  These  coefficients  are  structured  by  the  computer  as  the 
elements  of  matrix  [A].  In  general , the  coefficients  are  complex  and  are 
comprised  of  several  terms.  For  example,  the  coefficient  of  the 
current  is,  from  (77), 


CJX1 


mn 


4tt 


4tTU)£ 


1+jkR 


mn  1 


R 


1,11  1 1 , t „ \ 

o + (X-X"  )a  , 

2 v m ' mnl 


mnl 


Gmnl  ^ 


(112) 


However,  each  coefficient  reduces  to  a single  complex  number  pair. 

For  NXTOTH(l)  = 2,  NXT0TH(2)  = 2,  NYTOTH(l)  = 2,  and  NYT0TH( 2)  = 2, 


and 


p - 2 + 2+  2+  2=8, 

m = 2, 

n = 2, 

6p  = 6 x 8 = 48, 

12mn  = 12  x 2 x 2 = 48. 


Thus  [A]  is  a 48-by-48  element  matrix. 

Although  this  part  of  the  program  appears  to  be  lengthy,  it  actually 
contains  a basic  code  that  is  repeated  three  times,  once  each  for  Planes  1,  2, 
and  3. 

Part  3 (lines  238  through  244).  Subroutine  CQRDCT  is  shown  in  section 
A. 2.  This  is  a compact  form  of  a public-domain  program  called  CQRDC,  which 
stands  for  Complex  QR  Decomposition  [14].  CQRDCT  was  written  especially  for- 
this  project  by  P.  N.  Swarztrauber  [15].  Consult  reference  [14]  for  a 
description  of  this  subroutine. 

Notice  that  the  symbol  A is  used  for  two  different  matrices,  neitner  of 
which  is  the  matrix  [A]  discussed  above.  First,  it  defines  a doubly 
subscripted  complex  array  that  contains  matrix  [A]  of  (100).  Secono  , : 
defines  the  upper  triangular  matrix  [R]  of  (102)  used  in  the  solution  o' 

(107). 
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Part  4 (lines  245  through  309).  The  measurement  data  is  read  from  the  data 
file  (line  288),  converted  to  field  strength  by  dividing  the  voltage  readings 
by  £ , and  corrected  for  the  incident  field  per  (90).  The  three  components  of  the 
scattered  field  are  structured  as  the  elements  of  matrix  [E]. 

Part  5 (lines  310  through  359).  Subroutine  CQRSLT  is  shown  in  section  A. 3. 
This  is  a compact  form  of  a public-domain  program  called  CQRSL,  which  stands  for 
Complex  QR  Solver  [14].  CQRSLT  was  written  especially  for  this  project  by  P.  N. 
Swarztrauber  [15].  Consult  reference  [14]  for  a description  of  this  subroutine. 

The  solution  to  (107),  and  thus  the  solution  to  (100),  is  written  out  as  the 
12mn  elements  of  matrix  [3].  These  are  designated  by  the  symbol  XS. 

6.2  FIELD  Program  (FORTRAN) 

The  FORTRAN  program  used  to  solve  (77),  (78),  and  (79)  for  the  field 
strengths  E,  given  the  induced  currents,  3,  is  given  in  Appendix  A,  section 
A. 2.  This  program  employs  much  that  is  used  in  the  program  listed  in  section  A.l, 
and  like  that  program,  it  was  written  for  greatest  proof-readi ng  ease  rather  than 
for  greatest  efficiency. 

The  form  of  the  program  is  based  upon  the  same  five  assumptions  listed  in 
Section  6.1  with  the  following  exceptions: 

a.  Input  data  is  the  current  matrix  [3],  which  is  stored  in  a computer  file 
in  a specific  ordered  sequence  that  corresponds,  poi nt-for-poi nt,  wi th  the 
projections  of  the  measurement  data  locations. 

b.  The  locations  of  the  computed  field  strengths  can  be  at  any  arbitrary 
points  chosen  by  the  user. 

The  discussions  of  this  program  is  divided  into  the  following  four  parts: 

1.  Lines  1 through  54  contain  the  same  basic  information  discussed  in  section 

6.1. 

2.  Lines  55  through  127  read  the  input  data  from  the  computer  file  that  contain 
matrix  [3],  and  prints  out  the  matrix  elements  for  verification. 

3.  Lines  128  through  289  compute  the  coefficients  matrix  [A],  and  then  solve 
(100)  for  the  elements  of  field  strength  matrix  [E], 

4.  Lines  290  through  349  (the  end  of  the  program)  add  the  incident  field  per 
(90),  and  print  out  the  computed  fields  strength  at  the  selected  field 
point(s). 
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A brief  discussion  of  this  program  follows: 


Part  1 (lines  1 through  54).  This  part  of  the  code  is  identical  to  its 
counterpart  in  the  preceding  section,  except  that  the  coordinates  of  the  point 
at  which  the  user  desires  to  compute  the  field  strength  are  defined  as  XFP, 
YFP,  AND  ZFP.  They  may  be  changed  at  will. 

Part  2 (lines  55  through  127).  The  input  data  are  those  values  of  3 
that  are  selected  according  to  their  coordinate  locations.  These  in  turn, 
are  determined  by  the  indices  IXH  and  IXY  as  discussed  in  section  6.1.  The 
elements  of  the  current  matrix  [3]  that  is  formed  from  the  input  data  are 
defined  by  the  symbol  XS. 

The  coordinates  of  the  current  locations  in  Planes  1,  2,  and  3 are 
calculated  by  the  same  code  discussesd  in  section  6.1.  They  are  designated  by 
the  letters  XP,  YP,  and  ZP. 

Part  3 (lines  128  through  289).  Constants  amni , bmni  , cmni , Gmni , Rmni , 
and  S.j , and  the  current  coefficients  CJX,  CJY,  and  CJZ  are  calculated  by  the 
same  code  discussed  in  section  6.1.  A major  difference  in  the  two  programs 
occurs  at  this  point;  namely,  this  program  computes  the  elements  of  matrix  [A] 
and  solves  for  the  elements  of  matrix  [E]  point-by-point  without  the  need  to 
construct  matrix  [A]  in  toto.  This  greatly  reduces  the  central  memory  size 
requirements  as  compared  with  the  program  of  section  6.1.  At  this  point,  the 
incident  field  contribution  has  not  yet  been  added  in. 

Part  4 (lines  290  through  349).  The  incident  field  strengths  are 
calculated  by  the  same  code  discussed  in  section  6.1,  and  used  to  correct  the 
computed  fields  per  (90). 


6.3  CURRENT  Program  (HPL) 

The  HPL  program  used  to  solve  (77),  (78),  and  (79)  for  the  induced 
currents  in  Planes  1,  2,  and  3 is  given  in  Appendix  B,  section  B.l.  It  is 
divided  into  two  parts  so  that  a small  desktop  computer  can  he  used.  Part  1 
computes  the  elements  of  the  coefficients  matrix  [A],  and  performs  the  OR 
factorization  given  by  (106)  and  (107).  Part  2 computes  the  elements  of  the 
current  matrix  [3]. 

This  program  is  based  on  the  FORTRAN  program  discussed  in  section  6.1. 

It  uses  the  same  computational  approach,  and  the  same  assump,  ions  non  > ov-" 
for  assumption  a (pertaining  to  the  input  data).  In  the  HPL  P'oora  , 
coordinates  of  the  data  points  that  are  to  be  used  in  the  co  'nut a*,  o a e 
first  manually  entered  in  accordance  witli  prompting  instructions  o 1 o \ i *.o  > . o \ 
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the  computer.  Later,  the  field  data  values  are  manually  entered,  again 
following  the  computer's  prompting.  Having  obtained  the  data  in  this  way,  the 
computer  calculates  the  solutions  to  (77),  (78),  and  (79)  by  the  same  methods 
discussed  in  section  6.1.  These  current  values  and  their  coordinates  are 
stored  on  the  HPL  program  tape  for  subsequent  use  by  the  FIELD  program 
described  in  section  6.4,  below. 

A detailed  discussion  of  the  design  of  this  program  will  not  be  given 
because  (a)  familiarity  with  HPL  programming  is  assumed  on  the  part  of  the 
reader,  and  (b)  the  general  design  of  the  program  is  closely  similar  to  the 
corresponding  FORTRAN  program.  Refer  to  the  discussion  given  in  section  6.1 
for  information  on  these  programs. 

Variable  assignments  are  given  in  section  B.2.  Instructions  for 
executing  the  HPL  CURRENT  program  are  given  in  section  B.3. 

6.4  FIELD  Program  (HPL) 

The  HPL  program  used  to  solve  (77),  (78),  and  (79)  for  the  field  strength 
at  any  arbitrary  location  within  the  building,  after  the  currents  are 
obtained,  is  also  given  in  Appendix  B,  section  B.l.  It  is  based  upon  the 
FORTRAN  program  discussed  in  section  6.2.  It  uses  the  same  computational 
approach,  and  the  same  assumptions  apply.  In  the  HPL  program,  the  coordinates 
of  the  field  points  are  manually  entered  in  accordance  with  prompting 
instructions  provided  by  the  computer.  The  current  values  and  their 
coordinates  are  obtained  from  the  CURRENT  program  (see  section  6.3,  above). 

The  general  design  of  this  program  is  closely  similar  to  the 
correspond! ng  FORTRAN  program.  Refer  to  the  discussion  given  in  section  6.2 
for  information  on  these  programs. 

Variable  assignments  are  given  in  section  B.2  (same  as  for  HPL  CURRENT 
program).  Instructions  for  executing  the  HPL  FIELD  program  are  given  in 
section  B.3. 
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6.5  HPL  File  Structure  on  Magnetic  Tape 
Both  the  HPL  CURRENT  and  the  HPL  FIELD  programs  are  on  a single  magnetic 
tape.  They  are  structured  such  that  (a)  each  program  may  be  run  independently 
of  the  other,  or  (b)  the  two  programs  can  be  run  in  uni nterrupted  sequence. 
However , the  CURRENT  program  must  have  been  run  before  the  FIELD  program  can 
be  run  alone  because  the  latter  requires  computed  information  from  the 
former . Once  the  CURRENT  program  has  been  run  and  its  output  stored  on  the 
taoe,  the  FIELD  program  can  be  run  any  number  of  times  to  compute  the  field 
strengths  at  any  number  of  arbitrary  field  points.  Each  of  these  field 
computations  will  pertain  to  the  one  set  of  currents  computed  by  the  CURRENT 
program,  of  course. 

HPL  CURRENT  program,  Part  1,  is  stored  in  file  #1  of  the  magnetic  tape. 
CURRENT  program.  Part  2,  and  the  FIELD  program  stored  in  files  #2  and  #3, 
respectively.  Computed  currents  and  their  coordinates  from  the  CURRENT 
program  are  stored  in  file  #4. 

6.6  Print  Options  for  HPL  Programs 
The  two  HPL  programs  are  provided  with  a means  to  either  permit  or 
prevent  the  printing  out  of  certain  numerical  information  that  is  available 
within  the  computer.  This  is  accomplished  by  the  use  of  flags.  The  user  can 
determine  whether  or  not  to  print  this  information  by  modifying  the  programs 
before  running  them.  Some  of  the  printable  information  is  useful  primarily  as 
an  aid  in  diagnosing  the  cause  of  problems  and  would  not  normally  be  printed 
out  in  a routine  run.  Other  information,  such  as  the  computed  currents  J,  may 
or  may  not  be  needed  in  hard-copy  form.  Printing  out  any  information 
lengthens  the  time  required  to  run  the  programs. 

The  print  option  applies  to  the  following  information.  The  procedure  for 
setting  the  print  option  flags  is  described  below. 

1.  Flag  1 controls  the  printing  of  the  six  constants. 


a. 

OMEGA 

= 2*PI*FREQ,  where  FREQ  = test  frequency 

b. 

MU 

= permeability  of  free  space 

c. 

EPSILON 

= permittivity  of  free  space 

d. 

KAY 

= OMEGA*|/(MU*EPSILON) 

e. 

0MS4P 

= 0MEGA*/MU/(4*PI) 

f. 

S04PE 

= 1/(0MEGA*4*PI*EPSIL0N) 
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2. 


Flag  2 controls  the  printing  of  the  A coefficients. 

3.  Flag  3 controls  the  printing  of  the  quantities  u and  v defined  in 
section  6.1.  In  the  HPL  programs,  u and  v are  designated  by  MA  and  NA, 
respecti vely . 

4.  Flag  4 controls  the  printing  of  the  diagonal  of  the  matrix  [R]. 

5.  Flag  5 controls  the  printing  of  Ez,  defined  by  (90). 

6.  Flag  6 controls  the  printing  of  the  residual  of  the  least  squares 

approximation  involved  in  the  solution  of  (100). 

7.  Flag  7 controls  the  printing  of  the  currents  J. 

Items  1,  2,  3,  4,  and  6,  above,  are  primarily  useful  for  diagnostic 
purposes,  and  are  not  normally  printed  out.  Ez  (item  5)  may  be  of  interest  in 
certain  cases.  The  current  J (item  7)  may  be  printed  out  to  obtain  an  idea  of 
the  distribution  of  the  induced  currents  in  the  walls  of  the  rectangular 
box.  These  currents  are,  of  course,  the  output  of  the  CURRENT  program,  and 
they  are  input  data  to  the  FIELD  program.  Thus  they  may  be  printed  out  when 
it  is  desired  to  verify  that  the  same  set  of  values  appear  in  both  programs. 

In  the  CURRENT  program,  statements  6 through  12  of  Part  1 pertain  to 
these  flags.  In  the  FIELD  program,  statements  6,  7,  and  8 apply.  To  prevent 
printing,  the  syntax  following  the  colon  is  "sfg";  to  permit  printing,  the 
syntax  is  "cfg".  The  program  listing  given  in  section  B.l  show  the  flags  set 
to  prevent  printing  of  all  the  items  listed  above. 
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7.  NUMERICAL  RESULTS  FOR  A PRACTICAL  BUILDING 


The  result  for  a building  at  a California  site  is  hereby  presented. 
Appropriate  input  data  are  given  below. 
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D 
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Door  size 


*t 
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G 
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AX 


X 
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y0 


= building  width  = 8.763  m. 

= building  depth  = 24.232  m. 

= building  height  > 5 m.  (The  ceiling  plane  may  be  neglected). 

= 4.267  m by  3.658  m.  (el ectri cal ly  large). 

= door  symmetrical  point  = 3.048  m.  (The  door  is  not 
symmetrically  situated  on  Plane  3,  however). 

= horizontal  distance  between  the  transmitter  and  the 
door  = 9.906  m. 

= height  of  the  transmitter  above  ground  = 1.219  m. 

= input  power  to  the  transmitter  = 45  W. 

= transmitting  antenna  gain  = 0.3,  including  the  ground  effect. 

= operating  frequency  of  the  transmitter  = 148.1  MHz 
= correspond!- ng  wavelength  = 2.026  m. 

= height  above  ground  of  the  feed  point  of  the  measuring  antenna 
1.4478  m.  or  1.9050  m.  (two  antenna  heights) 

= Ay  = size  of  the  grid  = 0.6096  m.  az  = 0.4572  m. 
r distance  of  Position  (1,1)  from  Plane  2 = 0.9144  m 
I distance  of  Position  (13,1)  from  Plane  5 = 0.5334  m. 

= distance  of  Position  (1,1)  from  Plane  3 = 1.0668  m. 
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zQ  = distance  from  the  floor  plane  to  the  edge  of  the  first  grid  on 

Planes  2 and  5=0 

rn  = distance  between  the  transmitter  and  reference  antenna  = 9.684  m. 

There  are  16  rows  and  13  columns.  The  actual  measured  data  are  tabulated  as 
fol 1 ows : 


Table  4 


+x  vy  *y  VZ  *7 


Position 

in  volts 

in  deg. 

in  volts 

in  deg. 

in  volts 

in  deg 

(i.i) 

0.040 

-171.4 

0.030 

-46.3 

0.030 

107.3 

(2,1) 

0.034 

-145.6 

0.013 

-25.0 

0.090 

116.6 

(3.1) 

0.009 

-97.4 

0.022 

168.5 

0.146 

110.1 

(4,1) 

0.014 

41.7 

0.024 

-104.6 

0.132 

101.2 

(5,1) 

0.009 

68.5 

0.041 

-73.0 

0.088 

116.9 

(9,1) 

0.010 

81. 

0.014 

-116.3 

0.010 

-90.4 

(10,1) 

0.009 

95.2 

0.009 

-175.4 

0.011 

-157.5 

(11,1) 

0.009 

56.3 

0.010 

135.8 

0.015 

18.9 

(12,1) 

0.006 

39.7 

0.005 

42.4 

0.011 

3.0 

(13,1) 

0.003 

-24.5 

0.011 

-14.1 

0.012 

-158.5 
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(1.2) 

0.025 

101.8 

0. 

(2.2) 

0.044 

118.7 

0. 

(3,2) 

0.026 

87.7 

0. 

(4,2) 

0.010 

-20.5 

0. 

(5,2) 

0.006 

-131.9 

0. 

(1,2) 

0.008 

40.7 

0. 

(10,2) 

0.009 

54.7 

0. 

(11,2) 

0.007 

116.7 

0. 

(12,2) 

0.006 

139.8 

0. 

(13,2) 

0.005 

96.3 

0. 

(1,3) 

0.011 

-23.7 

0. 

(2,3) 

0.029 

39.5 

0. 

(3,3) 

0.018 

10.2 

0. 

(4,3) 

0.015 

-111.7 

0. 

(5,3) 

0.015 

-158.4 

0. 

(9,3) 

0.019 

-34.4 

0. 

(10,3) 

0.020 

-27.0 

0. 

(11,3) 

0.006 

-70.6 

0. 

(12,3) 

0.008 

-153.5 

0. 

(13,3) 

0.008 

-82.3 

0. 

-110.6  0.043  -24.0 

87.1  0.044  -3.2 

101.2  0.112  1.6 

165.4  0.131  -9.7 

-154.8  0.099  5.4 

-100.6  0.012  -133.7 

54.6  0.027  133.2 

132.1  0.024  42.4 

-35.5  0.005  22.0 

-20.6  0.020  178.1 

167.7  0.048  -112.5 

43.1  0.083  -122.9 

29.4  0.127  -111.3 

22.5  0.140  -106.9 

59.9  0.096  -112.2 

161.2  0.039  -174.4 

81.1  0.027  120.4 

130.5  0.027  -35.7 

-72.3  0.022  -103.7 

-37.6  0.025  166.7 


,030 

,034 

,056 

,025 

,023 

010 

007 

009 

002 

009 

034 

039 

047 

007 

Oil 

003 

013 

Oil 

002 

016 
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(1,4) 

0.015 

-167.6 

0. 

(2,4) 

0.031 

-74.7 

0. 

(3,4) 

0.038 

-105.4 

0. 

(4,4) 

0.026 

-175.5 

0. 

(5,4) 

0.029 

110.0 

0. 

(9.4) 

0.005 

-118.3 

0. 

(in, 4) 

0.010 

-19.3 

0. 

(11,4) 

0.014 

-4.7 

0. 

(12,4) 

0.003 

9.2 

0. 

(13,4) 

0.004 

-122.6 

0. 

(1,5) 

0.017 

95.4 

0. 

(2,5) 

0.016 

-133.1 

0. 

(3,5) 

0.023 

-177.4 

0. 

(4,5) 

0.018 

115.1 

0. 

(5,5) 

0.026 

68.7 

0. 

(9.5) 

0.027 

134.6 

0. 

(in, 5) 

0.008 

133.7 

0. 

(11,5) 

0.007 

16.4 

0. 

(12,6) 

0.012 

43.8 

0. 

(13,5) 

0.006 

62.8 

0. 

119.0 

0.046 

130.2 

-15.8 

0.098 

127.6 

-75.7 

0.120 

147.0 

-128.6 

0.114 

161.7 

69.9 

0.094 

148.9 

117.6 

0.034 

105.2 

81.7 

0.023 

35.1 

-147.8 

0.040 

-73.2 

-129.7 

0.019 

-119.0 

-50.9 

0.030 

104.5 

73.1 

0.042 

23.9 

179.5 

0.098 

32.2 

180.0 

0.114 

40.3 

150.6 

0.104 

41.5 

113.6 

0.099 

37.5 

00 

• 

0.044 

41 .3 

11.9 

0.019 

-0.8 

-115.9 

0.029 

-131.3 

-148.7 

0.018 

148.5 

-145.9 

0.045 

66.7 

041 

021 

039 

025 

007 

012 

015 

008 

010 

012 

026 

008 

034 

041 

017 

Oil 

014 

006 

017 

012 
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(1.6) 

0.022 

-14.8 

0. 

(2,6) 

0.011 

119.4 

0. 

(3,6) 

0.021 

102.1 

0. 

(4,6) 

0.020 

29.1 

0. 

(5,6) 

0.027 

-38.0 

0. 

(9,6) 

0.016 

80.7 

0. 

(10,6) 

0.008 

-125.3 

0. 

(11,6) 

0.011 

-127.4 

0. 

(12,6) 

0.004 

59.7 

0. 

(13,6) 

0.007 

47.0 

0. 

(1,7) 

0.019 

-131.2 

0. 

(2,7) 

0.013 

128.5 

0. 

(3,7) 

0.012 

30.1 

0. 

(4,7) 

0.026 

-51.7 

0. 

(5,7) 

0.028 

-89.0 

0. 

(9,7) 

0.021 

-29.5 

0. 

(10,7) 

0.012 

-115.4 

0. 

(11,7) 

0.016 

-166.2 

0. 

(12,7) 

0.006 

170.5 

0. 

(13,7) 

0.004 

56.4 

0. 

0.1 

0.061 

-92.5 

63.0 

0.103 

-76.1 

102.5 

0.126 

-70.9 

78.4 

0.112 

-76.4 

29.2 

0.094 

-67.2 

-9.6 

0.016 

-51.7 

-45.3 

0.009 

-129.5 

35.2 

0.037 

180.0 

150.0 

0.013 

120.8 

-165.8 

0.036 

3.5 

-79.9 

0.054 

163.5 

-71.5 

0.078 

-179.0 

3.0 

0.110 

178.0 

-13.1 

0.123 

175.0 

-57.2 

0.098 

176.0 

-96.1 

0.018 

-102.8 

-121.4 

0.005 

-34.0 

94.2 

0.024 

116.0 

61.3 

0.010 

-53.9 

-139.3 

0.054 

-64.9 

020 

012 

046 

056 

030 

014 

020 

008 

Oil 

015 

008 

014 

027 

043 

038 

020 

023 

013 

020 

003 
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(1,8) 

0.025 

130.4 

0. 

(2,8) 

0.018 

60.4 

0. 

(3,8) 

0.020 

-21.5 

0. 

(4,8) 

0.021 

-58.6 

0. 

(5,8) 

0.012 

140.3 

0. 

(8,8) 

0.020 

-103.0 

0. 

(10,8) 

0.017 

160.9 

0. 

(11,8) 

0.021 

132.7 

0. 

(12,8) 

0.004 

152.9 

0. 

(13,8) 

0.003 

-101.8 

0. 

(1.9) 

0.032 

31.2 

0. 

(2,9) 

0.019 

2.6 

0. 

(3,9) 

0.010 

-66.9 

0. 

(4,9) 

0.002 

-115.3 

0. 

(5,9) 

0.020 

62.6 

0. 

(9,9) 

0.026 

-158.9 

0. 

(10,9) 

0.025 

108.5 

0. 

(11,9) 

0.025 

78.5 

0. 

(12,9) 

0.013 

-48.0 

0. 

(13,9) 

0.014 

-100.2 

0. 

158.6 

0.077 

70.3 

-163.7 

0.109 

71.7 

-118.1 

0.130 

82.4 

-116.7 

0.131 

89.4 

-146.5 

0.092 

80.9 

-134.3 

0.015 

-22.6 

-166.6 

0.025 

8.6 

94.9 

0.034 

56 . 5 

19.9 

0.004 

102.9 

-38.2 

0.036 

-120.5 

105.0 

0.070 

-55.0 

129.1 

0.102 

-47.5 

172.5 

0.119 

-18.8 

146.0 

0.113 

-13.0 

142.8 

0.102 

-29.9 

137.8 

0.017 

-152.9 

126.8 

0.017 

-112.8 

52.6 

0.003 

-11.3 

-42.3 

0.029 

126.1 

-112.1 

0.028 

159.9 

,021 

,028 

,032 

,037 

,043 

020 

028 

015 

025 

018 

025 

019 

020 

038 

059 

018 

024 

013 

022 

019 
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(1,10) 

0.027 

-55.6 

0.018 

24.4 

0.069 

-136.3 

(2,10) 

0.017 

-88.4 

0.020 

7.9 

0.110 

-134.8 

(3,10) 

0.004 

145.0 

0.018 

31.5 

0.117 

-122.7 

(4,10) 

0.009 

-158.2 

0.040 

59.1 

0.105 

-118.6 

(5,10) 

0.029 

-60.9 

0.051 

63.4 

0.097 

-133.5 

(9,10) 

0.016 

139.0 

0.034 

102.2 

0.023 

171.3 

(10,10) 

0.017 

21.2 

0.023 

65.1 

0.027 

-138.4 

(11,10) 

0.019 

-26.3 

0.022 

-53.7 

0.024 

-89.6 

(12,10) 

0.013 

-96.2 

0.018 

-122.3 

0.011 

7.4 

(13,10) 

0.006 

-145.3 

0.021 

-179.0 

0.024 

88.0 

In  addition,  at  x = Dc  = 3.048  m and  y = hm  = 1.448  m (door  aperture), 

Vz  = 0.119v  and  = - 40.2  deg. 

A few  important  conclusions  may  be  drawn  from  the  above  table. 

(1)  On  the  average,  Vz  is  the  dominant  component  as  it  is  expected. 

(2)  The  values  of  Vz  are  consi sitently  larger  at  (i,  j)  with  i = 3,4  and 
5 for  all  j's,  since  these  positions  are  in  the  lit  region  directly 
illuminated  by  the  source.  The  values  of  Vz  are  relatively  weaker 
at  (i,j)  for  i = 9 through  13  near  plane  5,  since  those  positions 
are  in  the  shadow  region. 

(3)  For  i = 3,4,  and  5 (lit  region),  values  of  Vz  vary  between  a high  of 
0.146  volt  and  a low  of  0.088  volt,  indicating  a standing-wave  type 
of  field  distribution  inside  the  building. 

(4)  Phases  of  the  component  fields  change  rather  irregularly. 
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In  view  of  the  subject  presented  in  Section  3,  one  particular  component 
field  strength  or  polarization  with  respect  to  the  narrow  slot  which  may 
appear  on  an  object  to  be  placed  inside  the  building  may  have  more  penetration 
effect  than  the  other  components  or  polarizations.  For  this  reason, 
measurements  of  the  component  field  strengths  rather  than  the  total  field 
amplitude  (|EX|^  + |Ey|^  + |EZ|2)^  should  prove  more  meaningful. 

A second  set  of  measured  data  with  a new  measuring  antenna  height 
(hm  = 1.905  m)  is  tabulated  as  follows: 

Table  5 


Gri  d 

Position 

Vx 

in  volts 

in  deg. 

vy 

in  volts 

in  deg. 

vz 

in  volts 

^z 

in  deg. 

tun 

0.055 

179.0 

0.055 

-49.7 

0.059 

87.5 

(2,1) 

0.063 

-165.9 

0.037 

176.3 

0.113 

115.6 

(3,1) 

0.015 

171.6 

0.074 

176.3 

0.157 

111.9 

(4,1) 

0.028 

63.1 

0.034 

-154.6 

0.153 

101.0 

(5,1) 

0.030 

91.4 

0.046 

-61.2 

0.119 

121.4 

(9,1) 

0.025 

146.3 

0.008 

-8.5 

0.022 

-54.0 

(10,1) 

0.020 

143.8 

0.012 

-18.8 

0.016 

-152.3 

(11,1) 

0.014 

139.6 

0.017 

-99.3 

0.026 

66.5 

(12,1) 

0.010 

114.5 

0.009 

175.6 

0.017 

26.9 

(13,1) 

0.008 

62.9 

0.020 

95.2 

0.020 

-137.1 
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(1.2) 

0.037 

99.3 

(2,2) 

0.063 

118.6 

(3,2) 

0.030 

98.7 

(4,2) 

0.017 

-9.2 

(5,2) 

0.012 

-73.3 

(4,2) 

0.022 

43.1 

(10,2) 

0.013 

52.0 

(11.2) 

0.009 

155.6 

(12,2) 

0.012 

179.0 

(13,2) 

0.006 

75.1 

(1,3) 

0.030 

-17.5 

(2,3) 

0.037 

22.1 

(3,3) 

0.030 

-19.4 

(4,3) 

0.045 

-79.7 

(5,3) 

0.026 

-98.3 

(9,3) 

0.038 

-18.3 

(10,3) 

0.028 

-21.4 

(11,3) 

0.011 

-58.8 

(12,3) 

0.011 

-128.1 

(13,3) 

0.008 

-79.1 

042 

-117.1 

0.060 

-22.5 

034 

124.1 

0.083 

-8.2 

064 

114.5 

0.14,0 

-3.6 

024 

156.8 

0.152 

-5.3 

019 

-138.2 

0.124 

2.2 

003 

-41  .5 

0.027 

-137.8 

007 

-18.7 

0.035 

143.3 

010 

-114.1 

0.037 

54.7 

003 

132.4 

0.011 

5.7 

017 

61.7 

0.031 

169.0 

044 

170.5 

0.068 

-129.6 

029 

60.7 

0.096 

-124.6 

046 

21.5 

0.140 

-106.9 

025 

-8.3 

0.154 

-100.2 

010 

46.0 

0.122 

-114.5 

006 

179.8 

0.048 

-142.7 

008 

162.6 

0.030 

125.4 

015 

-173.0 

0.041 

-19.7 

008 

80.1 

0.031 

-79.1 

024 

34.7 

0.032 

177.2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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(1,4) 

0.020 

-132.7 

0. 

(2,4) 

0.035 

-64.5 

0. 

(3,4) 

0.038 

-94.9 

0. 

(4,4) 

0.035 

-164.3 

0. 

(5,4) 

0.032 

154.4 

0. 

(4,4) 

0.011 

-92.3 

0. 

(10,4) 

0.008 

-13.9 

0. 

(11,4) 

0.008 

53.2 

0. 

(12,4) 

0.007 

86.4 

0. 

(13,4) 

0.003 

-79.4 

0. 

(1,5) 

0.021 

120.2 

0. 

(2,5) 

0.011 

-135.5 

0. 

(3,5) 

0.025 

-178.3 

0. 

(4,5) 

0.046 

128.8 

0. 

(5,5) 

0.040 

107.6 

0. 

(4,5) 

0.033 

162.9 

0. 

(10,5) 

0.016 

148.9 

0. 

(11,5) 

0.012 

100.3 

0. 

(12,5) 

0.020 

86.1 

0. 

(13,5) 

0.009 

80.4 

0. 

133.0 

0.071 

134.4 

-43.3 

0.123 

138.7 

-67.0 

0.140 

155.3 

-89.0 

0.132 

166.5 

-6.0 

0.123 

159.1 

128.9 

0.039 

116.9 

157.3 

0.028 

39.5 

-152.7 

0.050 

-52.9 

-90.1 

0.024 

-103.5 

15.8 

0.036 

128.5 

96.5 

0.064 

24.2 

179.5 

0.110 

35.5 

-162.0 

0.127 

49.4 

171.1 

0.123 

50.5 

134.2 

0.123 

50.4 

72.0 

0.050 

56.5 

56.0 

0.020 

11.6 

-131.6 

0.035 

-118.3 

-101.9 

0.025 

164.1 

-31.2 

0.054 

87.1 

044 

012 

052 

034 

008 

012 

014 

022 

005 

023 

027 

015 

040 

038 

019 

022 

010 

018 

017 

Oil 
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(1.6) 

0.017 

13.1 

(2,6) 

0.010 

95.8 

(3,6) 

0.029 

82.4 

(4,6) 

0.029 

42.0 

(5,6) 

0.031 

-11.2 

(9,6) 

0.016 

110.6 

(10,6) 

0.013 

-113.6 

(11,6) 

0.011 

-142.8 

(12,6) 

0.014 

100.4 

(13,6) 

0.008 

91.7 

(1,7) 

0.020 

-109.8 

(2,7) 

0.011 

173.3 

(3,7) 

0.010 

42.7 

(4,7) 

0.027 

-39.8 

(5,7) 

0.033 

-75.1 

(9,7) 

0.020 

-18.6 

(10,7) 

0.021 

-100.2 

(11,7) 

0.019 

-127.4 

(12,7) 

0.002 

163.4 

(13,7) 

0.008 

46.5 

0.022 

37.2 

0.087 

-76.1 

0.020 

103.5 

0.118 

-70.6 

0.047 

116.6 

0.130 

-61.4 

0.048 

98.6 

0.120 

-63.1 

0.028 

50.5 

0.107 

-59.5 

0.024 

29.4 

0.020 

-53.9 

0.012 

3.5 

0.015 

-127.9 

0.009 

152.4 

0.043 

-170.3 

0.007 

-142.9 

0.015 

130.1 

0.015 

-89.4 

0.046 

17.8 

0.013 

-44.7 

0.078 

178.1 

0.014 

-14.9 

0.098 

176.1 

0.034 

32.5 

0.122 

-168.5 

0.041 

3.9 

0.130 

-167.9 

0.042 

-47.9 

0.111 

-171.4 

0.024 

-70.5 

0.019 

-91.0 

0.013 

-128.2 

0.003 

2.1 

0.023 

120.3 

0.027 

118.0 

0.015 

90.9 

0.018 

-29.4 

0.011 

-71.0 

0.055 

-53.0 

(1,8) 

0.030 

152.6 

0. 

(2,8) 

0.015 

146.7 

0. 

(3,8) 

0.014 

-84.4 

0. 

(4,8) 

0.014 

-136.6 

0. 

(5,8) 

0.024 

154.0 

0. 

(9,8) 

0.019 

-101.6 

0. 

(10,8) 

0.013 

172.6 

0. 

(11,8) 

0.013 

134.2 

0. 

(12,8) 

0.008 

-10.5 

0. 

(13,8) 

0.007 

-32.5 

0. 

(1,9) 

0.024 

55.3 

0. 

(2,9) 

0.016 

-2.6 

0. 

(3,9) 

0.009 

-58.2 

0. 

(4,9) 

0.005 

141.9 

0. 

(5,9) 

0.028 

82.7 

0. 

(9,9) 

0.019 

-153.1 

0. 

(10,9) 

0.028 

98.8 

0. 

(11,9) 

0.017 

64.6 

0. 

(12,9) 

0.017 

-77.2 

0. 

(13,9) 

0.015 

-111.1 

0. 

-149.6 

0.083 

86.0 

-127.2 

0.109 

84.7 

-79.8 

0.127 

97.1 

-86.7 

0.121 

101.9 

-115.5 

0.093 

92.1 

-109.3 

0.018 

-1.2 

-158.4 

0.027 

19.5 

93.4 

0.032 

63.9 

45.8 

0.003 

-178.8 

-61.7 

0.038 

-111.6 

162.2 

0.083 

-33.5 

-177.7 

0.112 

-31.3 

-170.5 

0.118 

-15.3 

172.2 

0.109 

-5.9 

161.8 

0.100 

-23.4 

-178.4 

0.016 

-135.9 

125.2 

0.018 

-106.3 

36.5 

0.002 

1.4 

-8.1 

0.028 

144.8 

-114.4 

0.036 

164.0 

018 

021 

028 

038 

049 

034 

022 

028 

025 

019 

014 

017 

023 

032 

041 

017 

013 

021 

018 

013 
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-126.4 


(1,10) 

0.025 

(2,10) 

0.015 

(3,10) 

0.004 

(4,10) 

0.006 

(5,10) 

0.029 

(9,10) 

0.011 

(10,10) 

0.021 

(11,10) 

0.012 

(12,10) 

0.020 

(13,10) 

0.011 

-55.6  o.mo 

-67.2  0.007 

41.8  0.012 

-16.5  0.033 

-30.1  0.041 

107.7  0.035 

45.5  0.015 

-15.8  0.022 

-120.5  0.009 

-164.6  0.021 


23.6  0.077 

3.8  0.111 

108.5  0.107 

103.5  0.079 

86.0  0.087 

118.0  0.030 

71.5  0.028 

-44.6  0.025 

-102.3  0.011 

163.7  0.028 


-127.5 

-117.1 

-114.6 

-121.9 


178.3 
-129.1 

-69.1 

37.6 

102.3 


Note  that  the  general  conclusions  made  for  the  data  presented  in  Table  4 
with  hn  = 1.448  m.  hold  good  also  for  the  data  presented  in  Table  5 with  h^  = 
1.905  m.  The  highest  value  for  Vz  is  now  0.157  volt  (against  0.146  volt  for 
Table  4)  at  position  (3,1).  This  result  is  expected  from  the  behavior  of  the 
ground  reflection  coefficient  Ry  given  in  (92).  Thus,  if  possible, 
measurements  should  be  made  with  the  feed  point  of  the  measuring  antennas 
farther  away  from  the  floor.  Of  course,  this  is  limited  by  the  supporting 
pole  and  the  convenience  of  making  measurements. 
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Based  on  the  measured  data  presented  in  Tables  4 and  5,  and  the  hybrid 
approach  represented  by  (77)  to  (79),  we  obtain  the  following  results. 

(1)  When  the  number  of  patches  to  be  used  is  eight  on  Planes  1,  2,  and  3 
represented  by  m^)  = n^  = m^)  = n^  = m^)  = n(3)  = 2,  the  minimum  number 
of  locations,  p,  where  the  measured  data  are  needed  is  also  eight  in 
accordance  with  (98)  or  (99).  Thus,  we  choose  the  measurement  positions 


(1.1).  (2.D 

, (1.2) 

and  (2,2)  with  two 

measuri nc 

1 antenna 

heights  so 

that  the 

equal  sign  in  (98)  or  (99)  holds.  The 

speci fic 

results  are: 

for  Plane  1 

(floor) 

where  z'  = 0, 

position 

(i.i) 

Jx  = 0.5750  + 

j0.3598 

Jy 

= -0.4200  - 

jO.4349 

position 

(2,1) 

Jx  = 0.4628  + 

j 1.0371 

Jy 

= -0.3880  + 

j 0.3937 

position 

(1,2) 

Jx  = 0.2582  - 

jO.3361 

Jy 

= -0.2460  - 

j 0.586 5 

position 

(2,2) 

Jx  = -0.9534  - 

j 0.5706 

Jy 

= 0.7870  + 

j 0.2328, 

(10 

for  Plane  2 

(front) 

where  x'  = 0, 

position 

(1,1) 

Jy  = 1.3713  - 

j 0.4874 

Jz 

= -0.9528  + 

,10.2119 

position 

(2,1) 

Jy  = -0.1948  + 

j 0.1866 

Jz 

= -0.6499  - 

j 0.0266 

posit i on 

(1,2) 

Jy  = -0.9331  + 

jO.1377 

Jz 

= 0.1453  + 

j 0.8761 

position 

(2,2) 

Jy  = 0.1786  + 

j 0.0755 

Jz 

= -0.2295  - 

j 0.87 28, 

(1 

and  for  Plane  3 (right)  where  y'  = 0, 

position 

(1,1) 

Jx  = -2.1408  + 

j 0.0044 

Jz 

= -1.2517  + 

jO. 3303 

position 

(2,1) 

Jx  = -0.1024  - 

j 0. 1476 

Jz 

= -0.9308  - 

j 0.3842 

position 

(1,2) 

Jx  = 1.5798  + 

j0.2727 

Jz 

= 0.0751  - 

j 1.2033 

posi ti on 

(2,2) 

Jx  = 0.0568  - 

,10.1989 

Jx 

= 0.0481  + 

j 1.2338 

(10: 

(2)  When  = n(^) 

= m^2)  = m^)  = 3 

and  n^ 

= n(3)  . 

2,  we  require  p > 14. 

Actually  we  i 

choose  the  measurement  positions  (1, 

1),  (2,1) 

, (3,1),  (1 

,2),  (2,2), 
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(3,2),  (1,3),  (2,3),  and  (3,3),  with  two  measuring 
positions.  Thus,  the  greater  sign  in  (99)  holds, 
for  Plane  1 where  z‘  = 0, 


antenna  heights,  a total  of  18 
The  1 east-square-f it  results  are: 


position 

(i,i) 

Jx 

= 0.0496 

+ 

jO.0705 

Jy 

= 

0.0705 

- 

j 0.0432 

position 

(2,1) 

Jx 

= -0.0494 

- 

j 0.1247 

Jy 

= 

0.1369 

+ 

jO.0110 

position 

(3,1) 

Jx 

= 0.0274 

+ 

j 0.0885 

Jy 

= 

-0.0898 

- 

j 0.0077 

position 

0,2) 

Jx 

= -0.1276 

- 

j 0.0836 

Jy 

= 

-0.0908 

+ 

j 0.0498 

position 

(2,2) 

Jx 

= 0.1340 

+ 

j 0.1464 

Jy 

= 

-0.1433 

- 

j 0 . 0202 

position 

(3,2) 

Jx 

= -0.0712 

- 

j 0.1124 

Jy 

= 

0.0780 

+ 

j 0 . 0154 

position 

(1,3) 

Jx 

= 0.0627 

+ 

j 0.0519 

Jy 

= 

-0.0725 

- 

j 0.06 10 

position 

(2,3) 

Jx 

= -0.0536 

- 

j 0.0808 

Jy 

= 

0.0972 

+ 

j 0.0545 

position 

(3,3) 

Jx 

= 0.0201 

+ 

j 0.0589 

Jy 

= 

-0.0353 

- 

j 0 . 033 1 , 

Plane  2 where  x1  : 

= 0, 

position 

(1,1) 

Jy 

= 0.0350 

- 

j 0.0439 

Jz 

= 

0.0174 

- 

j 0 . 0146 

position 

(2,1) 

Jy 

= -0.0316 

+ 

j 0.0370 

Jz 

= 

0.0364 

- 

j 0.0289 

position 

(3,1) 

Jy 

= 0.0447 

- 

J0.0075 

Jz 

= 

0.0275 

- 

jO.0140 

position 

(1,2) 

Jy 

= -0.0170 

+ 

j 0.0 185 

Jz 

= 

0.0033 

- 

j 0.0340 

position 

(2,2) 

Jy 

= 0.0137 

- 

j 0 . 0129 

Jz 

= 

0.0339 

+ 

j 0 . 0263 

position 

(3,2) 

Jy 

= -0.0270 

+ 

j 0.0049 

Jz 

= 

-0.0204 

- 

jO.0117, 

for  Plane 

3 where 

y'  - 

0, 

position 

(1,1) 

Jx 

= 0.0505 

- 

j 0.005 7 

Jz 

= 

0.0254 

+ 

j 0.006 7 

position 

(2,1) 

Jx 

= -0.0710 

+ 

j 0.0250 

Jz 

= 

0.0001 

- 

j 0.0 125 

position 

(3,1) 

Jx 

= 0.0321 

- 

jO.0213 

Jz 

= 

-0.0046 

4- 

j 0 . 0012 

posi ti on 

0,2) 

Jx 

= -0.0446 

- 

j 0.0084 

Jz 

= 

-0.0581 

+ 

j 0 . 0352 

position 

(2,2) 

Jx 

= 0.0638 

- 

j 0.0048 

Jz 

= 

0.0138 

- 

j 0.0600 

position 

(3,2) 

Jx 

= -0.0298 

+ 

j 0.0080 

Jz 

= 

-0.0143 

+ 

j 0 . 03 15. 

(103) 


(104) 


(105) 
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(3)  When  we  use  m^)  = = m^)  = m^)  = 4 and  n^  = n^)  = 2,  we  require 

p > 16.  Actually  we  choose  the  measurements  positions,  (1,1),  (2,1),  (3,1), 
(4,1),  (1,2),  (2,2),  (3,2),  (4,2),  (1,3),  (2,3),  (3,3),  (4,3),  (1,4),  (2,4), 
(3,4),  and  (4,4)  with  two  measuring  antenna  heights,  a total  of  32 
positions.  Thus,  the  greater  sign  in  (99)  also  holds  here.  The  specific 
results  are:  (also  in  the  sense  of  least  square-fit) 

for  Plane  1 where  z'  = 0, 


posi ti on 

(1.1) 

= 

-0.0642 

- 

JO .0120 

Jy 

= 

0.0770 

+ 

jO.0337 

position 

(2.1) 

Jx 

= 

0.1257 

- 

j 0.0 152 

Jv 

= 

-0.1334 

+ 

j 0.0348 

position 

(3,1) 

Jx 

= 

-0.0948 

+ 

jO.0174 

Jy 

= 

0.1441 

- 

jO.0577 

position 

(4,1) 

Jx 

0.0456 

+ 

j 0.0060 

Jv 

= 

-0.0852 

+ 

j 0.0102 

position 

(1,2) 

Jx 

= 

0.1129 

- 

j 0.0859 

Jy 

= 

-0.0985 

- 

j 0.0197 

position 

(2,2) 

Jx 

= 

-0.1797 

+ 

j 0.1288 

Jy 

0.1543 

- 

j 0 . 0446 

position 

(3,2) 

Jx 

= 

0.1447 

- 

j 0.0659 

Jy 

= 

-0.1803 

+ 

jO.0552 

position 

(4,2) 

Jx 

= 

-0.0770 

- 

jO. 0015 

Jy 

0.1098 

+ 

j 0.0063 

posi ti on 

(1,3) 

Jx 

= 

-0.1266 

+ 

j 0.1000 

°y 

= 

0.0750 

- 

jO.0255 

position 

(2,3) 

Jx 

= 

0.1949 

- 

j 0.1147 

Jy 

= 

-0.1379 

+ 

j 0.0678 

position 

(3,3) 

Jx 

= 

-0.1585 

+ 

j 0.0349 

Jy 

= 

0.1776 

- 

j 0.0469 

posi ti on 

(4,3) 

Jx 

= 

0.0717 

+ 

j 0.0195 

Jy 

= 

-0.1044 

- 

j 0.0243 

posi ti on 

(1,4) 

Jx 

= 

0.0913 

- 

j 0 . 03 14 

Jy 

= 

-0.0681 

+ 

jO.0279 

position 

(2,4) 

Jx 

= 

-0.1295 

+ 

j 0.0254 

Jy 

= 

0.1222 

- 

jO. 0554 

position 

(3,4) 

Jx 

= 

0.0932 

+ 

j 0.0184 

Jy 

= 

-0.1312 

+ 

j 0.0240 

position 

(4,4) 

Jx 

= 

-0.0277 

- 

j 0.0291 

Jy 

= 

0.0611 

+ 

j 0.0242, 
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for  Plane  2 where  x'  = 0 


position 

(i.i) 

Jy 

= 

0.0143 

- 

j 0 . 0057 

Jz 

= 

-0.0125 

- 

j 0.0083 

position 

(M) 

Jy 

= 

-0.0419 

- 

j 0.0099 

Jz 

= 

-0.0105 

- 

j 0.0079 

posit i on 

(3,1) 

Jy 

= 

0.0522 

+ 

j 0 . 0234 

Jz 

= 

0.0185 

+ 

j 0 . 01 49 

position 

(4,1) 

Jy 

= 

-0.0159 

- 

j 0.0256 

Jz 

-0.0123 

- 

j 0.0180 

position 

(1,2) 

Jy 

= 

-0.0009 

- 

j 0 . 0027 

Jz 

= 

0.0548 

+ 

j 0.0570 

position 

(2,2) 

Jy 

= 

0.0091 

+ 

j 0.0 144 

Jz 

= 

-0.0625 

- 

j 0.0301 

position 

(3,2) 

Jy 

= 

-0.0182 

- 

jO.0143 

Jz 

= 

0.0282 

+ 

j 0 . 0 1 7 9 

position 

(4,2) 

Jy 

= 

0.0062 

+ 

jO.0136 

Jz 

-0.0004 

- 

j 0.0042, 

and  for  Plane 

1 3 where  y 

1 

0 

position 

(i.i) 

Jx 

= 

0.0110 

- 

j 0 . 0056 

Jz 

= 

0.0650 

- 

j 0.0023 

position 

(2,1) 

Jx 

= 

0.0116 

+ 

j 0.0274 

Jz 

= 

-0.0477 

+ 

j 0.0469 

position 

(3,1) 

Jx 

= 

-0.0410 

- 

j 0.0080 

Jz 

= 

0.0291 

- 

j 0.0498 

position 

(4,1) 

Jx 

= 

0.0218 

- 

j 0.0 133 

Jz 

= 

0.0016 

+ 

j 0.0209 

posit i on 

(1,2) 

Jx 

= 

-0.0141 

- 

jO.0102 

Jz 

= 

-0.0973 

- 

j 0.07 57 

position 

(2,2) 

Jx 

= 

0.0085 

- 

jO.0141 

Jz 

= 

0.1167 

+ 

j 0.0304 

posi ti on 

(3,2) 

Jx 

= 

0.0148 

+ 

jO.OllO 

Jz 

= 

-0.0770 

- 

j 0 . 0 2 2 6 

position 

(4,2) 

Jx 

= 

-0.0083 

+ 

j 0 . 0057 

Jz 

= 

0.0036 

+ 

jO . 0169. 

Results 

for  other 

cases 

with  m(-*-) 

= n(D  = m(2) 

= m1 

(3) 

= 5,  n( 

: 2) 

=n(3)  = 

(107) 


(108) 

and 


m(D  = „(D  . m(2)  , mO)  . 6>  n(2)  . n ( 3)  . 2 are  also  available  but  not  listed 
here. 


Comparing  the  results  in  (100)  with  those  corresponding  ones  in  (103)  and 
(106),  we  find  that  JY  and  J in  Plane  1 vary  substantially  without  clear 
indications  of  convergence  to  certain  values  when  the  total  number  of  patches  used 
is  increased  up  to  4 x 4 x 2.  This  is  not  unexpected  because  the  incremental 
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interval  of  ax  = AY  = 0.6096  m (approximately  0.3x)  is  too  large.  Should  a smaller 

increment  in  the  order  of  ax  = AY  = 0.2  m (or  O.lx)  be  used,  a tendency  of 

convergence  for  JY  and  J may  be  established.  Same  comments  also  apply  to  and  J 

x y y l 

in  Plane  2 and  to  Jx  and  Jz  in  Plane  3. 

When  the  equivalent  surface  current  densities  obtained  in  (100)  through  (108), 
however  imperfect  as  they  are,  are  used  in  (77)  to  (79)  for  computing  the  field 
strengths  at  spaces  inside  the  building,  the  results  are  still  quite  interesting  as 
seen  from  the  summary  presented  below. 

With  the  current  densities  in  (100)  to  (102),  based  on  eight  patches  (2  x 2 x 

2),  substituted  in  (77)  to  (79),  we  obtain  the  field  strength  components  at  x = 

0.9144  m (xQ),  y = 1.0668  m (yQ),  and  z = 1.4478  m (hm)  as  follows: 

Ex  = 0.0714  7-171.4°  , Ey  = 0.0536  /-46.3°  , Ez  = 0.0536  /1Q7.20  V/m. 

Since  the  effective  length  of  the  dipoles  being  used  in  the  measurement  is 

approximately  a = i = % = 0.5527  m,  the  received  voltage  components  becomes 
x y z 

vx  = 0.0395  7-171.4°  , Vy  = 0.0296  7-46.3°  , Vz  = 0.0296  /1Q7.2°  v, 

which  are,  of  course,  the  same  as  those  appearing  at  position  (1,1)  in  Table  1 

because  position  (1,1)  is  one  of  the  data  points.  Similarly,  the  same  results  at 

positions  (1,2),  (2,1)  and  (2,2)  in  Table  1 and  those  in  Table  2 with  x = 0.9144  m, 

y = 1.0668  m and  z = 1.9050  m (a  second  measurement  height)  are  also  recovered. 

The  voltage  components  computed  at  x = 1.2192  m,  y = 1.3716  m and  z = 1.6764  m, 

which  are  not  one  of  the  original  data  points,  with  different  number  of  patches  are: 

(1)  With  2 x 2 x 2 or  the  results  in  (100)  through  (102), 

Vx  = 1.5958  /160.0°  , Vy  = 1.0780  7-14.8°  , Vz  = 0.0584  /129.7°  v; 

(2)  With  3 x 3 x 2 or  the  results  in  (103)  through  (105), 

vx  = °*0663  /162,7°  , Vy  = 0.0037  /27.5°  , Vz  = 0.0781  /86.2°  v; 
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(3)  With  4 x 4 x 2 or  the  results  in  (106)  through  (108), 

Vx  = 0.0531  7-175.1°  , Vy  = 0.0246  /78.9°  , Vz  = 0.0776  /83.7°  v. 

For  other  different  numbers  of  patches  such  as  5x5x2,  6x6x2  and  7 x 6 x 2 of 
which  the  current  densities  are  also  available  though  not  explicitly  listed,  the 
results  of  voltage  components  are  also  presented  below  for  the  analysis  purpose. 


Number  of 
Patches 

Vx 

in  volts 

il; 

Tx 

in  deg. 

vy 

in  volts 

*y 

in  deg. 

vz 

in  volts 

*z 

in  deg. 

5x5x2 

0.0433 

155.8 

0.0245 

125.6 

0.0811 

81.1 

6x6x2 

0.0539 

173.2 

0.0124 

-66.7 

0.0828 

90.3 

7x6x2 

0.0407 

165.1 

0.0017 

31.2 

0.0807 

84.1 

Clearly,  the  voltage  components  or  their  field  strengths  counter  parts  do  converge 
to  reasonable  levels  when  the  number  of  patches  used  in  the  computation  is 
increased.  Furthermore,  the  z component  is  still  the  dominant  one  as  it  should  be 
even  with  the  number  of  patches  as  small  as  3 x 3 x 2.  The  y-component  is  the 
weakest  as  it  should  be,  since  the  direction  of  propagation  is  along  the  y-axis. 

This  simple  analysis  indicates  the  method  proposed  in  this  report  for  computing  the 
field  strength  inside  a building  based  on  the  limited  number  of  measurement  data  may 
be  satisfactory  to  a certain  degree.  With  large  computing  facilities  available  to 
accommodate  a larger  number  of  patches,  the  accuracy  in  computation  and  the  ability 
to  predict  the  field  strength  distribution  inside  a building  similar  to  the  geometry 
shown  in  figure  9 can  be  improved. 
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8.  CONCLUDING  REMARKS 


The  problem  of  determining  the  field  distribution  inside  an  empty  building  or 
cavity  with  apertures  due  to  an  external  source  has  been  formulated  for  two 
different  configurations.  One  of  them  requires  that  the  aperture  take  the 
form  of  a horizontal  narrow  slot  contained  in  an  infinite  conducting  plane, 
that  the  cavity  wall  is  lossy,  that  the  cavity  is  in  free  space,  and  that  the 
wave  from  the  external  source  is  normally  incident  on  the  slot.  For  this 
case,  the  problem  is  solved  by  an  application  of  the  equivalence  principle  and 
the  generalized  network  presentation.  It  is  found  that  at  natural  resonances 
the  field  amplitude  in  the  slot  aperture  is  inversely  proportional  to  the 
square-root  of  the  wall  conductivity  and  that  a field  of  high  strength  can 
exist  inside  the  cavity. 

The  other  configuration  is  large  buildings  and  apertures  situated  on  a 
lossy  but  homogeneous  ground.  For  this  latter  case,  the  field  distribution 
inside  the  building  is  determined  by  an  analytical  means  together  with  a 
limited  number  of  measurement  data.  This  approach  enables  an  engineer  to  make 
a quick  on-site  check  with  a moderate  computing  facility  available,  even 
though  the  accuracy  is  not  very  good.  To  improve  the  accuracy  in  the  final 
computed  results,  a larger  computing  facility  will  be  required. 
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APPENDIX  A.  FORTRAN  PROGRAMS 
A.l  CURRENT  Program 


PROGRAM  CURRENT  ( IN  PUT  * OUTPUT* TAPE5- INPUT, T APE6-0UTPUT* TAPE7 ) 
C THIS  PROGRAM  READS  FIELD  DATA  IN  VOLTS 
C AND  COMPUTES  CURRENTS  IN  PLANES  1*  2,  AND  3 

COMPLEX  0MS4P*S0*PE*GMN1*GNL2*  6ML3*  S1*S2*S3*AMN1*ANL2* AML  3 

1 *BMN1*BNL2*BML3*CMN1*CNL2*CML3*CJX1,CJX3,CJY1,CJY2 

2 *E(24  >*CJZ2*C  JZ3,  A,DIAG(2<»)*C1*RSUBV*  ESUBZ,XS(2<») 
INTEGER  JPVT(24)*NXT0TH(2)*NYT0TH(2)*IXH(13*2)*IYH(16*2) 
COMMON  A(24*24) 

REAL  MU*  K A Y*  L SUBZ 
NPRNT1-0 


NPRNT21 

NPRNT31 

NPRNT4’ 

NPRNT5* 


C NPRNT1 
C NPRNT2 
C NPRNT3 
C NPRNT^ 
C NPRNT5 


CONTROLS 
CONTROLS 
CONTROLS 
CONTROLS 
CONTROLS 
MDIMA-24 
NXTOTHUi  -2 
NXTOTH(2)-2 
N YTOTH ( 1 ) *2 
NYTOTH<2)-2 
I XH ( 1 * 1 ) *1 
I XH ( 1 * 2 ) -1 
I X H ( 2 * 1 ) *2 
IXH<2*2»«2 
IYHC1* 1)*1 
I YH ( 1 » 2 ) * 1 
I YH  ( 2 > 1 ) ■ 2 


PRINTING  OF  0MEGA*MU*EPS*KAY*0MS4P*S04PE 
PRINTINT  OF  XP,YP,ZP  AND  COEFFICIENTS 
PRINTING  OF  M A » N A 
PRINTING  OF  DIAGONALS 
PRINTING  OF  RESIDUALS 


IYH(2*2 )»2 
DX-  • 6096 
DY-.6096 
DZ  *572 
DXDY-DX+OY 
DXDZ»DX*DZ 
D YDZ-DY*DZ 


XST  ART  91** 
YSTART  *1 #0668 
ZSTART»1.**78 
FREQ-1 81 E 8 
P I-*. * AT  AN ( 1 • ) 


OMEGA-2. *PI*FREQ 


68 


45 


50 


55 


60 


65 


70 


75 


80 


05 


90 

95 

i 


C 

c 


c 

c 

c 

c 


MU*4.E-7*PI 
EPS-1. E-9/(36.*PI) 

LS-1 

LF-2 

KAY-OMEGA*SQRT(MU*EPS ) 

0MS4P»CMPLX(0.,0MEGA*MUm.*PI ) ) 

S04PE  — CMP  LX  <0.,  .25/  (0MEGA*PI*EPS ) ) 

IF  (NPRNT1.GT.0)  WRITE  (6,15)  OME  G A , M U,  E P S , K A Y , OMS  4P»  S 04  P E 
15  F ORMAT ( *0  THE  VARIABLES  OMEGA, MU, EPS# KAY, 0MS4P , S04PE  ARE 
1LISTED  BELOW* / ( 8E15 • 6 ) ) 


IF— 2 


DO  100  IZ-LS,LF 
Z-(FL0AT(IZ)-1.)*DZ+ZSTART 
NXTQT-NXTQTH ( I Z ) 
NYTOT-NYTOTH(IZ) 

DO  100  IR -1 , NYTOT 


I Y= I YH { I R , I Z ) 

Y-(FLOAT( IY)-1. )*DY+YSTART 
DO  100  IC-1,NXT0T 
IX-IXH( IC, IZ  ) 

X«{ FLOAT ( IX)-1.)*DX*XSTART 
IF-IF+3 

IF(NPRNT2.GT.0) WRITE (6,1)  X,Y,Z  _TCI 

F OR  MAT ( /*1THE  FOLLOWING  COEFFICIENTS  CORRESPOND  TO  THE  FIEL 
1P0INT  */*  X-*F8.4*,  Y-*F8.4*,  Z**F8.4/) 


COMPUTE  COEFFICIENTS  IN  THE  FLOOR  OR  Z-0  PLANE 
(SUBSCRIPT  1,  BOTTOM  XY-PLANE,  FLOOR) 


JU  — 1 

NYT-NYT0TH(1 ) 

DO  101  IN-1, NYT 
N- I YH ( IN, 1 ) 

YP- ( FLOAT ( N )-l. )*DY*YSTART 
NXT-NXTOTH(l) 

DO  101  I M - 1 , N XT 
M-IXH( I M, 1 ) 

XP-(FL0AT(M)-1. )*DX+XSTART 


JU-JU+2 

R MN1 -SORT ( (X-XP)**2+(Y-YP)**2+Z**2 ) 

GMN1-CMPLX(C0S(KAY*RMN1),-SIN(KAY*RMN1) ) /RMN1 
Sl»CMPLX(KAY*KAY*RMNl*RMNl-3.,-3.*KAY*RMNl ) /RMN1**4 

AMN 1 - ( X— X P ) *S  1 
BMN1- ( Y— YP )*S1 
CMN1-Z*S1 

IFCNPRNT2.GT.0) WRITE (6,31)  RMN 1 ,GMN1 , SI , 

1AMN1,BMN1,CMN1  ..  ..  .,c1n 

FORMAT ( *0  RMN1-*F10.6*,  GMNl -*2 F 1 0 . 6* , S1»*2F10.6/ 
1*  AMN1-*2F 10.6* , BMN1«*2F10.6*,  C MN 1 * *2 F 1 0 . 6 ) 
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100 

105 

110 

115 

120 

125 

130 

135 

140 

145 


ooor>  ooo  o o o or>o 


COMPUTE  COEFFICIENTS  OF  JX  AND  JY  CORRESPONDING  TO  E SUB  X 

CJX1«-0MS4P*GMN1*DXDY-S04PE*(CMPLX( 1 .,KAY*RMN1 ) /RMN1**2 
1 4-(X-XP)*AMNl)*GMNl*DXDY 

CJY1— SQ4PE*(Y-YP)*AMN1*GMN1*DXDY 
A ( IF, JU) *C JX1 
A( IF,  JU+1)»C JY1 

IF(NPRNT2.GT.Q)WRITE(6,2)  XP , YP , C J X 1 , C J Y1 

2 FORM AT ( *0  COEFFICIENTS  OF  JX  AND  JY  CORRESPONDING  TO 

IE  SUB  X AT  XP»*F8.4*,  YP-+F8.4*  AND  ZP-O.  AR  E * / ( 4E 1 5 . 6 ) ) 

COMPUTE  COEFFICIENTS  OF  JX  AND  JY  CORRESPONDING  TO  E SUB  Y 
CJX1--S04PE*(X-XP)*BMN1*GMN1*DXDY 

CJY1— 0MS4P*GMN1*DXDY-S04PE*(CMPLX( 1 KAY+RMN1) /RMN1**2 
1 ♦(Y-YP)*BMN1)*GMN1*DXBY 

A(IF+1,JU) -CJX1 
A(IF+1, JU*1)«CJY1 

IF(NPRNT2.GT.0) WRITE (6,3)  XP,YP,CJX1,CJY1 

3 FOR  MAT ( *0  COEFFICIENTS  OF  JX  AND  JY  CORRESPONDING  TO 

IE  SUB  Y AT  XP-*F8.4*,  YP«*F8.4*  AND  ZP«0.  AR E*  / ( 4E1 5 . 6 ) ) 

COMPUTE  COEFFICIENTS  OF  JX  AND  JY  CORRESPONDING  TO  E SUB  Z 

CJX1»-S04PE*( X-XP )*CMN1*GMN1*DXDY 
C JY1»-S04PE*( Y-YP) *CMN1*GMN1*DXDY 
A(  IF*2,JU> -C JX1 
A(IF+2»JU*1)-CJY1 

IF(NPRNT2.GT.O) WRITE (6,4)  XP,YP,CJX1,CJY1 

4 FORMAT ( *0  COEFFICIENTS  OF  JX  AND  JY  CORRESPONDING  TO 

IE  SUB  Z AT  XP-*F8. 4* , YP«*F8.4*  AND  ZP-O.  ARE* / ( 4E 15. 6 ) / ) 
101  CONTINUE 

COMPUTE  COEFFICIENTS  IN  THE  SECOND  OR  X-0  PLANE 
(SUBSCRIPT  2,  NEAR  YZ-PL  ANE  > FRONT  WALL) 

DO  102  L*L  S, LF 

ZP- ( FLOAT ( L )-l. )*DZ*Z START 

NYT-NYTOTHCL ) 

DO  102  IN-1, NYT 
N * I YH ( I N, L ) 

YP-  (FLOAT (N)-l. )*DY*YSTART 
JU-JU+2 

RNL2-SQRT( X**2M Y-YP )**2+( Z-ZP )**2 ) 

GNL2«CMPLX(C0S(KAY*RNL2),-SIN( KAY*RNL2) ) /RNL2 

S2-CMPLX(KAY*KAY*RNL2*RNL2-3.,-3.*KAY*RNL2) /RNL2**4 

ANL2-X*S2 

BNL2-( Y-YP )*S2 

CNL2* ( Z-ZP ) *S2 

IF (NPRNT2.GT.0) WRITE (6,32)  RNL 2 , GNL 2 , S2 , 

1ANL2,BNL2,CNL2 
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150 

155 

160 

165 

170 

175 

180 

185 

190 

195 
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32  FORMAT (*0  RNL2«*F10.6*,  GNL2  **2F10 • 6* , S2«*2F10.6/ 

1*  ANL2»*2F10.6*,  BNL2»*2F10.6*,  C NL 2 * *2 F 1 0 . 6 ) 

COMPUTE  COEFFICIENTS  OF  JY  AND  JZ  CORRESPONDING  TO  E SUB  X 

C JY2»-S04PE* ( Y-YP ) *ANL2*GNL2*DYDZ 
C JZ2«-S04PE*  ( Z-ZP)  *ANL2*GNL2*DYDZ 
A{ IF,JU)*CJY2 
A(IF, JU+1 1-CJZ2 

IF (NPRNT2.GT.0) WRITE (6,5)  YP, ZP,C  JY2, C JZ2 

5 FORMAT ( *0  COEFFICIENTS  OF  JY  AND  JZ  CORRESPONDING  TO 

IE  SUB  X AT  XP*0.,  YP«*F8.4*  AND  ZP«*F8.4*  ARE*  / ( 4E 15. 6 ) ) 

COMPUTE  COEFFICIENTS  OF  JY  AND  JZ  CORRESPONDING  TO  E SUB  Y 

C JY2«-0MS4P*GNL2*DYDZ-S04PE*(CMPLX( 1 .,KAY*RNL2)/RNL2**2 
1 ♦( Y-YP  )*BNL2 )*GNL2*DYDZ 

CJZ2«-S04PE*(Z-ZP)*BNL2*GNL2*DYDZ 
A ( I F *1 , J U ) *C  J Y2 
A ( IF*1,JU+1)-CJZ2 

IF (NPRNT2.GT.0) WRITE  C6,6)  YP,ZP,CJY2,CJZ2 

6 F ORMAT { *0  COEFFICIENTS  OF  JY  AND  JZ  CORRESPONDING  TO 

IE  SUB  Y AT  XP*0.,  YP-+F8.4*  AND  ZP»*F8.4*  ARE*/(4E15.6) ) 


COMPUTE  COEFFICIENTS  OF  JY  AND  JZ  CORRESPONDING  TO  E SUB  Z 
C JY2  — S04PE*  ( Y-YP  )*CNL2*GNL  2*0 YDZ 

CJZ2*-0MS4P*GNL2*DYDZ-S04PE*(CMPLX( 1 .,KAY*RNL2  ) /RNL2**2 
1 ♦CZ-ZP)*CNL2S*GNL2*DYDZ 

A(IF*2, JU) -C JY2 
A(IF*2>JU*1)*CJZ2 

IF(NPRNT2.GT.O) WRITE (6,7)  YP,ZP,CJY2,CJZ2 
7 F OR MA T ( *0  COEFFICIENTS  OF  JY  AND  JZ  CORRESPONDING  TO 
IE  SUB  Z AT  XP-O.,  YP  *♦ F8 • 4*  AND  ZP»*F8.4*  ARE*  / < 4E 15. 6 ) / ) 
102  CONTINUE 

COMPUTE  COEFFICIENTS  FOR  THE  THIRD  OR  Y«0  PLANE 
(SUBSCRIPT  3,  NEAR  XZ-PLANE,  RIGHT  WALL) 

DO  103  L*LS, LF 

ZP« (FLOAT  CL )-l. )*DZ*ZSTART 

NXT*NXTOTH ( L ) 

DO  103  IM*1,NXT 
M*I XH( IM, L ) 

XP«  (FLOAT (M)-1.)*DX4-XSTART 
JU-JU+2 

RML3*SQRT ( (X-XP)**2*Y**2*( Z-ZP )**2 ) 
GML3«CMPLX(C0S(KAY*RML3),-SIN(KAY*RML3) ) /RML3 
S3«CMPLX(KAY*KAY*RML3*RML3-3.,-3#*KAY*RML3) /R ML  3**4 

AM  L 3 * ( X— X P ) * S 3 
BML3-Y*S3 
CML3«( Z-ZP )*S3 
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200 

205 

210 

215 

220 

225 

230 

235 

240 

245 

250 


ooo  ooo  ooo 


IFCNPRNT2.GT.0) WRITE (6,33)  R ML  3, GML 3 , S 3 , 

1 AML3, BML3 , CML3 

33  FORMAT ( *0  R Ml  3- *F  10  • 6* » GML 3-*2F 10. 6*>  S3«*2F10.6/ 

1*  AML3«*2F10.6*,  BML3-+2F10 .6*,  C ML  3 «*2F 1 0 . 6 ) 

COMPUTE  COEFFICIENTS  OF  JX  AN0  JZ  CORRESPONDING  TO  E SUB  X 

C JX3— 0MS4P*GML3*DXDZ-S04PE*(CMPLX( 1. ,KAY*RML3) /RML3**2 
1 +( X-XP)*AML3)*GML3*DXDZ 

CJZ3--S04PE*(Z-ZP)*AM13*GML3*DXDZ 
ACIF, JU)«CJX3 
A ( I F , J U+l  ) *C  J Z3 

IF(NPRNT2.GT.O) WRITE (6,8)  XP,ZP,CJX3,CJZ3 

8 FORMAT ( *0  COEFFICIENTS  OF  JX  AND  JZ  CORRESPONDING  TO 

IE  SUB  X AT  XP-4F8.4*,  YP-0.  AND  ZP-+F8.4*  AR E * / ( 4E 1 5 . 6 ) ) 

COMPUTE  COEFFICIENTS  OF  JX  AND  JZ  CORRESPONDING  TO  E SUB  Y 

CJX3— S04PE*( X-XP )*BML3*GML3*DXDZ 
C JZ3  — S04PE*( Z-ZP)*BML3*GML3*DXDZ 
A(IF*1,JU)-CJX3 
AdF+l,  JU*1)-CJZ3 

IFCNPRNT2.GT.0) WRITE (6,9)  XP,ZP,CJX3,CJZ3 

9 FORMAT ( *0  COEFFICIENTS  OF  JX  AND  JZ  CORRESPONDING  TO 

IE  SUB  Y AT  XP-4F8.4*,  YP-O.  AND  ZP-*F8.4*  AR E * / ( 4E 15 . 6 ) ) 

COMPUTE  COEFFICIENTS  OF  JX  AND  JZ  CORRESPONDING  TO  E SUB  Z 
CJX3*-SQ4PE*(X-XP)*CM13*GML3*DXDZ 

C JZ3»-0MS4P*GML3*DXDZ-S04PE*(CMPLX( 1 ,»KAY*RML3)/RML3**2 
1 ♦( Z-ZP)*CML3)*GML3*DXDZ 

A( IF  + 2, JU) -CJX3 
A ( IF *2 » JU  + 1 ) *C  JZ 3 

IF(NPRNT2.GT.O) WRITE (6 >10)  XP,ZP,CJX3,CJZ3 
10  FORM AT ( *0  COEFFICIENTS  OF  JX  AND  JZ  CORRESPONDING  TO 

IE  SUB  Z AT  XP-+F8.4*,  YP-O.  AND  ZP-+F8.4*  AR E * / ( 4E 1 5 . 6 ) / ) 
103  CONTINUE 
100  CONTINUE 
NA-JU+1 
MA-IF+2 

IF(NPRNT3.GT.0)WRITE<6,13)  MA,NA 
13  FORMAT ( *0  MA-*I5*  NA-*I5/) 

CALI  CQRDCT(MDIMA,MA,NA,A,JPVT,DIAG) 

IF (NPRNT4.GT.0) WRI TE (6,12 ) (DIAG(I),I*1,NA) 

12  F0RMAT(*0DIAG0NAl*//(8E15«6)  ) 

Cl -CMP  LX ( .376681,-. 004938) 

DT0R-PI/180. 

FTOM-. 304801 
KAY-177. 72*DT0R 
X PP-10. ♦F  TOM 
YPP--32 • 54  FTOM 
WR I TE ( 6, 17 ) 
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255 

260 

265 

270 

275 

280 

285 

290 

295 

300 
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17  FORMAT!/*  INPUT  DATA  AND  E SUB  Z FROM  THE  ANTENNA  AT  THE 
IF  IE  LD  POINTS*/ /3(3X*AMP*4X*PHASE*1X) > 1 1 X*  E SUB  Z*/) 

C READ  THE  OBSERVED  FIELD  DATA 
C 

KX-0 
K Y*  1 
K Z * 1 
K YF  *16 
IF  — 2 

DO  200  IZ*LS,LF 
Z*(FLOAT( I Z ) -1  • )*DZ+ZSTART 
IF  (IZ.EQ.2  ) GO  TO  201 
ZD-.75*FT0M 
TC*80.*DT0R 
GO  TO  202 

201  ZD«2.25*FT0M 
TC-99.+DT0R 

202  NYTOT-NYTOTHUZS 
DO  200  I N* 1 > N YTOT 
I Y*IYH(IN>IZ) 

Y»( FLOAT ( I Y ) — 1 • )*BY+Y$TART 
NXTOT-NXTOTH ( IZ  J 
DO  200  IM-1, NXTCT 
I X-I XH( IM> I Z ) 

X*(FLOAT( I X ) - 1 • )*DX+XSTART 
I F « I F ♦ 3 

217  KX«KX*1 

IF(KX.LE.IB)  GO  TO  218 
KX-1 
K Y-KY+1 

IF (KY.LE.KYF  ) GO  TO  218 
KY-1 
K Z -KZ* 1 
K YF  *1 1 

IF(KZ.GT.2)  GO  TO  219 

218  READ ( 5» 14  ) AMPX, PHSX, AMPY, PHSY, AMP Z, PHS Z 
14  FORMAT ( 6E10.3 ) 

IF((IX.NE.KX).QR.<XY.NE.KY).OR.(IZ.NE.KZ) ) GO  TO  217 

REMOVE  DIRECT  ANTENNA  CONTRIBUTION  FROM  THE  OBSERVED  FIELD  DATA 

L SUB  Z*  • 56 
AMPX*AMPX /LSUBZ 
AMPY*AMPY/LSUBZ 
AMPZ*AMPZ /LSUBZ 

R SUB T*  SORT ( ( XPP-X ) **2* ( YPP-Y ) **2  ) 

THETA-90. *DTOR-AT AN (ZD/RSUBT) 

RSUBV* ( COS ( THETA) -Cl )/ (COS ( THETA )*C1 ) 

R* SORT (R SUB T**2*ZD**2) 

ARG*KA Y*R  *TC 

ESUBZ-(1.*RSUBV)*18.515939/R*CMPLX(C0S(ARG) »- S IN ( ARG  ) ) 

WRITE (6* 16  ) AMPX,PHSX,AMPY,PHSY,AMPZ,PHSZ,ESUBZ 
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305 

310 

315 

320 

325 

330 

335 

340 

3 4 5 

350 

355 


16  F0RMAT(F8.6,F8.2»F8.6,F8.2,F8.6,F8.2,2E15.6) 

E ( I F ) -AMPX*CMPLX(C0S (DT0R*PHSX), SIN ( DT0R*PHSX  ) ) 

E(  IF  + 1) -AMPY*CMPLX(COS(DTOR*PHSY),SIN(DTOR*PHSY) ) 

E ( IF +2) «AMPZ*CMPLX ( C OS ( DTOR * PHSZ ) , S I N ( DTOR * P HS Z ) )-ESUBZ 
200  CONTINUE 

219  CALL  CQRSLT(MDIMA,MA,NA,A,JPVT,DIAG, E, XS,RESID) 
RESID-RESID/FLOAT(MA) 

IF(NPRNT5.GT.Q) WRITE<6,21)  RESIO 

21  FORMAT ( /♦  RESIDUAL-+E15.6) 

WRITE (6»18) 

18  FORMAT ( /♦  CURRENT  IN  THE  ZP«0  PLANE*/4X*XP*6X*YP*15X 
1*J  SUB  X*24X+J  SUB  Y*/) 

JU--1 

NYTOT  -NYTOTH ( 1 ) 

NXTOT-NXTOTHU  ) 

DO  301  IN-1, NYTOT 
N- 1 YH ( IN , 1 ) 

YP-CFLOAT  (N)-l*  )*DY+YSTART 
DO  301  IM»1,NXT0T 
M« I XH ( I M, 1 ) 

XP-(FL0AT(M)-1. )*DX+XSTART 
JU-JU+2 

WRI TE ( 6, 22  ) XP,YP,XS(JU),XS ( JU+1 ) 

WRXTE(7,23)  X S ( J U ) , X S ( JU  + 1 ) 

22  F0RMAT«2F8.4,4E15.6) 

23  FQRMATC 4E15.6) 

301  CONTINUE 
WRXTEl 6,19) 

19  F OR  MAT  ( /*  CURRENT  IN  THE  XP«0  PL  ANE* /4  X*  YP*6  X*  ZP*15  X 
1*J  SUB  Y*24X* J SUB  Z*/) 

DO  302  L »L  S , L F 
ZP«(FL0AT(L)-1.)*DZ+Z START 
NYTOT-N YTOTH ( L ) 

DO  302  IN-1, NYTOT 
N-I YH( IN, L ) 

YP*( FLOAT (N)-l. )*DY+YS TART 
JU-JU+2 

WRITE (6, 22)  YP,ZP,XS(JU),XS< JU+1) 

WR I TE ( 7 , 2 3 ) XS(JU),XS( JU+1) 

302  CONTINUE 
WRITE (6,20  ) 

20  FORMAT ( /*  CURRENT  IN  THE  YP-0  PLANE*/4X*XP*6X*ZP*15X 
1*J  SUB  X*2  4 X*  J SUB  Z*/) 

DO  303  L *L  S , L F 

ZP- ( FLOAT (L)-l. )*DZ+ZSTART 

NXTOT-NXTOTH ( L ) 

DO  303  IM«1,NXT0T 
M - 1 XH  ( IM  , L ) 

XP- ( FLOAT <M)-1. ) ♦DX  + XST  ART 
JU-JU+2 

WRITE (6, 22)  XP,ZP,XS(JU),XS(JU+1) 

WRITE(7,23)  XS ( JU), X S ( JU+1) 

303  CONTINUE 
STOP 

END 
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SUBROUTINE  C QR DC T ( M D IM A, M , N , A , J P VT , D I AG ) 

SUBROUTINE  C QR  DC  T COMPUTES  THE  QR  FACTORIZATION  OF  THE  MATRIX  A 
WHICH  CAN  BE  INPUT  TO  SUBROUTINE  CQRSLT  IN  ORDER  TO  COMPUTE  A 
LEAST  SQUARES  SOLUTION  TO  THE  EQUATIONS  AX-B. 

INPUT  PARAMETERS 

MDIMA  THE  FIRST  OP  ROW  DIMENSION  OF  A 
M NUMBER  OF  ROWS  IN  A 

N NUMBER  OF  COLUMNS  IN  A 

A A DOUBLY  SUBSCRIPTED  COMPLEX  ARRAY  WITH  DIMENSIONS  (MDIMA, N) 

THAT  CONTAINS  THE  MATRIX  WHOSE  OR  DECOMPOSITION  IS 
TO  BE  COMPUTED. 

OUTPUT  PARAMETERS 

A CONTAINS  THE  UPPER  TRIANGULAR  MATRIX  R AND  INFORMATION 

THAT  CAN  BE  USED  TO  EFFICIENTLY  DETERMINE  THE  MATRIX 
Q.  THIS  ARRAY  IS  INPUT  TO  SUBROUTINE  CQRSLT  WHEN  A 
LEAST  SQUARES  SOLUTION  IS  COMPUTED. 

J P VT  A SINGLY  SUBSCRIPTED  INTEGER  ARRAY  OF  LENGTH  N THAT 

CONTAINS  INTERMEDIATE  RESULTS  THAT  ARE  INPUT  TO  SUBROUTINE 
CQRSLT  IF  A LEAST  SQUARES  SOLUTION  IS  COMPUTED. 

D I AG  A SINGLY  SUBSCRIPTED  COMPLEX  ARRAY  OF  LENGTH  N THAT 

CONTAINS  INTERMEDIATE  RESULTS  THAT  ARE  INPUT  TO  SUBROUTINE 
CQRSLT  IF  A LEAST  SQUARES  SOLUTION  IS  TO  BE  COMPUTED. 

INTEGER  JPVT(N) 

COMPLEX  A(MDIMA,N),DIAGCN), AHOLD 
DO  8 1*1, N 
D I AG ( I ) *0. 

J P VT ( I ) * I 
8 CONTINUE 
DO  1 J * 1 * N 
C WRITE ( 6,753 ) J 

753  FORMATt*  J-+I5) 

MMJ«M-J*1 
EMA  XN-O. 

DO  2 K * J » N 

ENORM*VNORM(MMJ, A( J,K) ) 

I F < ENORM.LE.EMAXN)  GO  TO  2 
JPVTH-K 
EMAXN*ENORM 
2 CONTINUE 

I F ( EMAXN. EQ.O.  ) RETURN 

JPVT(J)*JPVTH 

DO  3 I *1 » M 
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AHOLD*  A ( I » J ) 

A(I»J)*A(  I, JPVTH) 

A ( I > J P VTH ) *AHOLD 

3 CONTINUE 

C COMPUTE  THE  Q VECTOR 

VNA*VN0RM(1, A( J,  J ) ) 

IF ( VNA.NE .0. ) GO  TO  5 
D I AG ( J ) “-EMA  XN 
GO  TO  6 

5 DXAG(J) «-EMAXN*A(J#J ) / VN  A 

6 A(J,J)*A(J,J )-DI AG ( J ) 
QNORM*VNORM(MMJ, A ( J, J ) > 

DO  7 I ■ J * M 

A(I# J ) *A( I, J ) /ONORM 

7 CONTINUE 
IF(J.EQ.N)  RETURN 
JP1-J+1 

DO  4 K * J P 1 » N 

CALL  QTX(MMJ#A(J»J)*A(J>K)*A(J*K)) 

4 CONTINUE 
1 CONTINUE 

RETURN 

END 


SUBROUTINE  CQRSLT(MDIMA,M,N,A,JPVT,DIAG>B,X,RESID) 

SUBROUTINE  CORSLT  COMPUTES  THE  LEAST  SQUARES  SOLUTION  OF 
AX-B  USING  THE  OR  DECOMPOSITION  THAT  IS  COMPUTED  BY 
SUBROUTINE  CQRDCT. 

INPUT  PARAMETERS 

MDIMA  THE  FIRST  OR  ROW  DIMENSION  OF  A 
M THE  NUMBER  OF  COLUMNS  OF  A. 

N THE  NUMBER  OF  COLUMNS  OF  A. 

A A TWO  DIMENSIONAL  COMPLEX  ARRAY  THAT  CONTAINS  THE 

INFORMATION  FOR  THE  OR  DECOMPOSITION  OF  A THAT  WAS 
COMPUTED  BY  SUBROUTINE  CQRDCT. 

J PVT  A SINGLY  SUBSCRIPTED  INTEGER  ARRAY  OF  LENGTH  N THAT 

CONTAINS  RESULTS  THAT  WERE  COMPUTED  BY  SUBROUTINE  CQRDCT. 

D I A G A SINGLY  SUBSCRIPTED  COMPLEX  ARRAY  OF  LENGTH  N THAT 

CONTAINS  RESULTS  THAT  WERE  COMPUTED  BY  SUBROUTINE  CQRDCT. 

B A SINGLY  SUBSCRIPTED  COMPLEX  ARRAY  OF  LENGTH  M THAT 

CONTAINS  THE  RIGHT  HAND  SIDE  OF  THE  EQUATIONS  AX*B. 
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OUTPUT  PARAMETERS 

X A SINGLY  SUBSCRIPTED  COMPLEX  ARRAY  OF  LENGTH  N THAT 

CONTAINS  THE  LEAST  SQUARES  SOLUTION  OF  THE  EQUATIONS  AX-B. 

RESID  THE  L 2 NORM  OF  THE  RESIDUAL  VECTOR  R-B-AX. 

COMPLEX  A(MDIMA,N),DIAG(N)*B(M) ,X<N)#SUM, XHOLD 
INTEGER  JPVT(N) 

DO  1 I *1  * N 
X C I ) -0. 

IF(VN0RM(1#DIAG(I  ) ) • EQ • 0 • I GO  TO  1 
NM-I 

1 CONTINUE 
DO  2 J *1 > NM 
MMJ  *M-J+1 

CALL  QTX(MMJ>A(J>J)>B(J)>B(J  ) ) 

2 CONTINUE 
LB-M-NM 
RESID*0. 

IF(LB.GT.O)  RESID* VN0RM(LB»B(NM+1) ) 

C BACKSUBSTITUTION 

X(NM)*B(NM)/DIAG(NM) 

NMM1 *NM-1 
DO  5 J*1»NMM1 
K * N M — J 
KP1*K*1 
SUM  *0 • 

DO  6 L*KP1>NM 
SUM*  SUM  + A ( K#  L )*X(L) 

6 CONTINUE 

X(K)*( B(K)-SUM) /DIAG(K) 

5 CONTINUE 

C UNSCRAMBLE  THE  SOLUTION 

DO  9 1*1*  NHM1 
J *NM— I 

J P VTH* J P VT ( J ) 

XHOLD  *X { J ) 

X ( J )*X( JPVTH) 

X( JPVTH)*XHOLD 
9 CONTINUE 
RETURN 
END 


77 


r>  o o 


FUNCTION  VNOR  M { M»  X ) 


FUNCTION  VNORM  COMPUTES  THE  NORM  OF  THE  M-VECTOR  X. 

COMPLEX  X ( M ) 

XMAX-O. 

DO  1 I * 1 * M 

XMAX*AMAX1(XMAX>ABS{ R E A L ( X{ I ) ) )>ABS(AIMAG( X ( I) ) ) ) 

1 CONTINUE 
VNORM-O. 

IFtXMAX.EQ.O. ) RETURN 
SUM«0. 

DO  2 I *1  > M 

SUM-SUM+CREAL (X(I) ) /XMAX)**2+( AIMAG( X C I ) ) /XMAX)+*2 

2 CONTINUE 

VNORM*XMAX*SQRT ( SUM ) 

RETURN 

END 


SUBROUTINE  QTX ( M# 0, X , OX) 

C COMPUTE  0 TIMES  X WHERE  0 IS  AN  ELEMENTARY  REFLECTOR  AND 

C X IS  AN  ARBITRARY  M-VECTOR. 

COMPLEX  Q(M) > X ( M ) >QX(M)»QHX 
OHX-O. 

DO  1 I *1  * M 

OHX«QHX*CONJG(Q(I ) )*X(I) 

1 CONTINUE 
QHX  *QHX  + OHX 
DO  2 I-1>M 
QX(I)«XCI)-OHX*Q(I) 

2 CONTINUE 
RETURN 
END 
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A.  2 FIELD  Program 


PROGRAM  FIELD  ( I N PUT  , GUTP  UT  > TA  p E 5 * I N P UT  » T A P E 6 » OUTP  'JT  ) 

C THIS  PROGRAM  READS  COMPUTED  CURRENTS  AS  DATA  FROM  FILE  123CXXX 
C AND  COMPUTES  FIELD  STRENGTHS  AT  A SELECTED  FIELD  POINT 

COMPLEX  0MSAP,S0APE»GMN1,GNL2»GML3»S1>S2,S3,AMN1,ANL2,AML3 

1 »BMN1,BNL2*BML3,CMN1,CNL2>CML3»CJX1,CJX3>CJY1,CJY2 

2 >E(3)»CJZ2,CJZ3,C1»RSUBV>ESUBZ,XS(2A) 

INTEGER  NXT0TH(2)>NYT0TH(2),IXH(13,2),IYH(16,2) 

REAL  M U>  K A Y 

NPRNT1*0 

NPRNT2*0 

NPRNT3=0 

C NPRNT1  CONTROLS  PRINTING  OF  OM E G A , MU > E P S > K A Y, OM S A P » S OA p E 
C NPRNT2  CONTROLS  PRINTINT  OF  XP,YP,ZP  AND  COEFFICIENTS 
C NPRNT3  CONTROLS  PRINTING  OF  M A » N A 
N XTOTH ( 1 ) =2 
N X TO  TH ( 2 ) *2 
NYTOTH(l) -2 
N YTOTH ( 2 ) * 2 
IXH(1,1)=1 
I XH(1,2) -1 
I XH ( 2 * 1 ) *2 
I XH ( 2 » 2 ) s 2 
IYH( 1» 1) *1 
I YH ( 1 f 2 ) * 1 
I YH  ( 2 » 1 ) ■ 2 
I YH ( 2 > 2 ) * 2 
X=0.9144 
Y*1 .0668 
Z«1.^78 
DX*  • 6096 
DY= .6096 
DZ*  .4572 
DXDY«DX*DY 
DXDZ=DX*DZ 
DYDZ  = D Y*D  Z 
XSTART*.91A4 
YSTART*1. 0668 
ZSTART*l.<fA78 
FREQ= 1 .481 E8 
PI  = A,*ATAN(1.  ) 

TPI  -PI  + PI 
OMEGA*TP I*FREQ 
MU=A.E-7*PI 
E P S = 1 • E~9  / ( 36  • * P I ) 
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l S * 1 
l F * 2 

KAY®0MEGA*SQRT(MU*EPS) 

0MS4P«CMPLX(0.,0MEGA*MU/(4.*PI ) ) 
S04PE=-CMPLX(0.,.25/(0MEGA*PI*EPS) ) 

IF(NPRNT1»GT.0) WRITE (6,15)  OMEGA, MU, EPS, KAY, 0MS4P,S04PE 
15  FORMAT ( *0  THE  VARIABLES  0M E G A , MU, E P S , K A Y, 0M S4 P , S 04 P E ARE 
1LISTED  BEL0W*/(8E15.6  ) ) 

READ  CURRENTS  IN  THE  ZP=0  PLANE,  J SUB  X AND  J SUB  Y 


JU  = -1 

NYTDT^NYTOTH (1  ) 

N X T0T  = N X TO  TH ( 1 ) 

DO  301  IN*1, NYT0T 
DO  301  IM=1»NXT0T 
J U = J U + 2 

READ(5,22)  XS(JU) ,XS( JU+1) 

22  FORMAT (4E15.6) 

301  CONTINUE 

READ  CURRENTS  IN  THE  XP*0  PLANE,  J SUB  Y AND  J SUB  Z 

DO  302  L * L S , L F 
NYTOT  sNYT0TH { L ) 

DO  302  IN*1, NYTOT 
J U®  J U+  2 

READ(5,2 2)  X S ( J U ) , X S ( J U+ 1 ) 

302  CONTINUE 

READ  CURRENTS  IN  THE  YP®0  PLANE,  J SUB  X AND  J SUB  Z 

DO  303  L s L S , L F 
NXT0T«NXT0TH(L  I 
DO  303  IM»1,NXT0T 
JU»  JU*-2 

READ(5,22)  XS ( JU) , XS ( JU+1 ) 

303  CONTINUE 

WRITE  CURRENTS  AND  COORDINATES  FOR  CHECK 
WRI TE (6,1 8) 

18  FORM  AT ( /*  CURRENT  IN  THE  ZP-0  PLANE*/4X*XP*6X*YP*15X 
1*J  SUB  X*2  4 X*  J SUB  Y*/) 

J U=  -1 

NYTOT *NYTOTH  ( 1 ) 

NXTOT-NXTOTHd  ) 

DO  321  IN*1, NYTOT 
N = I YH ( IN, 1) 

YP* (FLOAT (N )-l • )*DY+YSTART 
DO  321  IM«1,NXT0T 
M « I X H ( IM, 1 ) 

XP«(FL0AT(M)-1.)*DX+XSTART 
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JU=JU+2 

WRITE (6*322)  XP,YP*XS(JU)*XS(JU+1) 

322  FORMAT (2F8,A*AE15*6) 

321  CONTINUE 

WRI TE (6*19) 

19  FORMAT(/*  CURRENT  IN  THE  X P *0  PLANE*  / A X* YP *6 X* Z P*1 5 X 
1*J  SUB  Y*  2 A X * J SUB  Z*/) 

DO  323  L*  L S * L F 
ZP*(FLOAT(L ) -1 • )*DZ+ZSTART 
N YTOT  *NYTOTH ( l ) 

DO  323  IN-1  * NYTOT 
N*IYH( IN* L ) 

YP  « (FLOAT  ( N ) -1 . )*DY-fYSTART 
JU-JU+2 

WRITE (6*322)  YP,ZP*XS(JU)*XS(JU+1) 

323  CONTINUE 
WRI TE (6,20) 

20  FORMAT ( / * CURRENT  IN  THE  YP*0  P L A N E * / A X * X P * 6 X * Z P * 1 5 X 
1*J  SUB  X*  2 A X * J SUB  Z*/) 

DO  3 2 A L*LS,LF 

ZP«(FLOAT (L ) -1 • )*DZ+ZSTART 

NXTOT  -NXTOTH ( L ) 

DO  32A  I M = 1 , NXTOT 
M * I XH ( IM, L ) 

XP*(FL0AT(M)-1« )*DX+XSTART 
JU-JU+2 

WRITE(6,322)  XP,ZP*XSIJU),XS(JU+1) 

3 2 A CONTINUE 
I F = 1 

E ( IF  ) *0. 

E ( IF  + 1) =0. 

E ( I F +2 ) *0  . 

IF(NPRNT2.GT.O) WRITE (6*1)  X,Y,Z 
1 FORMAT ( /+1THE  FOLLOWING  COEFFICIENTS  CORRESPOND  TO  THE 
1 F I F L D POINT  */*  X = *F8.A*»  Y=*F9.A*,  Z = *F8.A/) 

C 

C COMPUTE  COEFFICIENTS  IN  THE  FLOOR  OR  Z-0  PLANE 
C (SUBSCRIPT  1*  BOTTOM  XY-PLANE*  FLOOR) 

C 

JU  = -1 

N YT -NYTOTH ( 1 ) 

DO  101  IN*1*NYT 
N = I YH ( IN, 1 ) 

YP=(FL0AT(N)-1.)*DY+YSTART 
NXT»NXT0TH(1  ) 

DO  101  IM=1,NXT 
M = I XH ( IM, 1 ) 

XP«(FL0AT(M)-1.)*DX+XSTART 

JU-JU+2 

R MN1 s SOP  T ( ( X-XP ) **2+ ( Y-YP ) **2+Z**2 ) 

GMN1=CMPLX (C0S(KAY*RMN1) , — S I N ( KAY*RMN1 ) ) / RMN1 
Sl-CMPL  X( KA Y*K AY*RMNl*RMNl-3  .,-3.*KAY*RMN1)/RMN1**-, 
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AMN1«(  X-  X P ) * S 1 
BMN1=(Y-YP)*S1 
CMN1=Z*S1 

COMPUTE  COEFFICIENTS  OF  JX  AND  JY  C 0 P R E S PON D I N G TO  E SUB  X 

CJX1=-0MSAP*GMN1*DXDY-S04PE*(CMPLX(1.,KAY*RMN1)/RMN1**2 
1 +( X-XP ) *AMN1 )*GMN1*DXDY 

CJY1=-S04PE*(Y-YP)*AMN1*GMN1*DXDY 
E(IF)=E(IF)+CJX1*XS(JU)+CJY1*XS(JU+1) 
IF(NPRNT2.GT.0)WRITE(6,2)  XP,YP»CJX1,CJY1 
? F0RMAT(*0  COEFFICIENTS  OFF  JX  AND  JY  CORRESPONDING  TO 
IE  SUB  X AT  XP»*F8.4*,  YP  = *F8.4*  AND  ZP-0.  AR E * / ( A El  5 . 6 ) ) 

COMPUTE  COEFFICIENTS  OF  JX  AND  JY  CORRESPONDING  TO  E SUB  Y 

CJX1*-S0APE*(X-XP)*BMN1*GMN1*DXDY 

CJY1=-0MS4P*GMN1*DXDY-S04PE*(CMPIX(1.,KAY*RMN1)/RMN1**2 
1 ♦ ( Y-Y5 ) ARMN1 ) *GMN1*DXD Y 

E(IF+1)«E(IF+1)+CJX1*XS(JU)+CJY1*XS(JU+1) 

IF ( NPR NT2.GT.0) WRITE (6,3 ) XP,YP,CJX1,CJY1 

3 FORMAT ( $ 0 COEFFICIENTS  OF  JX  AND  JY  CORRESPONDING  TO 

IE  SUB  Y AT  XP-+F8.4*,  YP=*F8.4*  AND  ZP  = 0.  AR E*  / ( 4E15 .6  ) ) 

COMPUTE  COEFFICIENTS  OF  JX  AND  JY  CORRESPONDING  TO  E SUB  Z 

CJX1»-S0^PE*(X-XPJ*CMN1*GMN1*DXDY 
C JY1--S04RE* ( Y-YP)*CMN1*GMN1*DXDY 
E ( I F ♦ 2 ) * E ( IF  + 2)+CJXl*XS(JU)+CJYl*XS(  J U 1 ) 
IF(NPRNT2.GT.0)WRITE(6,^)  XP,YP,CJX1,CJY1 

4 F OR  M AT ( *0  COEFFICIENTS  OF  JX  AND  JY  CORRESPONDING  TO 

IE  SUB  Z AT  XP-+F8.**,  YP-*F8.A*  AND  ZP-O.  ARE  * /(4E 15. 6 ) / > 
101  CONTINUE 

COMPUTE  COEFFICIENTS  IN  THE  SECOND  OR  X«0  PLANE 
(SUBSCRIPT  2 , NEAR  Y Z-PL  ANE , FRONT  WALL) 

DO  102  L = L S , L F 
ZP»(FL0AT(L)-1.)*DZ+Z START 
NYT«NYTOTH(L  ) 

DO  102  IN-1, NYT 
N - I YH ( IN, L ) 

YP=  ( FLOAT  (N  )-l  . )*DY4-YSTART 
JU-JU+2 

RNL2-S0RT( X * *2+( Y-YP )**2+ ( Z-ZP )**2 ) 
GNL?«CMPLX(C0S(KAY*RNL2),-SIN(KAY*RNL2))/RNL2 
S2=CMPLX( KAY*KAY*RNL2*RNL2-3.,-3.*KAY*RNL2) /RNL2**^ 

A NL  2 * X*  S2 

BN L 2 * ( Y- Y P ) * S 2 

CNL2- ( Z-ZP )*S2 

COMPUTE  COEFFICIENTS  OF  JY  AND  JZ  CORRESPONDING  TO  E SUB  X 
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CJY2«-S04PE*(Y-YP)*ANL2*GNL2*DYDZ 
CJZ2=-S04PE*(Z“ZP)*ANL2*GNL2*DYDZ 
E(IF)»E(IF)+CJY2*XS(JU)+CJZ2*XS(JU+1) 
IF(NPRNT2.GT.0)WRITE(6,5)  YP,ZP>CJY2>CJZ2 

5 FORMAT ( *0  COEFFICIENTS  OF  JY  AND  JZ  CORRESPONDING  TO 

IE  SUB  X AT  XP*Q.,  YP=*F8.4*  AND  ZP  = *F8.4*  AR  E * / ( 4 El  5 . 6 ) ) 

C 

C COMPUTE  COEFFICIENTS  OF  JY  AND  JZ  CORRESPONDING  TO  E SUB  Y 
C 

CJY2=-0MS4P*GNL2*DYDZ-S04PE*(CMPLX(1.,KAY*RNL2)/RNL?**2 
1 +(Y-YP)*BN12 )*GNL2*DYDZ 

CJZ2«-S04PE*(Z-ZP)*BNL2*GNL2*DYDZ 
E ( IF  + 1) =E ( IF  41 )+C JY2*XS( JU) +CJZ2*XS ( JU  + 1) 
IF(NPRNT2.GT.0)WRITEC6»6)  YP,ZP,CJY2,CJZ2 

6 FORMAT (*0  COEFFICIENTS  OF  JY  AND  JZ  CORRESPONDING  TO 

IE  SUB  Y AT  XP*0«,  YPS*F8.4*  AND  ZP«*F8.4*  A P E * / ( 4E 1 5 . 6 ) ) 

C 

C COMPUTE  COEFFICIENTS  OF  JY  AND  JZ  CORRESPONDING  TO  E SUB  Z 
C 

CJY2»-S04PEM  Y-YP ) *CNL2*GNL2*DYDZ 

C JZ2=-0MS4P*GNL2*DYDZ-S04PE* (CMPLX( l.,KAY*RNL2)/RNL2**2 
1 ♦ ( Z-ZP)*CN12 )*GNL2*DYDZ 

E(IF+2)»E(IF+2)+CJY2*XS(JU)+CJZ2*XS(JU+l) 
IF(NPRNT2.GT.0)WRXTEC6*7)  YP,ZP,CJY2,CJZ2 

7 FORMAT ( *0  COEFFICIENTS  OF  JY  AND  JZ  CORRESPONDING  TO 

IE  SU8  Z AT  X P ®Q  » f YP»*F8.4*  AND  ZP=*F8.4*  ARE*  / ( 4E 1 5 . 6 ) / ) 

102  CONTINUE 
C 

C COMPUTE  COEFFICIENTS  FOR  THE  THIRD  OR  Y*0  PLANE 
C (SUBSCRIPT  3 » NEAR  XZ-PIANE,  RIGHT  WALL) 

C 

DO  103  L«LS»LF 
ZP=(FL0AT(L)-1.)*DZ+ZSTART 
N XT  *NXTOTH ( L ) 

DO  103  IMS1»NXT 
M*IXH(  IM,L ) 

XPS (FLOAT ( M ) -1 • )*DX+XSTART 
JU= JU+2 

RML3*S0RT(  (X~XP!**2+Y**2+(Z-ZP)**2  ) 
GML3«CMPLX(C0S(KAY*RML3)»-SIN(KAY*RML3))/RML3 

S3«CMPLX(KAY*KAY*RMl3*RML3-3.»-3.*KAY*RML3)/RHL3**4 

A M L 3 = ( X-X  P ) * S3 
B ML  3*  Y+  S3 
CML3«( Z~ZP  )*S3 

C COMPUTE  COEFFICIENTS  OF  JX  AND  JZ  CORRESPONDING  Tr  E SUB  x 

CJX3*-0MS4P*GML3*DXDZ-S04PE*(CMPLX(1  .#kay*rwl?) 

l 4(X-XP)*AML3)*GML3*DXDZ 

CJZ3*-S04PE* ( Z—ZP)*AML3*GML3*DXDZ 

E(IF)«E(IF)+CJX3*XS(JU)+CJZ3*XS(JU+1) 
IF(NPRNT2.GT.0)WRITE (6>8)  XP  >ZP»CJX?»CJ -3 
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255 

260 

265 

270 

275 

280 

285 

290 

295 

300 

305 


ooo  ooo  ooo  ooo 


8 F OR  MAT ( *0  COEFFICIENTS  OF  JX  AND  JZ  CORRESPONDING  TO 

IE  SUB  X AT  XP«*F8.4*,  YP*0.  AND  ZP»*F8.4*  A R E * / ( 4E 1 5 . 6 ) ) 

COMPUTE  COEFFICIENTS  OF  JX  AND  JZ  CORRESPONDING  TO  E SUB  Y 

CJX3=-S04PE*(X-XP)*BML3*GML3*DXDZ 
CJZ3=-S04PE*(Z-ZP)*BMt3*GMt3*DXDZ 
E(IF+1)*E(IF+1)+CJX3*XS(JU)+CJZ3*XS(JU+1) 
IF(NPRNT2.GT.0)WPITE(6»9)  XP,ZP»CJX3,CJZ3 

9 FORMAT (+0  COEFFICIENTS  OF  JX  AND  JZ  CORRESPONDING  TO 

IE  SUB  Y AT  XP«*F8.4*,  YP-0.  AND  ZP«*F8.4*  AR E* / ( 4E 1 5 . 6 ) ) 

COMPUTE  COEFFICIENTS  OF  JX  AND  JZ  CORRESPONDING  TO  E SUB  Z 
CJX3»-S04PE*(X-XP)*CML3*GML3*DXDZ 

C JZ3=~0MS4P*GML3*DXDZ-S04PE* (CMPLX(1.,KAY*RML3)/RML3**2 
1 ♦ ( Z-ZP)*CML3  )*GML3*DXDZ 

E(IF+2)=E(IF+2)+CJX3*XS(JU)+CJZ3*XS(JU+l) 
IF(NPRNT2.GT.O)WRITE(6,10)  XP,ZP,CJX3,CJZ3 
10  F0RMAT(*0  COEFFICIENTS  OF  JX  AND  JZ  CORRESPONDING  TO 

IE  SUB  Z AT  XP«*F8.4*,  YP«0.  AND  ZP«*F8.4*  AR E * / ( 4E 1 5 . 6 ) / ) 
103  CONTINUE 
100  CONTINUE 
N A * J U+  1 
M A= I F + 2 

IF(NPPNT3.GT.0)WRITE(6,13)  MA,NA 
13  FORMAT ( *0  MA«*I5*  NA-+I5/) 

ADD  DIRECT  ANTENNA  CONTRIBUTION  TO  E SUB  Z 

j: 

Cl*CMPLX(.376681»-.004938) 

DTOR=PI /180. 

F TOM* .304  801 
K AY*177. 72*DT0R 
X P P = 1 0 . * F T 0 M 
YPP=-32. 5*FT0M 
ZD=7-4.*FT0M 

TC»(41.5573*Z+19.8333)*DT0R 
RSUBT»SORT ( ( XPP-X)**2+(YPP-Y)**2 ) 
THETA=90.*DT0R-ATAN(ZD/RSUBT) 

R SUBV= (COS (THETA) -Cl)  / (COS (THETA )+Cl ) 

R ■ S QR  T ( RSUBT**2  + ZD**2  ) 

arg=kay*r+tc 

ESUBZ=(1»+RSUBV)*18«515939/R*CMPLX(C0S(ARG)>— SIN(ARG)) 
E(IF+2)«E(IF+2)+ESUBZ 

WRITE  E SUB  X,  E SUB  Y,  AND  E SUB  Z AT  THE  FIELD  POINT 
WR  I TE  ( 6 » 1 7 ) 

17  FORMAT ( f * E SUB  X,  E SUB  Y,  AND  E SUB  Z AT  THE  FIELD 
1P0INT*//*X*7X*Y*7X*Z*3X,3(5X*AMP*4X*PHASE*1X)/) 

AMPX  = C ABS ( E ( IF  ) ) 

AMPY  = CABS( E ( I F ♦ 1 ) ) 
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AMPZ=C ABS ( E ( I F + 2 ) ) 

PHSX=ATAN2 ( AIMAG< E( IF ) ),REAL (E (IF) ) ) /DTOR 
PHSY=ATAN2(AIMAG(E(IF+1)),REAL(E(IF+1)))/DT0R 
PMSZ=ATAN2(AIMAG(E(IF+2)),REAL(E(IF+2)))/DT0R 
WRITE ( 6, 16  ) X,Y,Z,AMPX,PHSX,AMPY,PHSY>AMPZ,PHSZ 
16  F0PMAT(3F8.4,F10.6»F8.2,F10.6,F8.2,F10.6.F8.2) 
STOP 
END 


APPENDIX  B.  HPL  PROGRAMS 
B.l  CURRENT  and  FIELD  Programs 


U : 
1 : 


4 : 


10 
1 1 
1 2 
1 3 

14 

15 
1 6 

17 

18 

19 

20 
21 


43 

44 

45 

46 

47 

48 

49 

50 


E u i 1 d i n g C h ar  ac  t e r i z at,  i o n F'  r o g r am  a " : 

E M T h e o r y d e u e 1 o p e d b y D r . M . ’ T . M a " : 

H F'  L P i - o g r a rn  a-  a-  u p e r i a-  e d b y M . G . Fl  r t h u r 
1 H at  i o n a 1 E u r e-  -au  o f S t an  d an  da-"  : 

Boulder,  CO  80303": 

CURRENT 


'FLAG  1 ( 
'FLAG  2 C 
FLAG  3 i: 


'FLAG  4 
"FLAG  5 
"FLAG  6 
"FLAG  7 
rad 

24; Ml n 
d i m A [ M 


0 N T R 0 L 3 P R I N T N G 0 F C 0 N 3 T A N T S " : a f g 
ONTROLS  PRINTING  OF  COEFFICIENTS": 
0 N T R 0 L S P R INTI H G 0 F M A A H D H A " : s f g 
PR  INTI NG  OF  D I AGONALS " : sf q 


1 

sf  g 


CONTROLS 
CONTROLS 
ONTROLS 
ONTROLS 


PRINTING  OF  E SUE  Z":sfg  5 
PR  I NT  I NG  OF  RES  I DUAL " : sf g 6 
P R INTI  N G Ci  F C U R R E N T S " : a f q 7 


d i m 
spc 
prt 
f xd 


N , 2 ] , E [ M , 2 ] , D C N , 2 ] , .T  [ N ] , I [ 30  ] 
W i C 1 ] 


C C 5 0 , 2 ] , R [ 6 0 ] , XCN, 21 


"CURRENT 

4 


PROGRAM" ; spc  2 


24  > I C 1 3 ] 
“ENTER"  : 


prt 

" C 0 0 R D I H A T E S 0 F 

FIELD  POINT 

3 : " 

24 

- P 

b e e p 

; ent 

"ENTER  XI. 

XI  = ? 

" , R [ 

31  ] 

“i  cr 

ii.  J 

b e e p> 

; e n t 

"ENTER  X2. 

X 2 = '? 

" , R [ 

38  3 

2 6 

b e e pi 

: e n t. 

" E H T E R V 1 . 

V 1 = ? 

")R[ 

34] 

2 7 

b e €'  p 

; e n t 

"ENTER  V 2 . 

Y2  = ? 

" , R [ 

39] 

cL  O 

b e e p 

; ent 

"ENTER  Zl. 

Z 1 = ? 

" , R[ 

37] 

2 9 

beep 

; ent 

"ENTER  Z2 . 

Z k‘  — ’? 

" , R [ 4 0 ] 

30 

"PRINT  ENTRIES": 

31 

prt 

ll  ^ — 11 

, RC31  ] , " X2= " 

, F:  [ 3 8 ] , 

" Y 1 = 

" , R [ 3 4 ] 

32 

prt 

II  17  ^ = ii 

, F [371  , " Z2= " 

, R [ 4 0 ] 

34 

-•  P 
dap 

"ARE 

THESE  NUMBER 

S C 0 F:  R E C T '?  ( 

Y N ) " 

35 

ent 

" " , W* 

3 6 

if  U*[l]= 

"Y"  or  M T [ 1 ] 

= " y " ; g t 

o + 8 

O ~7 

b e e p 

; €•  n t. 

"RE-ENTER  XI 

XI 

= ? " 

, R C 3 1 ] 

3 8 

beep 

; e n t 

"RE-ENTER  X2 

. X 2 

= ? " 

, R [ 3 8 ] 

3 9 

beep 

; e n t 

"RE-ENTER  VI 

Y 1 

= ? " 

, R [ 3 4 ] 

4 0 

b e e p 

; e n t 

"RE-ENTER  Y2 

. Y 2 

1 1 

, R [ 3 9 ] 

41 

beep 

: e n t 

"RE-ENTER  Zl 

Zl 

_ r-i  1 1 

, R [ 3 7 ] 

42 

hi  e e pi 

; a n t 

"RE-ENTER  Z2 

. Z2 

= ? " 

, R [ 4 0 ] 

gto  "PRINT  ENTRIES" 
"COEFFICIENTS"  : 
RC38I-RC31 ] > R [ 6 ] 

R [ 3 9 ] - R C 3 4 ] > R [ 7 I 
R [ 4 0 ] - R C 3 7 ] } R [ 8 3 
R [ 6 ] * R [ 7 ] } R [ 9 ] 

R [ 6 ] *R [ 8 ] }RC 10] 

R C 7 ] *R [ 8 ] > R C 1 1 ] 
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51:  1 . 481e8>R[ 14] 

52:  4 * at  n < 1 ) } P 
53s  2*F*RC14]>RC1S] 

54:  4e -7*P } R [ 173 

55:  1 e - 9 C36*P)>R[12] 

56:  1 > I C 15]  > I [ IS  3 > I [ 1 1 ] 

57:  2> I [ 14] > I C 17] > I [ 10] 

5S : R C 1 8 ] *\ < R C 1 7 ] *R [ 1 2 ] > > R [ l e ] 

59  : 0 > C C22,  1 ] ; R [ 1 8 ] *R  [ 1 7 ] ■"  <4*P  > > C [ 22 , 2 ] 

60  : 0 > C [ 24  , 1 ] ; - . 25,-  C R [ 1 8 ] *P*R  [ l 2 ] > } C C 24 , 2 ] 

61:  if  f 1 g 1 ; g t o +6 

62 : prt  “THE  VRR I RULES  OMEGA , MU , EPS  I LONKRY , GMS4P , S04PE " 

63:  prt  " RRE  LISTED  EEL0W";flt  6 
64:  s p c 

65  : prt  R C 1 8 ] , R C 1 7 ] , R [ 1 2 ] , R C 1 6 ] , C [ 22 , 1 ] , C [ 22 , 2 ] , C C 24 , 1 ] , C [ 24 , 2 ] 

66:  spc  2 

67:  dsp  "THIS  CALCULATION  TAKES  S MINUTES" 

68:  -2  > F 

69:  for  K = I C 1 1 ] to  I C 1 0 ] 

70:  < K — 1 !'  * R [ 8 ] + R [ 3 7 ] } Z 

71:  for  J = I [ 1 8 ] to  I [ 1 0 ] 

72:  < J-l >*RC7]+RC34] >Y 
73:  for  I = I C 1 5 ] to  I [14] 

74:  < I - 1 > R [ 6 ] + R [ 3 1 ] > X 
75:  F + 3 > F 

.76:  if  f 1 g 2 ; g t-  o + 7 

78:  prt  “THE  FOLLOWING  COEFFICIENTS  CORRESPOND  TO  THE  FIELD  POINT 
79:  spc 
80:  fxd  4 

8 1 : prt  " X = " , X , " Y= " , Y , " Z= " , Z 

82:  spc  2 


83: 

-1  > 

LI 

84: 

for 

N = I C 1 8 ] to  I Cl 

—> 

\ 

] 

85: 

< N- 

1 >*RC7]+RC34] >R 

[ 

32] 

86 : 

f or- 

M  = I C 1 5 ] to  IC1 

4 

] 

87: 

al - 

1 >*RC6]+RC 3 1 ] >R 

[ 

29] 

88: 

U + 2 

>U 

89: 

N.  ( ( 

X — R C 2 9 ] > 2 + < Y - R 

[ 

32] 

> 

'■x  ii! 

+ Z ■" 

i zL 

) } R [ 

4 

] 

9 0 : 

C Q :E 

<RC16]*RC24]  ')  s R 

[ 

24  ] 

> 

C [ 

20, 

1 

] ; -s 

i 

n 

( 

P 

[ 

i 

b 

] 

* 

R 

[ 

2 

4 ] 

R 

[ 

2 

4 

] 

> C 

[20,2] 

91  : 

■;  < f 

C 1 6 ] * R [ 2 4 ] ) A 2 - 3 

j 

R [ 

2 

4] 

4 } 

c 

[25, 

1 

] 

, 

- 

■ J 

* 

R 

[ 

1 

y 

] 

ir 

R [ 

2 

4] 

R 

[ 

2 

4] 

4 > C [ 2 5 , 2 ] 

92 : 

( X — 

R [ 2 9 ] j * C C 2 5 , 1 ] > 

c 

[2, 

1 

^ ; 

< X - 

R 

[29] 

j 

* 

c 

[ 

'cL 

5 

9 

2 

] 

> 

c 

[ 

9 

2 

] 

93 : 

C Y — 

R [32] > * C [25,  1 ] > 

c 

[5, 

1 

]; 

< Y - 

R 

[32] 

> 

c 

[ 

cL 

5 

j 

2 

] 

> 

L 

[ 

c- 

9 

2 

] 

94: 

z+c 

[25,  1 ] >C[  14,  1 ] ; 

7 

*CC 

2 

5, 

2]  > 

c 

[14, 

cL 

] 

95 : 

c 1 1 

C M U L x ( - C [ 2 2 , 1 

3 

, -c 

[ 

lL  izL 

, 2] 

, 

C [ 2 0 

9 

1 

] 

9 

c 

[ 

2' 

0 

5 

2 

] 

9 

CL 

i** 

, 1 

] 

, 

c 

[ 

■? 

2]  > 

9 6 : 

C [ 7 

, 1 ] *R [ 9 ] > C [ 7 , 1 ] 

5 

C [ 7 

J 

2] 

*R[ 

9 

] > C [ 

"7 

!» 

2 

1 

97: 

1 ” R 

[ 2 4 ] -'•■  2 + C X - F:  [ 2 9 ] 

> 

*C  [ 

2 

, 1 

] > C 

[ 

28,  1 

3 

98 : 

R [ 1 

6 ] * F‘  C 2 4 ] s R [ 2 4 ] 

2 

+ ( X 

- 

R [ 

29  ] 

> 

* C [ 2 

9 

2 

] 

> 

c 

[ 

2 

y 

9 

2 

] 

99 : 

c 1 1 

'CMUL' <C[24, 1 ] 

9 

C [ 2 

4 

J *=- 

3 , c 

[ 

28,1 

3 

* 

c 

[ 

cl' 

o 

2 

] 

j 

c 

[ 

2 8 

j 

1 ] 

, 

c 

[ 

2 

2]  > 

1 0 0 

c 1 

1 'CMUL'<C[28,1 

3 

, C[ 

9 

2], 

c 

[20, 

1 

] 

9 

c 

[ 

2 

0 

2 

] 

, 

c 

[ 2 

8 

, 1 

] 

, 

c 

[ 

2 y 

« 2 ^ 1 

101 

CC 

28,  1 ] * R [ 9 ] > C [ 2 8 

9 

1 ] ; 

c 

[ 2 

8 , 2 

] 

* R [ 9 

3 

> 

c 

[ 

2 

8 

9 

2 

] 

102 

C[ 

7,  1 ]-C[23,  1 ] >C[ 

{ 

, 1 3 

9 

CL 

i'"  » L 

] 

- C [ 2 

y 

9 

2 

] 

■«. 

c 

[ 

r’ 

•I 

] 

103 

— c 

[24,1 ] * ( Y - R [ 3 2 ] 

\ 

> C [ 

Q 

, 1 

] ; - 

c 

[24, 

cl 

3 

* 

< 

v' 

- 

R 

[ 

3 

2 

] 

> 

> C 

[ 

9 

2 

] 

1 04 

c 1 

1 ' OMUL ' < C [ 9 , 1 ] 

9 

C [ 9 

9 

2] 

, C [ 

2 

,n, 

c 

[ 

2 

, 

2 

] 

j 

c 

[ 

9 

J 

1 

] , 

C 

C9 

J 

2 

] 

} 

105 

c I 

1 'OMUL" <C[9, 1 ] 

n 

C [ 9 

9 

2 ] 

, C[ 

2 

0,1] 

9 

C 

[ 

2 

0 

, 

2 

] 

c 

[ 

9 

, 1 

] 

, c 

[ 

9 

, 

2 

] > 

106 

CC 

9 , 1 ] * R [ 9 ] / L [ 9 , 1 

] 

; C [ 

9 

9 

]*R 

[ 

9 ] > C 

c 

9 

, 

2 

3 

107 

CC 

7 , 1 ] > A [ F , IJ  , 1 ] ; C 

[ 

“7 

‘ !•  *— 

] 

> A 

[F, 

IJ 

,21 

108 

CC 

9,  1 ] >R[F,  IJ+1  , 1 ] 

* 

C [ 9 

9 

2] 

> AC 

F 

, U+ 1 

9 

2 

] 
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109 

110 
1 1 1 
1 12 
1 1 3 
114 
1 15 
116 
1 17 
1 1 8 
1 19 
120 
121 

123 

124 

125 

126 

127 

128 
129 
1 30 
131 

133 

134 

135 

136 

137 

138 


i f 
Pr 

pr 

pr 

sp 

pr 

®P 
- 1 


f 1 g 2 j g t o + 8 

c 

t "COEFFICIENTS  OF  JX  FIND  JV" 
t " C 0 R R E S P 0 N D I N G T 0 E S U E X A T " ; f x d 4 
t "XF=" , RC293 , "YP=" , RE  32  3 , "AND  2P=0  ARE" ; f 1 t 
c 

t C E 7 , 1 ] , C [ 7 , 2 ] , C [ 9 , 1 ] , C C 9 , 2 3 

r 


C 1 1 
C 1 1 


4 

, 1 ] 

* ( x - 

R [ 2 

Q 

3 

) 

} 

c 

[ 

r 

J 

1 

3 

c 

[2 

4 

9 

2 

3 

( 

y 

- 

R 

[ 

9 

3 

) } 

C[ 

f 

23 

CMC 

l <:  c 

[7, 

l 

3 

c 

[ 

“7 

9 

2 

3 

n 

C[ 

5 

, 1 

3 

•i 

c 

[ 

5 

9 

4 

3 

c 

[ 

\ 

9 

1 3 

, c 

[7 

,23 

j 

CNL 

L •"  < C 

[7, 

l 

3 

•I 

C 

[ 

"7 

j 

2 

3 

9 

C[ 

4 

0, 

1 

3 

, 

c 

[ 

2 

0 

, 

2 

3 

9 

c 

L 

*'  |l 

1 3 

, c 

[7, 

2 3 > 

1 

I * R 

[ 9 3 > 

C [ 7 

9 

1 

3 

% 

c 

[ 

r 

j 

2 

3 

* F: 

[ 

9 3 

> 

c 

[ 

r> 

9 

2 

3 

CNL 

L •-  < - 

C [ 2 

2 

, 

1 

] 

, 

- 

c 

[ 

2 

2 

9 

3 

, c 

[ 

2 

0 

9 

1 

3 

, 

c 

[ 

2 

0 

9 

2 

3 , 

C[ 

Q 
-•  , 

1 3 , 

C [ 9 

23 

1 

3*R 

[ 9 3 > 

C [ 9 

9 

1 

3 

5 

c 

[ 

9 

, 

Ci 

3 

*R 

[ 

93 

> 

c 

[ 

9 

9 

2 

3 

2 

4 ] •- 

2 + < Y 

-R  [ 

4' 

3 

\\ 

* 

c 

[ 

Cj 

9 

1 

3 > 

C 

[2 

8 

5 

1 

3 

] 

*R[ 

24  3/ 

R [ 2 

4 

3 

2 

+ 

f 

V 

- 

R 

[ 

3 2 

3 

|i  ir 

c 

[ 

cr 

9 

2 

3 

> 

c 

[ 

2 

C; 

9 

2 

3 

CNL 

l •" «:  c 

[24 

? 

1 

3 

n 

c 

[ 

2 

4 

, 

2 

3 , 

c 

[2 

y 

l 

3 

!» 

c 

[ 

4 

O 

,, 

2 

3 

j 

C[ 

2 y 

, 1 

3 , C 

[28 

23 

CNL 

L ■- 1 C 

[28 

9 

1 

3 

? 

c 

[ 

2 

8 

9 

2 

3 , 

c 

[2 

0 

n 

l 

3 

9 

c 

[ 

2 

0 

, 

2 

3 

? 

C [ 

2 y 

, 1 

3 , C 

[28 

23 

9 

13* 

R [ 9 3 

> C [ 

2 

C* 

9 

1 

3 

5 

c 

[ 

2 

o 

, 2 

3 

* R 

[ 

9 

3 

/ 

c 

[ 

2 

8 

9 

2 

3 

1 

3-C 

[28, 

1 3 > 

c 

[ 

9 

, 

1 

3 

!» 

c 

[ 

9 

9 

3 

-c 

[ 

2 

y 

j 

2 

3 

> 

c 

L 

9 

9 

2 

3 

C [ 7 , 1 3 } A [ F + 1 , U , 1 3 ; C [ 7 , 2 3 } A [ F + 1 , U , 2 3 
C [ 9 , 1 ] > fi  [ F+  1 , IJ+  1 , 1 3 ; C [ 9 , 2 3 > A [ F+ 1 , U+ 1 , 2 3 
l f 
sp 
pr 
pr 

Pr 
sp 

pr 
-P 


f 1 g 2 ; g t o + 8 

t '"COEFFICIENTS  OF  JX  AND  JV" 

t " C 0 R R E S P 0 N D I N G T 0 E S U E Y AT"; f x d 4 

t " X P = " , R [ 2 9 3 , " V P = " , R [3  2 3,  "AH  D ZP  = 0 ARE"; fit  6 

c 

t C C 7 , 1 ] , C [ 7 , 2 ] , C [ 9 , 1 ] , C C 9 , 2 1 


139 

-C[ 

2 

4,  1 3*<X-R 

[29 

3 

j 

C [ 

r' 

, 1 

3 ; - 

C [ 

24, 2 

3 

* 

-R 

[ 

cL 

9 

3 

;i 

} c 

[7 

9 

3 

140 

c 1 1 

•"CNUL-'  <C[ 

“7 

r q 

1 

3 

, c 

[7 

n 

23 

, C[ 

14 

,13, 

c 

[ 

14 

9 

3 

9 

c 

[ 

■J7 

, 1 

3 , 

C [ 

r 

j c- 

3 > 

141 

c 1 1 

C M U L <C[ 

II 

1 

3 

, c 

[7 

!i 

2 3 

, C [ 

2 0 

,13, 

c 

[ 

2 0 

9 ^ 

3 

9 

c 

[ 

i-' 

, 1 

3 , 

C[ 

■J7 

9 c- 

3 > 

142 

C [ 7 

9 

1 3 * F.:  [ 9 3 >C 

[7, 

1 

3 ; 

C[ 

r 

9 

3*R 

[ 9 

3 > C [ 

r 

23 

143 

-C  [ 

2 

4,  1 3*<Y-R 

[:] 

2 

3 

CL 

9 

, 1 
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c 19, 

1 

] , C 

C 

19 

9 

cL 

] , 

c 

C 2 

9 

1 

] 

, c 

C 

23, 

2] 

> 

293- 

c 1 

1 

CMUL' 

c c c 

2 

- n 

1 

] 

, c 

c 

n ^ 

] , 

c 

[24, 

1 

],c 

C 

24 

9 

ii! 

] , 

c 

C 2 

9 

1 

] 

, C 

C 

2 8 , 

2] 

> 

294  : 

CC 

•->  Q 

, 1 ]*RC 10] 

> 

CC 

2 

Q 

, 1 

] 

; c c 2 

O q 

ii! 

] * R C 

1 

0 ] > 

c 

C 2 

8 

9 

2] 

295: 

CC 

12 

, 1 ] - C C 2 8 , 

1 

] > 

c 

c 

12 

9 

1 ] ; C 

C 1 

L 

,2]- 

c 

C 28 

, 

2] 

c 

[ 1 

k 

9 

] 

296: 

CC 

8, 

1 ] > ACF+2, 

L 

, 1 

] 

<1 

CC 

y 

, 2 ] > 

RC 

F 

+ 2,1 

, 

2] 

297: 

CC 

1 2 

, 1 ] > A 

:f  + 2 

9 

u+ 

1 

9 

1 ] 

9 

c C 1 2 

, ■=- 

] 

> R C F 

+ 

2,L 

+ 

1, 

~> 

] 

298:  if  f 1 g 2 ; g t o +9 


? 8 9 : 

spc 

3 0 0 : 

pr t "COEFFICIENTS  OF  JN  AND 

J Z " 

301  : 

prt  "CORRESPONDING  TOE  SUE 

Z AT"; 

302: 

prt  "XP=" , RC29] , " YP=0,  AND" 

303: 

prt  " Z P = " , R C 3 6 ] , "ARE" ; f 1 1 6 

304 : 

spc 

3 0 5 : 

p r t CCS,  1 ] , C C 8 , 2 ] , C C 1 2 , 1 ] , L 

C 12, 2] 

306: 

spc 

307 : 

n ext-  M 

308 : 

n e x t L 

309 : 

n ext-  I 

310: 

n e x t.  J 

311: 

n e x t K 

312: 

U+ 1 > I C 1 6 ] 

313: 

F+2> I C 12D 

314: 

i f f 1 g 3 ; g t o + 3 

315: 

spc  2 

316: 

fxd  0: prt  " MR= " , I C 1 2 ] , " NR  = " 

, I C 1 6 ] 

317: 

c 1 1 CQRDCT  < I C 1 3 ] , I C 1 2 ] , I C 1 6 ] > 

318: 

f 1 t 6 

319: 

if  f 1 g 4 ; g t o +7 

3 2 0 : 

s P c 

321  : 

prt  "DIAGONAL" 

• 

s p c 

: 2 3 : 

f o r 1=1  t o I C 1 6 ] 

324  : 

prt  DC  I , 1 ] , DC  I , 2] 

“ r - * 

n e x t.  I 

326: 

. 3 7 6 6 8 1 > C C 1 6 , 1 ] ; - . 0 0 4 9 3 8 > C C 1 6 , 2 ] 

327: 

F'  x 1 8 0 > R C 5 ] 

3 2 8 : 

. 30480 1 > R C 15] 

3 2 9 : 

177. 7 2 * R C 5 ] > R C 16] 

330 : 

RC 15]*10>RC30] 

331  : 

RC 15] *— 32. 5 > R C 3 3 ] 

332 : 

ldf  2 

333 : 

e n d 

334: 

"CQRDCT" : 

335 : 

I >p4; J > p 5 ; K>p6 

336 : 

f o r 1 = 1 t o p 3 

337 : 

0 > D C 1 , 1 ] > D C I , 2 ] ; I > J C I ] 

O O p . 

next  I 

339: 

f o r J =1  t o p 3 

340 : 

p 2 - • J + 1 > p 7 
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341 

O > p 1 2 

34  2 

f o r K = J t.  o p 3 

343 

-■  VHCiRMR  •'  < p7 , J 

, K > > 

p 1 3 

34  4 

if  p 1 3 > p 1 2 ; K > 

P 8 ; p 

1 3 > p 1 

345 

n €'  x t K 

346 

if'  p 1 2 = 0 ; g t-  g 

"CQRDCTXT 

II 

34  7 

p 8 > J [ J ] 

3 4 8 

f or  1 = 1 to  p 2 

3 4 9 

R C I , J , 1 ] > p 1 0 ; 

R [ I , 

J , 2 ] > 

pll 

350 

R [ I , p 8 , 1 ] > R C I 

, J , 1 ] 5 fl  [ I 

5 P 8 , 2 ] 

351 

p 1 0 } fi  [ I , p 8 , 1 ] 

5 pl  1 

> R C I , 

P 8,2] 

' -J  S 

n e x t I 

•*;.  cr  -z> 

•"  V NOR  MR"'  < 1 , J, 

J > > P 

16 

354 

i f p 16#  0 ; g t c> 

" 5 ,i 

355 

-pi  2 3 DC J,  i ] ; 0 

> D [ J 

,2] 

3 5 6 

gt G "6" 

357 

“5" : -pl2*fi[ J, 

,T  , 1 ] 

p 1 6 .!• 

D [ J , 1 ] 

" 6 “ : R [ J , J , 1 ] - 

D [ J , 

1 ] > fl  [ 

J , J , 1 ] 

359 

■’  V N O R M R ' < p 7 , J 

, J > > 

P 1 4 

360 

f g r I = J t,  g p 2 

361 

RC  I , J,  1 ] p 1 4 > 

fit  I , 

J , l ] ; 

R [ I , J , 

3 6 2 

n e x t I 

363 

i f J = p 3 ; g % o " 

CQRDCTXT" 

364 

J+l > p 9 

365 

f g r K = p 9 t g p 

‘l! 

366 

c 1 1 Q T X R R R •"  < 

p7,  J 

, K > 

367 

n e x t K 

368 

n e x t J 

3 6 9 

"CQRDCTXT " : 

3 7 0 

p 4 5 1;  p 5 > J ; p 6 > 

K 

37  1 

r e t 

3 ? 3 
374 


37  7 

z>  r o 

3 7 '3 
3 8 0 
381 

3 82 

383 

384 

i-  cr 

i“i  ._i 

j 8 6 

3 8 7 


" OMUL  " : 

P 1 P 3 ~ p 2 p 4 > p 7 ; p 1 p 4 + p 2 p 3 > p 8 ; p 7 > p 5 ; p 8 > p 6 ; r e t 
" C D I V " : 

p 3 p 3 + p 4 p 4 } p 7 ; < p 1 p 3 + p 2 p 4 > 7 p 7 > p 1 0 

<.  p 2 p 3 - p 1 p 4 > p 7 > p 1 1 ; p 1 0 > p 5 ; p 1 1 > p 6 ; r e t 

" QTXRRR " : 

0 > p 1 0 > p 1 1 

f o r Q = p 2 t o p 2 + p 1 — 1 
R C Q , p 2 , 1 ] J'  p 1 2 ; - H [ Q , p 2,2]  ) p 1 3 

c 1 1 ■"  C M U L ■"  ( p 1 2 , p 1 3 , R [ Q , p 3 , 1 ] , R [ Q , p 3 , 2 3 , p 1 2 , p 1 3 > 

p 1 0 + p 1 2 .!■  p 1 0 ; p 1 1 + p 1 3 > p 1 1 
n €•  x t Q 

2p 1 0 > p 1 0 ; 2 p 1 1 > p 1 1 
f o r Q = p 2 t o p 2 + p 1 - 1 

c 1 1 •"  OMUL  < p 1 0 , p 1 1 , R [ Q , p2  , 1 ] , R [ Q , p2 , 2 ] , p 1 2 , p 1 3 > 

R [ Q , p 3 , 1 ] - p 1 2 } R [ Q , p 3 , 1 ] ; R[ Q , p 3,2]- p 1 3 > RCQ, p 3 , 2 ] 
next  Q 


389 

r e t 

3 90 

" VHORMfl " : 

391 

p 2 > p 5 } p 6 

392 

0 > p 7 

3 9 3 

f g r V = 1 t 

394 

max < p7  , at 

3 9 5 

p 5 + 1 > p 5 

3 9 6 

n e x t V 

39  7 

if  p 7 = 0 ; r 

398 

0>p8 
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399 : 

for  V = 1 to  p 1 

400 : 

p 8 + 1 Fl  C p 6 , p 3 , 1 ] 

4 01: 

pt'+l  }p6 

402 : 

n e x t V 

403: 

ret  p 7 £ \ p 8 

* 159 

0 : " B u i 1 di  ng  C h ar  ac  t e r i z at  l o n P r o q r .am  s " : 

1 : "EM  T h e o r y d e v e 1 o p e d b y D r . M . T . M a " : 

2 : " H P L P r o g r am  £ £ u p e r v i s e d b y M . G . R r t h u r " : 

3 : "Nat  i o n a 1 B u r e au  o f 8 1 an d ar  d s " : 

4 : " B o u 1 d €'  r , Col  o r ad  o " : 

5:  "CURRENT  Program:  Part  2,  File  2": 

6:  rad 
7:  2 4 > M > N 
8:  -2  > F 

9:  for  K=ICU]  to  I [10] 


10 

<K-1>*R[8]+RC37] >Z 

1 1 

if  K=2;gto  "201" 

12 

. 7 5 * R [ 15] >R[35] 

13 

8 0 * R [ 5 ] > R [ 2 8 ] 

14 

q t o " 2 0 2 " 

15 

"201" : 2. 25*R[ 15] >R[35 

] 

16 

9 9 * R [ 5 ] > R [ 2 8 ] 

17 

"202" : for  J= I [IS]  to 

I [ 1 

7] 

18 

< J-l >*R[7]+R[34]  >V 

19 

for  I = I C 1 5 ] to  I [14] 

20 

<1-1 >*RC6]+RE31 ] >X 

21 

F + 3 > F 

22 

" 2 1 8 " : f x d 5 

23 

£pc  2 

24 

prt  "EM  FIELD  DhTh  RTF  I ELD  POINT 

C o • i 

25 

- P 

26 

"ENTER" : 

27 

beep ; ent  "ENTER  R M P X . 

R M P X = ? 

" , R [ 1 ] 

28 

beep ; ent  "ENTER  PHSX. 

PHSX  = ? 

" , R [ 1 9 ] 

29 

b e e p ; e n t " E H TER  R M p V . 

RMPY  = ? 

" ? R [ 2 ] 

30 

beep; ent  "ENTER  PHSY. 

PHSY  = ? 

" , R [ 2 0 ] 

31 

beep; ent  "ENTER  RMPZ. 

RMPZ  = ? 

" , R [ 3 ] 

32 

beep; ent  "ENTER  PHSZ. 

PHSZ  = ? 

" , R [ 2 1 ] 

33 

"PRINT  ENTRIES": 

34 

prt  " R M P X = " , R[ 1 ] , " PHS 

•./  _ li 

, R [ 1 9 ] , " 

R M P Y = " , R [ 

35 

prt  " RMPZ= " , R C 3 ] , "PHS 

7 = 11 

, R [ 2 1 ] 

36 

spc  2 

37 

dsp  "ARE  THESE  NUMBER 

s c 

OR R EOT? 

< Y/N ) " 

38 

ent  " " ,W$ 

39 

if  W$Cl]s"Y"  or  WtCl] 

= 11 

" ; g t o +8 

40 

beep; ent  "RE-ENTER  RMPX. 

RMPX 

= ? " , R [ 1 ] 

P H 8 V = " , RC20] 
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41  : 

b € e p ; e n t " R E - E H TER  PH 8 X . 

P H S X = ? 

" , R [ 1 9 ] 

42: 

b € e p ; e n t " R E - E H TER  A M P V . 

AMF'Y  = ? 

" , R [ 2 ] 

43: 

b e e p ; e n t " R E - E N TER  P H S V . 

PHSY  = ? 

" , R [ 2 0 ] 

44: 

beepjent  "RE-ENTER  flMPZ. 

AMPZ  = ? 

" , R [ 3 ] 

45: 

beepjent  "RE-ENTER  PHSZ. 

PHSZ  = ? 

" , R [ 2 1 ] 

46: 

gto  "PRINT  ENTRIES" 

47: 

"CALCULATION" : 

48: 

R [ 1 I.-’.  5 6 > R [ 1 ] 

4 9 : 

RC2  ]••-'.  5 6 > R C 2 ] 

50: 

R [ 3 ] •■■■' . 5 6 > R C 3 ] 

51  : 

(.  R C 3 0 ] - X ) 2 + (.  R [ 3 3 ] - Y ) 2 ) 

> R [ 2 6 ] 

Cj  . 

9 0 * R [ 5 ] - at n < R [35] / R [ 2 6 ] ) > T 

■!'  I 

c os < T > - C [16,  1 ] > C [ 23 , 1 ] ; -C [ 1 6 , 2 ] > C [ 2 

3,2] 

54: 

c o s < T > + C [ 16,1] > C [28,  1 ] ; C [ 1 

6,2] > C [ 2 8 

, 2] 

cr  tr  . 

.j  -J  ■ 

c 1 1 ■'■CIUV-'  ( C [23,  1 ] , C [23,2] 

, C [ 23 , 1 ] , 

C [28,2] 

, C [23,  1 ] , C [23, 

56: 

< R [ 2 6 ] 2 + R [ 3 5 ] 2 > >R 

57 : 

R [ 1 6 ] * R + R [ 2 8 ] > R [ 4 ] 

5 • 

C 1 +C [ 23 , 1 ] > * 1 8 . 5 15939  > C [ 1 8 

, 1 ] 

59: 

C [ 23 , 2 ] * 1 8 . 5 1 5939 > C [18,2] 

6 0 : 

C[ 18,  1 ] / R > C [ 18,  1 ] ; C[  18, 21/ 

R } C [18,2] 

61  : 

c os < R [ 4 ] > > C C 28 , 1 ] ; -si n < R [ 4 ] > > C [28,2 

] 

6 2 : 

c 1 1 C M U L <C[  18,  1 ] , C [ 1 8 , 2 ] 

, C [ 23 , 1 ] , 

C [28,2] 

, C [ 1 8 , 1 ] , C [ 1 3 , 

6 3 : 

f X d 6 

64: 

i f f 1 g 5 ; g t o + 5 

65: 

spc 

6 6 : 

prt  "E  SUE  Z = " ; f 1 t 6 

6 7 : 

prt  C[  18,  1 ] , C[ 18, 2] 

6 8 : 
69 : 

-■  F-1 

R [ 1 ] * c os  < R [ 5 ] *R [ 1 9 ] > > E [ F , 1 ] ; R [ 1 ] *s i 

n < R [ 5 ] * 

R [ 1 9 ] > E [ F , 2 ] 

70 : 

R [ 2 ] *c os ( R [ 5 ] *R [ 20 ] ) > E [ F+ 1 

, 1 ] ; R [ 2 ] * 

s i n < R [ 5 

]*R[£0]  > E [ F + 1 

71  : 

R [ 3 ] *c  os  < R [ 5 ] *R [ 2 1 ] >-C [18, 

1 ] > E [ F + 2 , 

1 ] 

—>  •”«  . 
Ida 

R [ 3 ] *s i n < R [ 5 ] *R [ 2 1 ] > -C [18, 

2 ] > E [ F + 2 , 

2] 

“7  >•  a 

n e x t I 

74: 

n €■  x t J 

-?  cr  . 

1 a 

n e x t K 

76: 

dsp  "CALCULATION  IN  PROGRE 

SS  " 

i'O 

7 9 
80 
81 

8 2 
8 3 


" 2 1 9| : c 1 1 ' CQRSLT ' (I C 1 3 ] , I C 1 2 3 , I C 1 6 ] > 

R [ 2 £ ] I C 12  ] >RC22]  ; fit  6 
i f f 1 g 6 ; g t,  o + 3 

SRC 

prt  " R E S I D U fl  L = " , RC22] 

- 1 > U 

for  W = I [ 1 8 ] to  I [17] 


84: 

<N  - 1 

> * R [ 7 ] + R [ 3 4 ] > R [ 3 

2] 

35: 

f o r 

M = I [ 1 5 ] to  I [14] 

86 : 

01-1 

> * R [ 6 ] + R [ 3 1 ] > R [ 2 

9] 

87 : 

IJ  + 2 > 

IJ 

88: 

i f f 

1 g 7 ; g t o +6 

89: 

spc 

ci 

90 : 

prt 

"CURRENT  IN  THE 

*7 

P=0  PLANE 

HT"  ; 

91  : 

prt 

" X F'  = " , R[29]  , " Y P = 

9 

R [ 32 ] j f 1 t 

6 

92 : 

spc 

9 3 : 

prt 

" J SUE  X",X[IJ,U 

[ IJ  , 2 ] , " J 

SUE  Y 

94: 

F'l  6 X t 

M 

95: 

n e x t 

N 

96 : 

for 

L=I[U]  to  1110] 

97: 

<L-1 

> * R [ 8 ] + R [ 3 7 ] } R [ 3 

61 

94 


98: 

99 : 

10O 

101 

102 

103 

1 04 

105 

1 06 

107 

1 08 

109 

1 10 

1 1 1 

112 

113 

1 14 

115 

1 16 

1 17 

118 

1 19 

120 

121 

122 

123 

124 

125 

126 

127 

129 

130 

131 

132 

133 

134 

135 

136 

137 

1 39 

140 

141 

142 

143 

144 

145 

146 

147 

148 

1 49 

150 

151 

152 

153 

154 

1.55 


for  N = I [ 1 8 ] to  I [ 1 7 3 
< N- 1 > *R  [7 ] +R  C 34 ] > R [ 32  3 
U + 2 > U 

if  f 1 g 7 ; g t-  o +6 
spc  2 

prt  "CURRENT  IN  THE  XP =0  PLANE  AT" ;fxd  4 

prt  " V P = " , R [ 3 2 3 , " Z P = 11  , R [ 3 6 3 ; f 1 t 6 

spc 

prt.  1 J SUB  V " , X C LI , 1 ] , X C IJ  , 2 3 , " .J  SUB  Z " , X [ IJ+  1 , 1 ] , X C IJ+  1 , 2 ] 
next-  N 
next  L 

for  L= I [ 1 1 ] to  I [103 
( L - 1 )*R[83+R[373  > R [ 3 6 3 
for  M = I [ 15  3 to  I [143 
<M-1 )*R[63+R[31 3 >R[293 
IJ  + 2 > Li 

if  f 1 g 7 ; g t o +6 


spc 

2 

prt 

"CURRENT  IN  THE 

Y R = 0 PLANE 

HT  " ; f 

xd  4 

prt 

" XP= " , R[293 , " ZR= " 

, R [ 3 6 3 ; f 1 t 

6 

spc 

prt 

"J  SUB  X",X[ U,13, 

X [ U , 2 3 , " J 

S U B Z " 

, XUJ+1  ,13, 

XUJ+1 ,23 

next  M 
n e x t L 


r-cf  4,R[*3,X[*3 
1 dp  3 
e n d 

"CQRSLT"  : 

I > p 4 ; J > p 5 ; K>p6; L > p 7 
f o r 1 = 1 t o p 3 
0 ; X [I,13>  X [1,23 
i f V N 0 R M D < 1 , I > # 0 ; I > p 1 3 
n e x t I 

f o r J = 1 t o p 1 3 
P 2 - J + 1 > p 1 4 

c 1 1 Q T X A B E p 1 4 , J . J ) 

n e x t J 
p2-p 1 3 } p 1 5 
0 > R [22  3 

i f p 1 5 > 0 ; ■-  V N 0 R M E C p 1 5 , p 13+1  > > F:  [223 

c 1 1 •-  CD  I V < B [ p 1 3 , 1 3 , B [ p 1 3 , 2 3 , D [ p 1 3 , 1 3 , D [ p 1 3 , 2 3 , X [ p 1 3 , 1 3 , X [ p 1 3 , 2 3 

p 1 3- 1 > p 1 6 

f o r J = 1 t o p 1 6 

p 1 3 - J > K 

K + 1 } p 8 

0 > p 2 0 > p 2 1 

f o r L = p 8 t o p 1 3 

c 1 1 ■ CMIJL  ■-  <:  A [ K , L , 1 3 , A [ K , L , 2 3 , X [ L , 1 3 , X [ L , 2 3 , p22 , p23  > 

p 2 0 + p 2 2 / p 2 0 

p 2 1 + p 2 3 > p 2 1 
next  L 

B [ K , 1 3 - p 2 0 > p 2 2 ; B [ K , 2 3 - p 2 1 > p 2 3 

c 1 1 ■-  CD  I V •-  < p 2 2 , p 2 3 , D [ K , 1 3 , D [ K , 2 3 , X [ K , 1 3 , X [ K , 2 3 > 

n e x t J 

f o r 1 = 1 t o p 1 6 
pi 3- I > J 
■J  [ J 3 > p 1 2 

X [ J , 1 3 > p 1 0 ; X [ J , 2 3 > p 1 1 
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156:  X C p 1 2 , 1 ] > X [ J , 1 ] ; X [ p 1 2 , 2 ] > X C J , 2 ] 

157:  pi  0 > X C p 1 2 , 1 ] 5 p 1 1 > X [ p 1 2 , 2 ] 

158:  n e x t I 

159:  p 4 > I ; p 5 > J ; p 6 > K ; p 7 > L 

160:  r e t » 

161:  "OMUL": 

16  2:  p 1 p 3 - p 2 p 4 > p 7 ; p 1 p 4 + p 2 p 3 > p 8 ; p 7 > p 5 ; p 8 > p 6 ; r e t 

163:  "CD  IV": 

164:  p 3 p 3 + p 4 p 4 > p 7 ; < p 1 p 3 + p 2 p 4 > ••••'  p 7 > p 1 0 
165:  ( p 2 p 3 - p 1 p 4 > / p 7 > p 1 1 ; p 1 0 > p 5 ; p 1 1 > p 6 ; r * t 
166:  “ QTXflEE  " : 

167:  0 > p 1 0 > p 1 1 

1 6 8 : f o r Q = p 2 t o p 2 + p 1 - 1 

169:  iR  C Q , p 3 , 1 ] > p 1 2 ; - R C Q , p 3,2]  > p 1 3 

170:  ell  C i'1  LI  L ''  ( p 1 2 , p 1 3 , B [ Q , 1 ] , B C Q , 2 ] , p 1 2 , p 1 3 > 

171:  p 1 0 + p 1 2 > p 1 0 ; p 1 1 +p 1 3 > p 1 1 
172:  ne x t Q 

173:  2p 1 0 > p 1 0 ; 2 p 1 1 > p 1 1 
174:  f o r Q = p 3 t a p 3 + p 1 - 1 

175:  c 1 1 •"  C M IJ  L C p 1 0 , p 1 1 , H C Q , p 3 , 1 ] , fi  [ Q , p 3 , 2 D , p 1 2 , p 1 

176:  B C Q , 1 ] -p 1 2 > E C Q , 1 ] ; E C 0 , 2 ] -p 1 3 > E C Q , 2 ] 

177:  ne x t Q 
178:  ret 
179:  "VNORMB": 

180:  p 2 > p 5 > p 6 

181:  0 > p 7 

182:  for  V=1  to  pi 

1 8 3 : m ax  < p 7 , ab  s < B C p 5 , 1 ] > , ab  s ( B C p 5 , 2 ] > > > p 7 
184:  p5+ 1 > p5 

185:  next  V 
1 8 6 : i f p 7 = 0; r e t 0 

18  7:  0 > p 8 

188:  f o r V = 1 t o p 1 

189:  p 8 + ( B C p 6 , 1 ] p 7 ) 2 + ( I'  C p 6 , 2 ] p 7 ) 2 / p 8 

190:  p 6 + 1 > p 6 

191:  next  V 

192:  ret  p 7 * \ p 8 

193:  “VNORMB": 

194:  p 2 } p 5 > p 6 
195:  6 > p 7 

196:  f o r V =1  t o fj  1 

19  7:  m ax ( p 7 , abs ( D C p 5 , 1 ] > , ab s ( D [ p 5 , 2 ] > > > p 7 
198:  p5+l>p5 

199:  n e x t V 

200:  i f p 7=0; r e t 0 

201:  0 } p 8 

282:  f o r V =1  t o p 1 

2 0 3 : p8+  D [ p 6 , 1 ] •••■'  p 7 ) 2 + ( D C p 6 , 2 ] •••" p 7 ) 2 } p 8 
204:  p6+l)p6 

205:  ne x t V 
206:  ret  p7#\p8 
*1745 
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0 : "Building  C h ar  ac  t e r 1 z at  i o n P r o g r am  s " : 

1 : "EM  T h e o r y d e v e 1 o p € d b y B r . M . T . M a " : 

2:  "HPL  Programs  supervised  by  M.  G.  Arthur": 

3 : "National  Bur e au  o f S t an d ar  d s " : 

4 : "Boulder,  Col o r ad o " : 

5:  "FIELD  Program:  File  3": 

6:  "FLAG  1 CONTROLS  PRINTING  OF  CONSTANTS ": sfg  1 

? : "FLAG  3 CONTROLS  PRINTING  OF  MA  AND  NR": sfg  3 

8:  "FLAG  7 CONTROLS  PRINTING  OF  CURRENTS ": sfg  7 


9: 

r ad 

10: 

2 4 > M > N 

1 1 : 

dim  A [ M , N , 2 ] , B C M , 2 ] 

, D [ N , 2 ] , J [ N ] , I C 30] , C [ 50 , 2 ] , R C 60 ] , 

X [ H , 2 ] 

12: 

dim  W S C 1 ] 

13: 

spc  5 

14: 

prt  "FIELD  PROGRAM" 

; spc  2 

15: 

ldf  4,RC*],X[*] 

16: 

"ENTER" : 

17: 

f xd  4 

1 8 : 

2 4 > M > N 

19: 

prt  "FIELD  POINT 

POSITION: " 

20 : 

spc 

21  : 

beep; ent  "ENTER  X. 

X = ?",X 

o • 

beep; ent  "ENTER  Y. 

Y = ? " , Y 

2 3 : 

b e e p ; ent  "ENTER  Z . 

"7  — 11  "7 

24  : 

"PRINT  ENTRIES": 

OC  . 
C-  -J  • 

prt  " X= " , X , " Y=  " , Y , " 

Z = " , Z 

26: 

spc  2 

27:  dsp  "ARE  THESE  NUMBERS  CORRECT?  < Y--N  > " 
28:  ent  " " , W$ 

28 : if  W* [ 1 ] = " Y " or  [ 1 ] = " y " ; gt o +5 

30:  beep; ent  "RE-ENTER  X.  X = ?",X 

31:  beep; ent  " R E - E N TER  V . V = ? " , V 

32:  beep; ent  "RE-ENTER  Z.  Z = ?",Z 

33:  gt  o "PRINT  ENTRIES" 

34:  dsp  "CALCULATION  IN  PROGRESS" 

35:  RC38D-RC31 ] >R[S] 

36:  R[39]-RC34]>RC7] 

37:  RC40]-RC37] >R[8] 

38:  R[6]*RC7]  > F:  [ 9 ] 

39:  R[63*R[8]>R[10] 

40:  R [ 7 D *R C 8 ] > R [ 1 1 ] 


41:  1 . 481e8>R[ 14] 

42:  4*atn<n>P 
43:  2*P*R C 1 4 ] > R C 1 8 ] 

44:  4e-7*P>R[ 17] 

45:  le-9/<36*F'>  >R[  12] 

46:  1 > I [ 1 1 ] 

47:  2> I C 14] > I C 17] > I [ 10] 

48 : R [ 1 8 ] *\ < R C 1 7 ] *R [ 1 2 ] > > R L 1 6 ] 

49 : 0 > C C 22 , 1 ] ; R C 1 8 ] *R  C 1 7 ] ✓ < 4*P > > C C 22 , 2 ] 

50:  0 > C C 2 4 , 1 ] ; - . 25/CRC 1 8 ] *P*R L 1 2 ] > >CC24, 2] 
51:ifflgl;gto+6 

52:  prt  "THE  VARIABLES  OMEGA , MU , EPS  I LONKAY 
53:  prt  "ARE  LISTED  BELOW";flt  6 
54:  s p c 

55:  prt  R C 1 8 ] , R C 1 7 ] , R [ 1 2 ] , R C 1 6 ] , C C 22 , 1 ] , C C 22 
56:  s p c 


0 M S 4 P , S 0 4 P E 


2 ] , C C 2 4 , 1 ] , C C 2 4 , 
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cj  “7 

-1  } 

IJ 

58 

f o r 

N=  1 to  I [173 

5 9 

< H - 

1 

) R [ 7 3 + R [ 3 4 3 > R t 3 2 3 

60 

f o r 

M = 1 to  I C14] 

6 1 

< M - 

1 

> * R C 6 3 + R [ 3 1 3 > R C 2 9 3 

62 

U + 2 

> 

IJ 

63 

i f 

f 

1 g 7 ; g t o +7 

6 4 

s P c 

65 

prt 

"CURRENT  IN  THE  ZP=0  PLANE 

AT" 

; f xd  4 

6 6 

prt 

" X R = " , RC29  3 , " Y P = " , RC323 ; f 1 t 

6 

h>  i 
6 8 

pft 

".J  SUE  X"  , XCU,  1 3 , XCU,  2 3 , " J 

SUE 

Y" , X[U+1 

1 ^ , 

X [ U + 1 , 2 3 

69 

spc 

70 

n e x 

t. 

M 

71 

n e-  x 

t 

H 

7 

f o r 

L = I [ 1 1 3 to  I [ 1 0 3 

73 

<L- 

1 

* R.  [ 8 3 + R [ 3 7 3 > R [ 3 6 3 

74 

f o r 

N= 1 to  I [173 

*7  cr 
i ■_> 

< N - 

1 

> * R [ 7 3 + R [343  >R[323 

7 b 

IJ  + 2 

> 

IJ 

“7  “7 

i f 

f 

1 g 7 ; g t o + 7 

»*'  o 

spc 

79 

prt 

"CURRENT  IN  THE  XR=0  PLANE 

AT" 

; f xd  4 

8 0 

prt. 

" YF'  = " , R [ 32  3 , " ZP=  " , R [ 36  3 ; f 1 t 

£ 

81 

spc 

8 2 

prt 

"J  SUB  Y " , X [ IJ  , 13,  X [ IJ  , 2 3 , " J 

SUE 

Z"  , X[IJ+1 

1 3 , 

X[IJ+1 , 2 3 

'I'  "J 

spc 

8 4 

n e-  x 

t 

N 

85 

n i x 

t. 

L 

8 1 

f o r 

L = I [ 1 1 3 to  I [ 1 0 3 

‘d>  7 

(L- 

1 

> * R [ 8 3 + R [ 3 7 3 > R [363 

8 8 

f or 

M=1  to  I [14  3 

g 9 

< M - 

1 

) * F:  [ 6 3 + R [ 3 1 3 > R [ 2 9 3 

90 

U + 2 

> 

u 

91 

Q 

i f 

f 

1 g 7 ; g t o + 7 

9 3 

-•  P 
pr  t 

"CURRENT  IN  THE  YP=0  PLANE 

AT" 

;fxd  4 

94 

prt 

" XF'=  " , R [ 29  3 , " ZP=  " , R [ 36  3 ; f 1 t 

f. 

95 

spc 

96 

prt 

"J  SUE  X " , X [ IJ  , 13,  X [ IJ  , 2 3 , " J 

SUE 

Z"  , X[IJ+1 

1 ] , 

X[IJ+1 , 2 3 

9 7 

s p c 

98 

ri  e x 

t 

M 

99 

n s-  x 

t 

L 

100:  1 > F 

101:  0 > E [ F , 1 ] > E [ F , 2 ] 

102:  0 > £ [ F + 1 , 1 ] > E [ F + 1 , 2] 

103:  0 > E [ F + 2 , 1 ] > E [ F + 2 , 2] 

104:  - 1 > U 

105:  for  N=1  to  I C 1 7 ] 

106:  < N- 1 > *R [ 7 ] +R  C 34 ] > R C 32 ] 

107:  for  M=1  to  IC143 
108:  < M-  1 > * R C 6 ] +F:  C 3 1 ] } R [29  3 

1 0 9 : IJ  + 2 > U 

110:  < < X -RC293  > 2 + < Y -RE323  > 2 + Z 2 ) JRC243 

111:  c os  < R [ 1 6 ] *R C 24  3 ) /R [ 24  3 > C C 20 , 1 3 ; -s i n < R C 1 6 3 *R [ 24  3 > /R C 24  3 > C [ 20 , 2 3 
112:  <<RC163*RC243>  '-2-3  > /R  C243  --4  > CC25,  13;  -3*R  [ 1 6 3 *R  C 24  3 ■■■  RC243  -4  > C [ 25 , 2 3 

113:  (.  X-R  C 29  3 > *C  C 25 , 1 3 > C C 2 , 1 3 ; C X-R  C 29  3 > *C  C 25 , 2 3 > C C 2 , 2 3 
114:  < Y-R [ 32 3 > *C C 25 ’ 1 3 > C C 5 , 1 3 ; ( Y-R C 32 3 >*C t 25 , 2 3 > C C 5 \ 2 3 

115:  Z*C C 25 , 1 3 > C [ 1 4 , 1 3 ; Z*C C 25 , 2 3 > C C 1 4 . 2 3 
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116 
1 17 
1 18 
119 

128 

121 

122 

123 

124 

125 

126 

127 

128 
129 
1 38 

131 

132 

133 

134 

135 

136 

137 

138 

139 
1 4 0 

141 

142 

143 

144 

145 

146 

147 


160 

161 

162 

163 

164 

165 

166 
167 
163 

169 

170 

171 

172 


c 1 1 
C C 7 , 
1 ■ R C 
RC  16 
c 1 1 
c 1 1 
C C 2 8 
CC  7, 
-CC2 
c 1 1 
c 1 1 
CC9, 
e 1 1 
c 1 1 
BCF, 
BCF, 
-CC2 
c 1 1 
c 1 1 
CC7, 
c 1 1 
CC9, 
1 - • R C 
RC  16 
c 1 1 
c 1 1 
C C 28 


OMUL  - C-CC22,  1 ] , 
1 3 *R  C 9 3 > C C 7 , 1 ] ; C 
2 4 3 2 + ( X - R C 2 9 3 > * 
] * R C 2 4 ] ••••'  R [ 2 4 ] 2 + 
OMUL  " (CC 24,  1 3 , C 
'CMUL'  < C C 2 8 , 13, C 
, 1 3*RC93  >CC28,  1 3 
1 3-CC28,  1 3 >CC7,  1 
4,  1 3*<Y-RC323  > >C 
'CMUL  (CC9, 1 ] , CC 
'CMUL' ( C C 9 , 13, CC 
1 3*RC93 >CC9,  1 3 ; C 
CMUL  " (CC7,  1 3 , CC 
CMUL' (CC9, 13 , CC 
1 3 +C  C 2tf , 13+CC30, 
23+CC28, 2 3+C C 30 , 
4, 1 3*<X-RC293 ) >C 
'CMUL' CCC7,  1 3 , C C 
'CMUL' (CC7, 1 ] , CC 


C 2 8 , 1 ] , CC23, 2 
C 2 8 , 1 3 , C C 2 8 , 2 


> C C 9 , 2 ] 

3 , c C 9 , 2 3 ) 

, 1 3 , C C 9 , 2 3 > 


] > 


3 > 
] > 


-c  C 22 , 2 ] , c C 20 , 1 3 , C C 28 , 2 3 , C C 7 , 13 , C C 7 , 

C 7 , 2 ] * R C 9 3 > C C 7 , 2 ] 

C C 2 , 1 3 > C C 28 , 1 3 
<jf-RC293  )*CC2, 23  > C C28, 23 
^ ^ ^ j ^ 3 , L C 2 8 , 1 3 , U C 2 8 , 2 3 , C 
C 28 , 2 3 , C C 20 , 1 3 , C C 20 , 2 3 , C 
; L [28, 2 3 * R C 9 3 } C [28, 2 3 
3 ; CC7, 2 3-CC28, 23  > C C 7 , 23 
C 9 , 1 3 s -C  C 24 , 2 3 * < Y-R C 32  3 > 

9 , 2 3 j L C 2 , 1 3 , C C 2 , 2 3 , I J C 9 , 1 
3 i i ] j t C 2 U j 1 3 j L C 2 y j 2 3 ( U C 9 
C 9 , 2 3 * R C 9 3 > C C 9 , 2 3 

7.23,  X C IJ  , 13,  X C U , 2 3 , C C 28  , 13,  C C 28 , 2 3 > 

9 . 2 3 , X C IJ  + 1 , 13,  X C U + 1 , 2 3 , C C 3 0 , 1 3 , C C 3 0 , 2 

1 3 > B C F , 1 3 

2 3 > B C F , 2 3 

C 7 , 1 3 ; -C  C 24 , 2 3 * < X-R  C 29  3 > } C C 7 , 2 3 

7 . 2 3 , C C 5 , 1 3 , C C 5 , 2 3 , C C 7 , 1 3 , C C 7 , 2 3 ) 

7 . 2 3 ,  C C 20 , 1 3 , C C 20 , 2 3 , C C 7 , 1 3 , C C 7 , 2 3 > 

C 7 , 2 3 * R C 9 3 > C C 7 , 2 3 

-C  C 22 , 2 3 , C C 20 , 1 3 , C C 20 , 2 3 , C C 9 , 1 3 , C C 9 , 2 3 ) 


3 ) 


1 3*RC93  >CC7,  1 3 ; C 
'CMUL' (-CC22, 13, 

1 3 *R  C 9 3 > C C 9 , 1 3 ; C C 9 , 2 3 *R C 9 3 > C C 9 , 2 3 
24  3 ~2+ ( Y-R  C32  3 ) *C C 5 , 1 3 } C C 28 , 13 
3 * R C 2 4 3 -' R C 2 4 3 ~2+  < Y-R  C 32  3 > *C  C 5 , 2 3 > C C 28 , 2 3 
' CMUL ' < C C 24 , 1 3 , C C 24 , 2 3 , C [28,13, C C 28 , 2 3 , C C 28 , 1 3 , C C 28 , 
CMUL ' < C C 28 , 1 3 , C C 28 , 2 3 , C C 20 , 1 3 , C C 20 , 2 3 , C C 28 , 13, C C 28 , 
, 1 3 *R  C 9 3 > C C 28 , 1 3 ; C C 28 , 2 3 *R C 9 3 >C  C 28 , 2 3 


C C 9 , 1 3 -C  C 28 , 1 3 > C C 9 , 1 3 ; C C 9 , 2 3 -C  C 28 , 2 3 > C C 9 , 2 3 
ell  CMUL  ' ( C C 7 , 1 3 , C C 7 , 2 3 , X DU,  13, X C U , 2 3 , C C 28 , 
ell  ' CMUL ' ( C C 9 , 1 3 , C C 9 , 2 3 , X C U+ 1 , 1 3 , X C U+ 1 , 2 3 , C 
B CF+ 1 , 1 3 +C  C 28 , 1 3 +C  C 30 , 1 3 > B C F + 1 , 1 3 
B C F+ 1 , 2 3 +C  C 28 , 2 3 +C  C 30 , 2 3 > B C F+ 1 , 2 3 


1 3 , CC28, 23  ) 
C30,  13, C C 3 0 , 


3 > 

: 3 > 


23  ) 


148 

- 

c 

C 2 

4 

»— »■ 
l l 

* 

-R 

c 

ill 

9 

3 ) 

} 

c 

C 7 

y 

1 

5 -CC2 

4 

y 

cl 

3 

* 

f’ 

-R 

C 2 

9 

3 

> > 

c 

C 7 

, 

3 

149 

c 

1 

1 

' i 

3MUL' ( 

CC 

1* 

n 

i 

c 

c 

•’  y 

cL 

3 

C C 1 4 , 

1 

3 

y 

c 

c 

14 

y *=:- 

3 , 

C 

C 

i'-  y 

1 

3 , 

c 

c 

C- 

3 > 

150 

c 

1 

1 

i 

3 M IJ  L < 

CC 

"7 

y 

i 

3 , 

c 

c 

f 9 

cL 

3 

C C20, 

1 

3 

y 

c 

c 

20 

q Cl 

3, 

C 

C 

«■’  y 

1 

3 , 

c 

[ 

•*’ 

3 ) 

151 

c 

c 

•’  y 

1 

] * R C 9 3 

> c 

C 

—> 

i 

y 

1 3 

; 

c 

C 7 

2 3 * R C 9 3 

3 

C 

c 

r 

y 

2 3 

152 

- 

c 

C 2 

4 

1 3 * ( Y 

-R 

c 

3 

cL 

3 ) 

> 

c 

C 9 

y 

1 3 ; -CC2 

4 

y 

ii! 

3 

* 

T 

-R 

C 3 

2 

3 

;i  y 

c 

C 9 

, 

ii! 

3 

153 

1 

1 

'i 

:mul ' ( 

CC 

9 

y 

1 

3 , 

c 

c 

9, 

3, 

C C 1 4 , 

1 

3 

y 

c 

c 

14 

q lil 

3, 

c 

C 

9, 

1 

3 , 

c 

c 

9 

3 ) 

154 

c 

1 

1 

3 M U L ' < 

CC 

9 

y 

1 

3 , 

c 

[ 

9 , 

cL 

3 

CC  20, 

1 

3 

y 

c 

c 

2 0 

y ^ 

3 , 

c 

C 

9, 

1 

3, 

c 

c 

9 

•Si 

3 ) 

155 

c 

c 

Q 

- 9 

1 

] *R  C 9 3 

> c 

[ 

9 

» 

1 3 

5 

c 

C 9 

y 

2 3 * R C 9 3 

> 

c 

c 

9 

2 3 

156 

c 

1 

l 

x i 

3 M U L ' < 

CC 

y 

1 

3, 

c 

c 

i , 

cL 

3 

X C IJ  , 1 

3 

y 

X 

C 

u 

y ^ 

3 , 

CC 

2 

a 

, 1 

3 

, C 

c 

2 

3 ) 

157 

c 

1 

i 

' i 

3MUL ' ( 

CC 

9 

y 

1 

3, 

c 

c 

9, 

cl 

3 

XCU+1 

y 

l 

3 

y 

cu 

+ 1 

y ^ 

3 

y 

CC 

3 

0, 

1 

3 

, ' 

:c 

30 

158 

B 

c 

F + 

l l 

+ 
l.  .1 

r — i 

«i!  o 

9 

1 

3 

+ c 

c 

3 

0. 

1 

3 

•BCF+2 

y 

l 

3 

159 

B 

[ 

F + 

2 

2 3+CC 

•_i  r, 
c L O 

9 

cL 

3 

+ c 

c 

0, 

2 

3 

BCF  + 2 

9 

cl 

3 

3 > 


n e x t M 
n e-  x t N 

for  L = I C 1 1 3 to  I C 1 0 3 
(L-l )*RC83+RC373 >RC363 
for  N=1  to  I C 1 73 
(N-l > * R C 7 3 + R C 3 4 3 >RC32  3 
U + 2 > IJ 

••  ( X A 2 + ( Y - R C 3 2 3 >'  • 
c o s ( R C 1 6 3 * F:  C 2 5 3 ) 

( ( R C 1 6 3 * R C 2 5 3 ) • 2 
X*CC26,  1 3 >CC3,  1 3 
Y - RC323  )*CC26,  1 
(2-RC363 )*CC26, 1 


2 + < Z - R C 3 6 3 ) -•  2 > > R C 2 5 3 

/ R [ 25  3 > C C 2 1 , 1 3 ; -si n < RC  163*RC 
- 3 ) / R [ 25  3 4 > CC26,  1 3 ; - 3 * R C 1 6 3 
; X * C C 2 6 , 2 3 > U C 3 , 2 3 

3 > C C6,13;  ( Y-R C 32 3 >*CC26, 23  >C 
3 > C C 1 5 , 1 3 ; ( Z - R C 3 6 3 > * C C 2 6 ,23  ) 


25  3 > RC253  C 
*RC253  RC253 

C 6 , 2 3 
C C 1 5 , 2 3 


21.23 
4 > C C 2 6 
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173 

174 

175 

176 

177 

1 79 
180 
181 
182 

183 

184 

185 

186 

187 

188 
1 89 
198 

191 

192 

193 

194 

195 
1 9 6 

1 97 

198 

199 

2 0 0 
201 
2 0 2 
2 0 3 

20  4 
2 0 5 
206 
207 
2 0 8 
2 0 9 
210 

21  1 

214 

215 

216 
217 
2 1 8 

219 

220 
221 
dd  2 
2 2 3 


229 

230 


-C [ 24 , 1 3 * < Y-R  C 32  3 > > C C 1 0 , 1 ] ; -C  C 24 , 2 ] * < Y-R C 32 3 > > C C 1 8 , 2 3 

c 1 1 CMUL  ■"  < C C 1 0 , 1 3 , C C 1 0 , 2 3 , C C 3 , 1 3 , Ct  3 , 2 3 , C C 1 0 , 1 3 , C C 1 0 , 2 3 > 

c 1 1 CMUL  < C [ 1 0 , 1 3 , C [ 1 0 , 2 3 , C C 2 1 , 1 3 , C C 2 1 , 2 3 , C C 1 0 , 1 3 , C [ 1 0 , 2 3 > 

C C 1 0 , 1 3 *R [ 1 1 3 > C C 1 0 , 1 3 ; C C 1 0 , 2 3 *R C 1 1 3 > C C 1 0 , 2 3 

-C  C 24 , 1 3 * < Z-R  C 36  3 > > C [ 1 1 , 1 3 ; -C  C 24 , 2 3 * < Z-R [ 36 3 > > C [ 1 1 , 2 3 

c 1 1 CMUL  < CCll,  1 3 , C [ 1 1 , 2 3 , C [ 3 , 1 3 , C C 3 , 2 3 , C C 1 1 , 1 3 , C C 1 1 , 2 3 ) 

c 1 1 •-  CMUL  < C [ 1 1 , 1 3 , C [ 1 1 , 2 3 , C [ 2 1 , 1 3 , C L 2 1 , 2 3 , C C 1 1 , 1 3 , C C 1 1 , 2 3 > 

C [ 1 1 , 1 3 *R C 1 1 3 > C C 1 1 , 1 3 ; C C 1 1 , 2 3 *R C 1 1 3 > C C 1 1 , 2 3 

c 1 1 ■-  CMUL  < C [ 1 0 , 1 3 , C E 1 0 , 2 3 , X C U , 1 3 , X C IJ  , 2 3 , C C 28  , 1 3 , C [ 28 , 2 3 > 


r 

1 

1 

CMUL  - 

C[ 

1 

1 

9 

1 

3 , 

c 

[ 

1 

1 

, 

23  , 

[IJ+1 

9 

1 

3 

9 

[ 

IJ  + 

1 

9 

2 

3 

9 

c 

[ : 

:0 

9 

1 

3 

9 

c 

[ 

50 

9 •:4 

3 

E 

[ 

F 

, 1 

3 + C [ 2 

o 

, l 

3 

+ 

c 

[ 

3 0 

, 

1 

3 

/ 

B 

[F, 

1 

3 

E 

[ 

F 

9 ^ 

3 + C [ 2 

y 

, 4‘ 

3 

+ 

c 

[ 

30 

, 

3 

3 

> 

B 

CF, 

3 

r 

1 

1 

CMUL" 

-c 

[ 

2 

JI 

1 3 

, 

- 

c 

[ 
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2,23 

, C [ 2 

1 

9 

1 

3 

9 

c 

[2 

1 

9 

cl 

3 

9 

c 

[ 

0 

, 

1 

3 

9 

c 

[ 

10 

9 ^ 

3 

c 

[ 

1 

0, 

1 3*R[ 

1 

1 3 

> 

f 

c 

1 

0 , 

1 

3 

9 

c 

[ 

10, 

3 *R  [ 

1 

1 

3 

/ 

C 

[ 

10 

, 

2 

] 

1 

/ 

R 

[2 

5 3 •-•2  + 

£ 

Y - 

R 

[ 

y 

cl 

3 > 

■fc 

c 

[ 

6 

, 

1 3 : 

L 

[28, 

1 

3 

R 

[ 

1 

6 3 

* R [25 

] 

/R 

[ 

•7 

5 

3 

cl 

+ 

<: 

Y 

R 

[32  3 

>*C[ 

6 

, 

•“» 

C. 

3 

/ 

c 

[2 

O 

9 

cl 

3 

c 

1 

1 

CMUL' 

»; 

C[ 

cl 

£ 

ji 

1 

3 , 

c 

[ 

cl 

£ 

9 

23  , 

c 

[24, 

1 

3 

9 

c 

[ 

LI 

4, 

'cl 

3 

9 

c 

[ 

'ci 

1 

] 

, 

c 

[ 

ci 

O 

9 *=- 

3 > 

c 

1 

1 

CMUL' 

C[ 

cl 

'r 

ji 

1 

3 , 

c 

[ 

2 

£ 

9 

23, 

c 

[21, 

1 

3 

9 

c 

[ 

2 

1 , 

cl 

3 

9 

c 

[ 

L' 

1 

3 

9 

c 

[ 

L! 

8 

9 ^ 

3 > 

c 

[ 

cl 

8 , 

1 3*R[ 

1 

1 3 

> 

c 

[ 

2 

8 , 

1 

3 

9 

c 

[ 

cl  O ^ 

3 * R [ 

1 

1 

3 

c 

[ 

d y 

9 

2 

] 

c 

[ 

1 

0, 

1 3 - C [ 

2 

9 

1 

3 

\ 

c 

[ 1 

0 

9 

1 

3 

9 

C[  1 

0 

,23- 

C 

[ 

d 

O 

9 
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c 

[ 

1 

0 

9 

ci 
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- 
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[ 

24 

.13*': 

7 

-R 

[ 

•“ 

b 

3 

|l  ’[• 
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[ 

1 

1 

9 
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C [24 

9 

'cl 

3 

* 

z 

-R 

c 

• j 

b 

3 

:> 

> 

C[  1 

1 

9 

d 

3 

c 

1 

1 

CMUL' 

( 

C[ 

1 

1 

ji 

1 

3 , 

c 

[ 

1 

1 

9 

23  , 

f 

[6,1 

3 

9 

C 
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6 

, 

23 

9 

c 

[ 

1 

l 

, 

13, 

c 
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314 

ent  11  ",W$ 

3 1 5 

if  W $ [ 1 3 = " V 

" or  W $ [ 1 3 = 

" y 

" ; g t o “ENTER" 

3 1 6 

beep; dsp  " C 

HLCULHTION 

COMPLETED" 

317 

end 

3 1 8 

" ATAN2 " : 

319 

■at.  n < ib  s < p 1 > 

/.ib  S-  < p2  > > > p 

ij;! 

320 

if  p2<0  arid 

pi >0; { - p 3 > 

P 3 

321 

it"  p2<0  ind 

p 1 < 0 ; p 3 - •(  > 

p 3 

322 

if  p2>0  ind 

p 1 < 0 ; p 3 * - 1 

>P 

>! 

■->  cL  • J 

r e t f:'  3 

32  4 

"CMIJL  " : 

3 2 5 

p 1 p 3 “ p 2 p 4 > p 

7 ; p 1 p 4 + p 2 p 3 

> P 

8; p 7 > p 5 ; p8 >p6 

326 

" C D I V " : 

3 2 7 

p 3 p 3 + p 4 p 4 > p 

7 ; ( p 1 p 3 + p 2 p 

4::' 

P 7 > p 1 0 

328 

p 2 p 3 — p 1 p 4 ) 

p 7 > p 1 1 ; p 1 0 

> P 

5 ; p 1 1 }■  p 6 ; r e t 

=" , R[33 , " P HSZ= " , RC21 3 


r e t 


*2682 1 
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B.2  Variable  Assignments 


Complex  Variable  Assignments 

The  complex  variables  are  stored  in  arrays  with  one  extra  dimension. 

Thus  the  complex  variable  can  be  stored  in  two  parts:  the  real  part,  stored 

in  the  top  (X[X,1]),  and  the  imaginary  part,  stored  in  the  bottom  (X[X,2]). 

All  complex  arrays  are  assigned  their  own  array  variable  while  the  rest  of  the 
complex  variables  are  assigned  in  the  array  variable  "C". 


FORTRAN  VARIABLE 
A ( M , N ) 

B(M) 

OIAG(N) 

XS(N) 


HPL  VARIABLE 
A[M,N,2] 


B[M,2] 
D[N ,2] 
X[N,2] 


C[50,2]pool  of  complex  variables 


AML3 

AMN1 

ANL2 


BML3 


CJX1 

CJX3 


BMN1 


BNL2 


CJZ3 


CJV1 


CJY2 

CJZ2 


CMN1 

CNL2 


Cl 

C2 


CML3 


ESUBZ 


C[l] 
C[  2] 

C [3] 
C[  4] 
C[5] 
C[  6] 

C [ 7] 
C[  8] 
C[  9] 
C[10] 
C[11] 
C[12] 
C[  1 3] 
C[  14] 
C[15] 
C[16] 
C[  1 7] 
C[18] 
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GML3 

C[  19] 

GMN1 

C[20] 

GNL2 

C[21] 

0MS4P 

C[22] 

RSUBV 

C[23] 

S04PE 

C[24] 

SI 

C[  25] 

S2 

C[26] 

S3 

C[  27] 

TEMP 

C[28] 

TEMP 

C[29] 

TEMP 

C[30] 
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Real  Variable  Assignments 

The  most  commonly  used  real  variables  are  assigned  simple  variables  and 
the  rest  of  the  real  variables  are  assigned  in  the  variable  array  "R". 


FORTRAN  VARIABLE 
PI 
R 

THETA 

X 

Y 

Z 

AMPX 

AMPY 

AMPZ 

ARG 

DTOR 

DX 

DY 

DZ 

DXDY 

DXDZ 

DYDZ 

EPS 

ESUBZ 

FREQ 

FTOM 

KAY 

MU 

OMEGA 

PHSX 

PHSY 

PHSZ 

RESID 

RML3 

RMN1 


HPL  VARIABLE 
P 
R 
T 
X 
Y 
Z 

R[60]pool 
R[l] 

R[2] 

R[3] 

R[4] 

R[5] 

R[6] 

R[7] 

R[8] 

R[9] 

R[10] 

R[  1 1] 

R[12] 

R[  13] 

R[  14] 

R[15] 

R[16] 
R[17] 
R[18] 
R[19] 
R[20] 
R[21] 

R[  22] 
R[23] 
R[24] 


real  variables 
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RNL2 

RSUBT 

TC 

XP 

XPP 

X START 

YP 

YPP 

YSTART 

ZD 

ZP 

ZSTART 

X2 

Y2 

Z2 


RC  25] 
R[26] 
R[28] 
R[29] 
R[  30] 
R[31] 
R[32] 
R[33] 
R[34] 
R[  35] 
R[36] 
R[37] 
R[  38] 
R[39] 
R[  40] 
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Integer  Variable  Assignments 

The  most  commonly  used  integer  variables,  and  variables  used  in  "for" 
loops,  are  assigned  simple  variables.  The  integer  array  "JPVT"  is  stored  in 
its  own  variable  array  "J"  while  the  rest  of  the  integer  variables  are  stored 
in  the  array  "I". 


FORTRAN  VARIABLE 

I,  IX 

J, IY 

K , IZ 
L 

M 

N 

IF 

JU 

JPVT 

II 

IX 

IY 

IYF 

IZ 

KX 

KY 

KYF 

KZ 

LF 

LS 

MA 

MDIMA 

MF 

MS 

NA 

NF 

NS 


HPL  VARIABLE 
I 
J 
K 
L 

M(also  used  to  dimension  arrays) 
N(also  used  to  dimension  arrays) 
F 
U 

J[N] 

I(30)pool  of  integer  variables 

I[l] 

I C 2 ] 

I [3] 

I[4] 

I[5] 

I[6] 

I[7] 

I [8] 

I C 9] 

i[io] 

i[n] 

I [ 1 2] 

I[13] 

I C 14] 

I C 1 5] 

I C 16] 

I[17] 

I [18] 
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B.3  Instructions  for  Running  HPL  Programs 


To  run  the  HPL  CURRENT  program,  follow  this  procedure: 

1.  With  the  power  to  the  computer  turned  on,  insert  the  magnetic  program 
tape  in  accordance  with  the  computer's  operating  instructions. 

2.  Press  LOAD  and  type  in  1,  or  type  ldf  1,  and  press  EXECUTE.  This  will 
load  CURRENT  program.  Part  1,  into  the  computer. 

3.  Make  any  program  changes,  such  as  flags,  as  needed  or  desired.  These  are 
discussed  in  section  6.3.  The  procedure  listed  here  is  based  on  no 
changes  to  the  program  as  shown  in  section  B.l. 

4.  Press  RUN.  This  will  cause  the  computer  to  print  the  heading  "CURRENT 
PROGRAM"  on  the  paper  tape.  Next,  the  computer  will  print  "COORDINATES 
OF  FIELD  POINTS:"  on  the  tape. 

5.  The  computer  will  beep  and  display  "ENTER  XI.  XI  = ?"  as  a prompt. 

6.  Type  in  the  first  X-coordinate  value  and  press  CONTINUE.  The  computer 

will  beep  and  display  "ENTER  X2.  X2  = ?." 

7.  Type  in  the  second  X-coordinate  value  and  press  CONTINUE. 

8.  Continue  in  this  manner  until  XI,  X2,  Yl,  Y2,  Zl,  and  Z2  have  all  been 

entered. 

9.  After  Z2  has  been  entered,  the  computer  will  print  out  the  values  of  XI 
through  Z2  that  were  entered.  Then  it  will  display  "ARE  THESE  NUMBERS 
CORRECT?  (Y/N)." 

10.  If  all  are  correct,  type  in  either  y or  Y and  press  CONTINUE.  If  any  are 
not  correct,  type  in  either  n or  N and  press  CONTINUE. 
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11.  If  you  have  typed  n or  N,  the  computer  will  beep  and  display  "RE-ENTER 
XI.  XI?"  as  a prompt.  Repeat  steps  6 through  10  until  the  numbers 
printed  out  in  step  9 are  correct. 

12.  When  you  type  y or  V,  the  computer  will  display  "THIS  CALCULATION  TAKES  3 
MINUTES."  Take  a break  while  the  computer  does  its  thing. 

13.  When  the  computer  completes  this  calculation  of  the  A coefficients  and 
has  performed  the  QR  factorization  on  the  [A]  matrix,  it  automatically 
loads  CURRENT  program.  Part  2 from  file  #2. 

14.  The  computer  then  prints  "EM  FIELD  DATA  AT  FIELD  POINTS:",  beeps,  and 
displays  "ENTER  AMPX.  AMPZ  = ?"  as  a prompt. 

15.  Type  in  the  amplitude  of  the  x-direction  component  of  the  measured  field 
at  data  point  (XI,  Y1 , Zl)  and  press  CONTINUE.  The  computer  will  beep 
and  display  "ENTER  PHSX.  PHSX  = ?"  as  a prompt. 

16.  Type  in  the  phase  of  the  X-direction  component  of  the  measured  field  at 
data  point  (XI,  Yl,  Zl)  and  press  CONTINUE.  The  computer  will  beep  and 
display  "ENTER  AMPY.  AMPY  = ?"  as  a prompt. 

17  Continue  in  this  manner  until  AMPX,  PHSX,  AMPY,  PHSY,  AMPZ,  AND  PHSZ  for 
data  point  (XI,  Yl , Zl)  have  been  entered. 

18.  After  PHSZ  has  been  entered,  the  computer  will  print  out  the  values  of 
AMPX  through  PHSZ  that  were  entered.  Then  it  will  display  "ARE  THESE 
NUMBERS  CORRECT?  (Y/N)." 

19.  If  correct,  type  in  either  y or  Y and  press  CONTINUE.  If  not  ^or,o^., 
type  in  either  n or  N and  press  CONTINUE. 

20.  If  you  have  typed  n or  N,  the  computer  will  beep  and  display  R 

AMPX.  AMPX  = ?"  as  a prompt.  Repeat  steps  15  through  19  until  tue 
numbers  printed  out  in  step  18  are  correct. 
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21.  When  you  type  y or  Y,  the  computer  again  prints,  "EM  FIELD  DATA  AT  FIELD 
POINTS:",  beeps,  and  displays  "ENTER  AMPX.  AMPX  = ?"  as  a prompt. 

22.  Repeat  steps  15  through  21  for  measured  field  data  at  data  points  (X2, 

Yl,  Zl),  (XI,  Y2 , Zl),  ( X2,  Y2,  Zl),  (XI,  Yl,  Z2),  (X2,  Yl,  Z2),  (XI,  Y2 
Z2),  and  (X2,Y2,Z2)  in  that  order. 

23.  After  the  measured  field  data  have  all  been  entered  correctly  for  data 
point  (X2,  Y2,  Z2)  and  you  have  pressed  y or  Y,  the  computer  will  display 
"CALCULATION  IN  PROGRESS."  This  calculation  requires  several  seconds. 

24.  When  the  computer  completes  this  calculation  of  the  currents  J,  it 
automatically  loads  and  runs  the  FIELD  program  stored  in  file  #3. 

To  run  the  HPL  FIELD  program,  follow  this  procedure: 

1.  If  the  FIELD  program  is  to  be  run  immediately  following  the  CURRENT 
program,  the  computer  will  have  printed  the  heading  "FIELD  PROGRAM"  on 
the  paper  tape  immediately  after  loading  the  FIELD  program  in  step  24, 
above.  Procede  to  step  4,  below. 

If  the  field  program  is  to  be  run  at  some  later  time  following  the 
running  of  the  CURRENT  program,  and  the  computer  memory  has  been  erased 
so  that  it  no  longer  contains  the  results  of  running  the  CURRENT  program, 
begin  with  step  2,  below. 

2.  With  the  power  to  the  computer  turned  on,  insert  the  magnetic  program 
tape  in  accordance  with  the  computer's  operating  instructions. 

3.  Press  LOAD  and  type  3,  or  type  Idf  3,  and  press  EXECUTE.  This  will  load 
the  FIELD  program  into  the  computer. 

4.  Make  any  program  changes  as  needed  or  desired.  These  are  discussed  in 
section  6.4.  The  procedure  listed  here  is  based  on  no  changes  to  the 
program  as  shown  in  section  B.l. 


110 


5.  Press  RUN.  This  will  cause  the  computer  to  print  the  heading  "FIELD 
PROGRAM"  on  the  paper  tape.  Next,  the  computer  prints  "FIELD  POINT 
POSITION:"  on  the  tape. 

6.  The  computer  will  beep  and  display  "ENTER  X.  X = ?"  as  a prompt. 

7.  Type  in  the  X-coordinate  of  the  arbitrary  field  point  at  which  the  field 
strength  is  to  be  computed  and  press  CONTINUE.  The  computer  will  beep 
and  display  "ENTER  Y.  Y = ?." 

8.  Type  in  the  Y-coordinate  of  the  arbitrary  field  point  and  press 
CONTINUE.  The  computer  will  print  out  the  values  of  X,  Y,  and  Z that 
were  entered.  Then  it  will  display  "ARE  THESE  NUMBERS  CORRECT?  (Y/N).' 

10.  If  correct,  type  in  either  y or  Y and  press  CONTINUE.  If  not  correct, 
type  in  either  n or  N and  press  CONTINUE. 

11.  If  you  have  typed  n or  N,  the  computer  will  beep  and  display  "RE-ENTER 
X.  X = ?"  as  a prompt.  Repeat  steps  7 through  10  until  the  numbers 
printed  out  in  step  9 are  correct. 

12.  When  you  type  y or  Y,  the  computer  will  display  "CALCULATION  IN 
PROGRESS."  This  calculation  requires  several  seconds. 

13.  When  the  computer  completes  this  calculation  of  the  field  at  point  (X, 

Z),  it  prints  "E  SUB  X,  E SUB  Y,  and  E SUB  Z at  the  FIELD  POINT."  Then 
it  prints  out  the  values  of  X,  Y,  Z,  AMPX,  PHSX,  AMPY,  PHSY , AMPZ,  AN?  -v- 

14.  Then  the  computer  displays  "DO  ANOTHER  FIELD  POINT?  (Y/N)." 

15.  If  yes,  type  in  either  y or  Y and  press  CONTINUE.  If  no,  type  in  e't'V” 
n or  N and  press  continue. 

16.  If  you  have  tyDed  y or  Y,  the  computer  will  print  "FIELD  POINT  -DSI  li  i: 
on  the  tape,  and  continue  from  step  6,  above.  If  you  have  Typed 

the  computer  displays  "CALCULATION  COMPLETED"  and  ends  the  ru”. 
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Figure  1.  A theoretical  model  of  penetration  of 
electromagnetic  waves  into  a cavity 
with  apertures. 
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Figure  2.  A solid  perfectly  conducting  box  without 
apertures  in  free  space. 
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Figure  3 


Normalized  induced  currents  on  Face  I of  Figure  2 

O 

for  a (0.5X)  cube  illuminated  by  a normal  incidence, 
plane  wave  H-  toward  +y  -direction  (a)  |Jj,  and 

(»)  Uz|.  [s] 
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Figure  4 


Normalized  |J  | on  Face 
for  a (O.IX)^  cube  with 


I of  Figure  2 
normal  incidence. 
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x A 


Figure  5.  Normalized  |JZ|  on  Face  I of  Figure  2 

for  a (0. IX)3  cube  with  normal  incidence. 
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Figure  6.  A plane  wave  incident  normally  on  a conducting 
cavity  with  a slot  aperture  (a)  side  view, 
and  (b)  front  view. 
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Figure  7.  A network  representation  for  the  model  in 
Figure  6. 
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Figure  8.  Electric  field  distribution  for  TE^  mode  in  a 
copper  cavity  at  natural  resonance,  (x  * z * 0). 
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Figure  9.  Idealized  rectangular  building. 
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